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Abstract 

We study the g-state Potts model with nearest-neighbor coupling v = e^"^ — 1 in the 
limit g, u — > with the ratio w = v/q held fixed. Combinatorially, this limit gives rise to 
the generating polynomial of spanning forests; physically, it provides information about 
the Potts-model phase diagram in the neighborhood of {q, v) = (0, 0). We have studied 
this model on the square and triangular lattices, using a transfer-matrix approach at 
both real and complex values of w. For both lattices, we have computed the symbolic 
transfer matrices for cylindrical strips of widths 2 < L < 10, as well as the limiting 
curves B of partition-function zeros in the complex tn-plane. For real w, we find two 
distinct phases separated by a transition point w = wq, where wq = —1/4 (resp. 
wq = —0.1753 ± 0.0002) for the square (resp. triangular) lattice. For w > wq we find a 
non-critical disordered phase that is compatible with the predicted asymptotic freedom 
as w; ^ -|-oo. For w < wq our results are compatible with a massless Berker-Kadanoff 
phase with central charge c = — 2 and leading thermal scaling dimension xt,i = 2 
(marginally irrelevant operator). At w = wq we find a "first-order critical point": the 
first derivative of the free energy is discontinuous at wq, while the correlation length 
diverges as w | (and is infinite at w = wq). The critical behavior at w = wq seems to 
be the same for both lattices and it differs from that of the Berker-Kadanoff phase: our 
results suggest that the central charge is c = —1, the leading thermal scaling dimension 
is xt,i = 0, and the critical exponents are v = 1/d = 1/2 and a = 1. 

Key Words: Potts model, g ^ limit, Fortuin-Kasteleyn representation, spanning forest, 
transfer matrix, conformal field theory, phase transition, Berker-Kadanoff phase, square 
lattice, triangular lattice, Beraha-Kahane-Weiss theorem. 
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1 Introduction 



1.1 Phase diagram of the Potts model 

The Potts model [1-3] on a regular lattice £ is characterized by two parameters: the 
number q of Potts spin states, and the nearest-neighbor coupling v = e^"^ — 1.^ Initially q is 
a positive integer and f is a real number in the interval [—1, +oo), but the Fortuin-Kastclcyn 
representation (reviewed in Section 2.1) shows that the partition function Zc^q^v) of the q- 
state Potts model on any finite graph G is in fact a polynomial in q and v. This allows us 
to interpret q and v as taking arbitrary real or even complex values, and to study the phase 
diagram of the Potts model in the real (g, f )-plane or in complex (g, f )-space. 

According to the Yang-Lee picture of phase transitions [4] , information about the possible 
loci of phase transitions can be obtained by investigating the zeros of the partition function 
for finite subsets of the lattice C when one or more physical parameters (e.g. temperature or 
magnetic field) are allowed to take complex values; the accumulation points of these zeros 
in the infinite-volume limit constitute the phase boundaries. For the Potts model, therefore, 
by studying the zeros of Zciq-iV) in complex (g, f)-space for larger and larger pieces of the 
lattice £, we can learn about the phase diagram of the Potts model in the real (g, w)-plane 
and more generally in complex (g, v)-space. 

Since the problem of computing the phase diagram in complex (g, v)-space is difficult, 
it has proven convenient to study first certain "slices" through (g, ■;;)-space, in which one 
parameter is fixed (usually at a real value) while the remaining parameter is allowed to vary 
in the complex plane. Thus, the authors and others (notably Shrock and collaborators) 
have in previous work studied the chromatic polynomial {v — —1), which corresponds to the 
zero-temperature limit of the Potts antiferromagnet^; the flow polynomial {v = — g), which 
is dual to the chromatic polynomial [8]; the g-plane behavior for fixed real v in both the 
ferromagnetic and antiferromagnetic regions [9-15]; and the f -plane behavior for fixed real 
g, notably either g = 2, 3, 4, . . . [9-22] , g ^ 1 [14, 15] or g = [14, 15, 23] . 

In this paper we will study yet another slice, namely the limit q,v ^ with w — 
v/q finite.^ Prom a combinatorial point of view, this limit corresponds to the generating 
polynomial of spanning forests (see Section 2.2). From a physical point of view, this limit 
corresponds to investigating the behavior of the phase diagram in a small neighborhood of 

■^Here we are considering only the isotropic model, in which each nearest-neighbor edge is assigned the 
same coupling v. In a more refined analysis, one could put (for example) difi'erent couplings V\,V2 on the 
horizontal and vertical edges of the square lattice, difii'erent couplings v\,V2,vz on the three orientations of 
edges of the triangular or hexagonal lattice, etc. 

^See [5] for an extensive list of references through December 2000; and see [6, 7] for more recent work. 

■^We stress that, despite the fact that w — > 0, this is not an "infinite-temperature" limit in any relevant 
physical sense, because v and q are simultaneously varying, and one must take account of the joint efi'ect 
of the two parameters. Indeed, it turns out that for q small it is w = v/q, and not v itself, that plays the 
role of an "inverse temperature". Thus, the radius of convergence of the small-u expansion is asymptotically 
proportional to \q\ when g ^ 0; that is why the small-w expansion in the spanning-forest model is convergent 
for small \w\ but not in general for all w (see Section 3.1 and Appendix A). In particular, phase transitions 
can and do occur at finite values of real or complex w, as we shall see in this paper. 

For the two-dimensional lattices considered here, it will turn out that there is no phase transition for 
positive real w. But this expresses a deep fact about the critical behavior for these lattices, namely the 
asymptotic freedom [24]. The situation is likely to be quite different for three-dimensional lattices, for which 
there may well exist a "ferromagnetic" phase transition at some finite positive real Wcrit- 
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the point {q, v) — (0, 0) — more precisely, to investigating those phase-transition curves 
that pass through (q, v) = (0, 0) with finite slope w.^ This limit takes on additional interest 
in light of the recent discoveries [24] that (a) it can be mapped onto a fermionic theory 
containing a Gaussian term and a special four-fermion coupling, and (b) this latter theory 
is equivalent, to all orders in perturbation theory in l/w, to the A?^- vector model at = — 1 
with /3 = —w, and in particular is perturbatively asymptotically free in two dimensions, 
analogously to two-dimensional o"-models and four-dimensional nonabelian gauge theories. 

Further motivation for this study comes from our ongoing work [25] on the phase diagram 
and renormalization-group flows of the Potts model on the square and triangular lattices. 
These phase diagrams have been actively studied (see e.g. [26] for an extensive set of ref- 
erences); but certain aspects of the phase diagram in the antiferromagnetic regime remain 
unclear, notably on the triangular lattice. Let us begin, therefore, by giving a brief summary 
of what is known and what is mysterious. We shall parametrize the interval < g < 4 by 

g = 4cos^(7r/5) with 2 < 5 < oo . (1.1) 



Square lattice: Baxter [27,28] has determined the exact free energy for the square-lattice 
Potts model on two special curves in the {q, t')-plane: 

V = ±y/^ (1.2) 

V = -2 ± ^4 - ? (1.3) 

These curves are plotted in Figure 1. Curve (1.2+) is known [27] to correspond to the 
ferromagnetic phase-transition point, whose critical behavior is by now well understood 
[29-32]: for < g < 4 the critical ferromagnet is described by a conformal field theory 
(CFT) [30-32] with central charge 



and thermal scaling dimensions''^ 

2(6-1) 



-T,„ = t^ [„= 1,2.3,...]. (1.5) 



^By a standard duality transformation [see (2.21)-(2.24) below], the spanning-forest model on the lattice 
C at parameter w is equivalent to the q = model on the dual lattice C* at parameter v = l/w. In particular, 
our results concerning the spanning-forest model on the square and triangular lattices can be immediately 
translated into results for the q = Potts model at fixed v on the square and hexagonal lattices. We leave 
these straightforward translations to the reader. 

^For a scaling operator Oi, wc denote by yi its renormalization-group eigenvalue, so that j/i > (resp. = 0, 
< 0) corresponds to a relevant (resp. marginal, irrelevant) operator. Then Xi = d — yi is the corresponding 
scaling dimension, where d denotes the system's spatial dimension. When = 2 (as in this paper), we can 
use the language of conformal field theory (CFT): if a scaling operator Oi has conformal weight hi (resp. hi) 
with respect to the holomorphic (resp. antiholomorphic) variable, then Xi = hi + hi is the scaling dimension 
and Si = hi — hi is the spin. 

^See [33, Appendix A.l] for a convenient summary of the critical properties of the two-dimensional 

ferromagnetic Potts model. The parameter x used there is related to 5 by x = —2/S. Note also that there is 
a typographical error in equation (A. 10) of [33], which should read A^^s = {[2{s — r) + sx]"^ — a;^)/(8[(2 -|-a;)]. 

In CFT language, the n-th thermal scaling dimension XT,n is twice the conformal weight of the operator 
(f'n+i.i obtained by repeated fusion of the fundamental thermal operator 02.1- Please note that for integer 
(5 > 4, the central charge (1.4) coincides with that of a unitary minimal model, for which some zero coefficients 
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Baxter [28] conjectured that curve (1-3+) with < q < 4 corresponds to the antifer- 
romagnetic critical point. For q = 2 this gives the known exact value [37]; for g = 3 it 
predicts a zero-temperature critical point (f^ = —1), in accordance with strong analytical 
and numerical evidence [38-45]; and for g > 3 it predicts that the putative critical point lies 
in the unphysical region v < —1, so that the entire physical region —l<v<0 lies in the 
disordered phase, in agreement with numerical evidence ior q — A [44] . 

Saleur [46, 47] pursued the investigation of the phase diagram and critical behavior in 
the antiferromagnetic (— 1 < f < 0) and unphysical (f < — 1) regimes. Firstly, he extended 
Baxter's conjecture by suggesting [47] that the critical antiferromagnetic Potts model (1.3+) 
with < < 4 is described by a CFT with central charge 

c = 2-6/(5 . (1.6) 

Furthermore, Saleur conjectured that the leading thermal scaling dimension xt in this CFT 
is given by 

XT = 4/5 , (1.7) 
so that the associated critical exponent is 

" " 2^ = W^) ■ ^'-^^ 

This conjecture agrees with the known value at g = 2 (namely, z/ = 1); on the other hand, 
it predicts the value u — S/A at q = 3, which is incorrect [44,45]. However, after discussing 
the representation of the 3-state model as a free bosonic field with g — 1/3, Saleur finds 
another operator with xt = 3/2 and u = 2, which is the correct answer [44,45]. Clearly, this 
latter operator (and not the initially predicted one) is the leading thermal operator in the 
square-lattice 3-state Potts antiferromagnet [45].^ 

It is not clear whether (1.7)/(1.8) should be expected to be correct for q ^ 3. We defer 
testing the validity of these expressions, as a function of q, to a subsequent paper [25]. 
However, in the limit — > along the curve (1.3+) , we shall find that there are at least 

appear in the usual Coulomb-gas fusion rules. This means that the operator (f>n+is will be present in the 
Kac table only for n ^ 1,2, .... 5 — 3. In particular, for q = 2 (rcsp. q = 3), the operator with exponent 
(1.5) should be a local observable in terms of the Potts spin variables only for n = 1 (resp. n = 1,2,3). 
However, nothing would seem to prevent it from being observable in terms of the Fortuin-Kasteleyn clusters 
(see Section 2.1 below), whose definition is nonlocal in terms of the spins. Indeed, it is conceivable that the 
n = 2 operator is indeed present in the two-dimensional Ising model and causes L~^/^ corrections to scaling 
in the Fortuin-Kasteleyn bond observables [34-36] . 

''In the Coulomb gas picture [29] employed by Saleur [47] , there are two basic types of operators: electric 

(or vertex) operators of electric charge e, and magnetic (or vortex) operators of magnetic charge m. The 
thermal scaling field associated to a general operator with electric charge e G Z/3 and magnetic charge 
m e 3Z is given by 

9 9 

e am 
'^ = Tg^ — 

where g is the coupling constant of the free bosonic field onto which the original model renormalizes. In the 
square-lattice 3-state Potts antiferromagnet, g = 1/3 [47]. Saleur initially conjectured the thermal operator 
to be an electric operator of charge e = 2/3, thus leading to a;^ = 2/3 and = 3/4. The alternate conjecture 
{xt = 3/2 and v = 2) comes from identifying the thermal operator with a magnetic one with charge m = 3. 
This result agrees with the identification (made in Ref. [45]) of the thermal operator as a vortex operator 
with the smallest possible topological charge. (The normalization conventions in Refs. [47] and [45] differ: 
in Ref. [45], one has xt = -^{a^ / K + Kf}'^) with the correspondence a = ne, /3 = 2m and K = ng/2.) 
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two thermal-type operators that are more relevant than (1.7), namely one with xt — and 
another with xt ~ 0.3-0.4 (sec Sections 7.6 and 7.10 for more details). 

Saleur [47] also investigated the meaning of the other two special curves, (1.2_) and 
(1.3_). He suggested that there exists a Berker-Kadanoff phase [48] — i.e. a massless low- 
temperature phase with algebraically decaying correlation functions — extending between 
the curves (1.3±) in the range < g < 4 [i.e. throughout the hatched region in Figure 1] 
except when g is a Beraha number i?„ = 4cos^(7r/n) [n = 2,3,4, . . ., corresponding to the 
pink vertical lines in Figure 1] , and that the critical behavior of this phase is determined by an 
attractive fixed point lying on the unphysical self-dual line (1.2_). He further conjectured 
that the model on the line (1.2_) with < g < 4 — and hence throughout the Berker- 
Kadanoff phase — is described by a CFT with central charge 

provided that 6 is not an integer. Finally, he conjectured that the leading thermal scaling 
dimension in this CFT is 

= ^ . (1.10) 

Since Xt > 2, the energy is an irrelevant operator in this phase (except at g = 0, where it is 
marginal), in accordance with the fact that there is an entire interval of critical points, all 
governed by a single renormalization-group fixed point. ^ 

Finally, (1.3_) is the dual of the antiferromagnetic critical curve (1.3+). Therefore, the 
transfer matrices for (1.3-i-) with cylindrical boundary conditions are identical up to multi- 
plication by a constant. If we assume (as seems likely) that the different endgraphs needed 
in the two cases do not lead to any zero amplitudes, the theories (1.3±) should therefore be 
completely equivalent; in particular, they should have the same central charge and the same 
thermal scaling dimensions. (However, a local operator in one theory could correspond to a 
nonlocal operator in the dual theory.) This equivalence is corroborated by the fact that in 
CFT, the complete operator content is linked to the modular-invariant partition function on 
the torus [31]; obviously, in this geometry the lattice coincides with its dual. 

Let us tentatively accept these conjectures and determine their implications for the limit 
g, ^; — > with w — v/q fixed. The self-dual curves (1.2±) pass through {q,v) — (0,0) with 
slope w ~ ±oo, while the antiferromagnetic critical curve (1.3+) passes through (g, v) = (0, 0) 
with slope w = —1/4. Thus, all values —1/4 < w < +oo are predicted to lie in the high- 
temperature phase and hence be noncritical, while all values — oo < w < —1/4 are predicted 
to lie in the Berker-Kadanoff phase and hence be critical with central charge c — —2 and 
leading thermal scaling dimension xt — 2 (corresponding to a marginally irrelevant operator) 
as given by (1.9)/(1.10) with 5 — 2. The transition between these two behaviors occurs at 

w = Wo{sq) = -1/4 , (1.11) 



*It is worth noticing [47] that there is a unified way of looking at the ferromagnetic and Berker-Kadanoff 
phases as continuations of one another. Let us parametrize (1.2±) by q = 4cos-^(7rM/2) and v = 2 cos(7r'«/2), 
with < u < 1 for the ferromagnetic critical curve and 1 < u < 2 for the Berker-Kadanoff curve. (We thus 
have u = 2/6 for the ferromagnetic phase, and u = 2 — 2/5 for the Berker-Kadanoff phase.) Then we can write 
the central charge as c = 1 — 3u^/ {2 — u), and the thermal scaling dimensions read XT,n = n{n + u)/{2 — u). 
Moreover, continuing these formulae to — 1 < u < gives the central charge and thermal scaling dimensions 
for the tricritical Potts model. (This variable u corresponds to the negative of the variable x employed 
in [33, Appendix A.l].) 
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where we expect a critical theory with central charge c = — 1 as given by (1.6) with 5 — 2. 
If (1.7)/(1.8) is correct, we should expect the thermal operator at w = —1/4 to be marginal: 
Xt = 2 and hence z/ = oo. But as we shall see (Sections 7.6 and 7.10), the prediction 
(1.7)7(1.8) is wrong, and a more likely scenario is xt = 0, so that u = 1/2. Finally, w = +oo 
is a ferromagnetic critical point with central charge c — —2 and leading thermal scaling 
dimension xt — 2 (corresponding to a marginally relevant operator), as given by (1.4)/(1.5) 
with 5 = 2. In a separate paper [24] we have shown that the w = +oo theory can be 
represented in terms of a pair of free scalar fermions, while the theory at finite w can be 
mapped onto a fermionic theory that contains a Gaussian term and a special four-fermion 
coupling; furthermore, this latter theory is perturbatively asymptotically free and is in fact 
equivalent, to all orders in perturbation theory in l/ty, to the A^- vector model at = — 1 
with (5 = —w. 

Triangular lattice: Baxter and collaborators [49-51] have determined the exact free en- 
ergy for the triangular-lattice Potts model on two special curves in the (5, ^;)-plane: 



These curves are plotted in Figure 2. The uppermost branch {y > 0) of curve (1.12) is known 
to correspond to the ferromagnetic phase-transition point, whose critical behavior is identical 
(thanks to universality) to that of the square-lattice ferromagnet [29-32]. The significance 
of the other two branches of (1.12) is not clear. The curve (1.13) is not critical in general, 
but it does contain the zero-temperature critical points at g = 2 [52] and q = A [53] (see [7] 
for further discussion and references). The existence of an antiferromagnetic critical curve 
for the triangular-lattice Potts model is at present not established; a fortiori its location, if 
it exists, is unknown. 

Consideration of the renormalization-group fiow for the triangular-lattice Potts model 
has led the present authors to hypothesize [25] that there exists an additional curve of 
RG fixed points — repulsive in the temperature direction — emanating from the point 
(g, v) = (0, 0) and extending into the antiferromagnetic region v < as g grows. In a 
subsequent paper [25] we shall present numerical estimates of the location and properties 
(e.g. critical vs. first-order) of this new phase-transition curve and discuss how it combines 
with the known antiferromagnetic critical points at g = 2,4 with v = —1, and with the 
first-order phase transition q = v ^ —0.79692 ± 0.00003 [15,54], to form a consistent 
phase diagram. One interesting issue is whether this curve (or part of it) might constitute 
a locus of critical points, in analogy with the case of the square lattice. For the time being 
we limit ourselves to the conjecture, in analogy with the square-lattice phase diagram, that 
the region lying in the range < 5 < 4 between the new phase-transition curve and the 
lower branch of curve (1.12) will constitute a Berker-Kadanoff phase whose critical behavior 
is determined by an attractive fixed point lying on the middle branch of curve (1.12). 

Note that Saleur [46, p. 669] expects universality both for the Berker-Kadanoff phase 
and for the critical theories forming its upper and lower boundaries. This would suggest, 
in particular, that the central charge in the Berker-Kadanoff phase of the triangular-lattice 
model might also be given by (1.9), and that the thermal scaling dimension might be given by 
(1.10). We have numerical evidence of the former result, which will be published elsewhere 
[25]. Moreover, if Saleur's conjecture is true, one would expect the lower branch of curve 



+ 3^2 - g = 
V = -I 



(1.12) 
(1.13) 
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(1.12) to be governed by the same critical continuum theory as the curve (1.3±) of the 
square-lattice Potts model. 

Let us denote by Wo(tri) the slope at g = of the new phase-transition curve. Then, 
if we consider the limit q,v ^ with w = v/q fixed, we find two different regimes: all 
values u'o(tri) < w < +oo are predicted to lie in the high-temperature phase and hence be 
noncritical, while all values —oo<w< u'o(tri) are predicted to lie in the Berker-Kadanoff 
phase and hence be critical with central charge c = —2 as given by (1.9) with S = 2. The 
transition between these behaviors occurs at w — u'o(tri), which may possibly constitute a 
critical theory of unknown type. 

Let us observe, finally, that the analytical results of [24] show that the conjectured uni- 
versality of the Berker-Kadanoff phase does hold at least in the limit q,v ^ with w — v/q 
fixed. Indeed, for all two-dimensional lattices, the Berker-Kadanoff phase ai —oo < w < wq 
is simply the c = —2 theory of a pair of free scalar fermions, perturbed by a four-fermion 
operator that is (in this phase) marginally irrelevant. 

Remark. It should be stressed that the Potts spin model has a probabilistic interpreta- 
tion (i.e., has nonnegative weights) only when g is a positive integer and v > —1. Likewise, 
the Fortuin-Kasteleyn random-cluster model [cf. (2.4) below], which reformulates the Potts 
model and extends it to noninteger g, has a probabilistic interpretation only when g > 
and > (or in the hmit considered here, w > 0). In all other cases, the model belongs 
to the "unphysical" regime (i.e., the weights can be negative or complex), and the ordinary 
statistical-mechanical properties need not hold. For instance, the free energy need not pos- 
sess the usual convexity properties; the leading eigenvalue of the transfer matrix need not be 
simple; and phase transitions can occur even in one-dimensional systems with short-range 
interactions.^ Nevertheless, there is a long history of studying statistical-mechanical models 
in "unphysical" regimes: examples include the hard-core lattice gas at its negative-fugacity 
critical point [56-66]; the closely related [56,57,62,63] problem of the Yang-Lee edge sin- 
gularity [4,67-73]; and the low-temperature A^- vector model at < 1, with application 
to dense polymers [74-85]. Indeed, the previously cited papers of Baxter [28,50,51] and 
Saleur [46, 47] deal in part with "unphysical" regimes in the Potts model. And though one 
must be especially careful in such studies, it generally turns out that the "unphysical" regime 
can be understood using the standard tools of statistical mechanics, appropriately modified. 
In particular, conformal field theory (CFT) seems to apply also in the "unphysical" regime, 
although there is (as yet) no rigorous understanding of why this should be the case: well- 
studied examples include the Yang-Lee edge singularity [72,73] and dense polymers [82-84]. 
Note that the "unphysical" nature of these models means that the corresponding CFT is 
non-unitary. 

Some aspects of the studies made in the present paper of the regime w < in the 
spanning-forest model — notably. Sections 7.5, 7.6, 7.9 and 7.10 — must therefore be un- 
derstood as relying imphcitly on such a conjectured extension of conformal field theory, 
analogously to the just-cited studies. On the other hand, the internal consistency of our 
results provides additional evidence for the validity of such an extension. 

^For a recent pedagogical discussion of the conditions under which one- dimensional systems with short- 
range interactions can or cannot have a phase transition, see [55]. 
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1.2 Outline of this paper 



The purpose of this paper is to shed hght on these phase diagrams in the neighborhood of 
{q, v) — (0, 0) by studying the ^ = Potts-model partition function in the complex w-plane 
for lattice strips of width m and length n, using a transfer-matrix method. For fixed width 
m and arbitrary length n, this partition function can be expressed via a transfer matrix of 
fixed size M x M (which unfortunately grows rapidly with the strip width m): 

Zmxniw) = tr[A{w)T{wr-'] (1.14a) 

M 

= Y,o^kHh{wr-' . (1.14b) 

k=l 

Here the transfer matrix T{w) and the boundary-condition matrix A{w) are polynomials 
in w, so that the eigenvalues {Xk} of T and the amplitudes {afc} are algebraic functions of 
w. We can of course use T{w) and A{w) to compute the zeros of the partition function 
for any finite strip m x n; but more importantly, we can compute the accumulation points 
of these zeros in the limit n oo, i.e. for the case of an semi-infinite strip [5-7,86-89]. 
According to the Beraha-Kahane- Weiss theorem [90-92], the accumulation points of zeros 
when n — > oo can either be isolated limiting points (when the amplitude associated to the 
dominant eigenvalue vanishes, or when all eigenvalues vanish simultaneously) or belong to 
a limiting curve B (when two dominant eigenvalues cross in modulus). As the strip width 
m tends to infinity, the curve B = B„i is expected to tend to a thermo dynamic-limit curve 
Boo, which we interpret as a phase boundary in the complex w-plane. In particular. Boo is 
expected to cross the real axis precisely at the physical phase-transition point wq. 

Our approach is, therefore, to compute the curves Bm up to as large a value of m as our 
computer is able to handle, and then extrapolate these curves to m = oo. One output of our 
study is a numerical estimate of Wq. For the square lattice, we find 

wo(sq) = -0.2501 ± 0.0002 , (1.15) 

in excellent agreement with the prediction Wo(sq) = —1/4. We are therefore justified in 
assuming, in our subsequent analysis, that Wo(sq) = —1/4 exactly. For the triangular lattice, 
we find 

wo(tri) = -0.1753 ± 0.0002 . (1.16) 

We do not yet know whether the exact value of u'o(tri) is given by any simple closed-form 
expression. Nor do we know whether there exists a simple exact formula for the location of 
the critical curve Boo in the complex ui-plane, for either the square or triangular lattice. 

In order to shed further light on the phase diagram of the square- and triangular-lattice 
systems, we have studied the free energy (and its derivatives with respect to w), the central 
charge c and the thermal scaling dimension xt as a function of w on the real ty-axis. We find 
that at w — Wq there is a first-order phase transition: the free energy has a discontinuous 
first derivative with respect to w at w = Wq. From our numerical work it is not clear whether 
the discontinuity in this first derivative is finite or infinite: the analysis seems to favor a finite 
limit, but a weak divergence such as f'{w) ~ log(-u; — Wq) or f'{w) ~ loglog('u; — Wq) is also 
possible. The two phases are characterized as follows: 

• w < Wq: In this case the system is critical and it can be described by a conformal field 
theory with central charge c = —2. The leading thermal scaling exponent is xt = 2, 
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so that the energy is a marginal operator. Both results agree with Saleur's prediction 
(1.9)/(1.10) for the Berker-Kadanoff phase. 

• w > wq: In this case the system is non-critical, i.e., the correlation length is finite. 

The behavior at w = Wq is rather special. On the one hand, it is a coexistence point of 
two different phases; on the other hand, it is itself a critical point, which belongs to a different 
universality class from that of the Berker-Kadanoff phase, namely the one corresponding to 
the — > limit along the antiferromagnetic critical curve (1.3+) . Thus, at least for the 
square lattice, we expect that this point will be described by a conformal field theory of 
central charge c = — 1 in accordance with Saleur's prediction (1.6). 

Our numerical conclusions concerning the behavior at w = Wq are drawn principally 
from the results on the square lattice (as in this case we know the exact location of Wq), 
but we expect them to hold also for the triangular lattice. On the square lattice we have 
clear evidence that Wo(sq) = —1/4 is indeed critical, and we can give rough estimates of the 
central charge and thermal scaling dimensions: 

1) Extrapolation of small-w expansions using differential approximants (Appendix A) 
shows that the second derivative of the free energy diverges for w I wq as {w — wq) 
This exponent is close to the theoretical prediction a = 1 for a first-order critical point 
(see below), i.e. f"{w) ~ (w — Wq)~^, and is consistent with it if one does not take 
the alleged error bar too seriously (the error estimates in series extrapolation have no 
strong theoretical basis). We have checked (Section 7.4) that our data from finite- width 
strips are consistent with this latter behavior, possibly modified by a multiplicative log- 
arithmic correction. 

2) The central charge atw — Wq is different from both the Berker-Kadanoff value (c = —2) 
and the noncritical value (c = 0). Indeed, we get estimates around c ~ —1.3, which 
seem to be tending roughly towards c = — 1 as the strip width grows (Section 7.5). 

3) The correlation length at w — Wq is infinite for all even strip widths, reflecting the fact 
that w = wq'is an exact endpoint of the limiting curve B for even widths. For odd strip 
widths, the correlation length at w = wq is flnite and seems to be tending to infinity 
as the strip width tends to infinity (Section 7.6). 

4) Because the correlation length at w = wq is infinite for all even strip widths, we 
conclude from standard CFT arguments [93] that the leading thermal scaling dimension 
xt,i equals 0, and hence that v = 1/2 (Section 7.6). The scaling law dv — 2 — a then 
yields the specific- heat exponent o; = 1, in agreement with the results from small-u; 
expansions. 

5) The second thermal scahng dimension xt,2 can be obtained from the second gap for 
even widths (Section 7.6). We obtain a value xt,2 ^ 0.47 that is clearly different from 
the Berker-Kadanoff value xt,i = 2. This result is close to the one obtained from the 
first gap for odd widths: xt,2 ^ 0.37. 

In conclusion, the phase transition at w — Wq is rather unusual. The point w — Wq is 
what Fisher and Berker [94, pp. 2510-2511] have called a first-order critical point: namely, 
it is both a first-order transition point (the first derivative of the free energy is discontinuous 
at w = Wq) and critical (the correlation length is infinite for w < Wq and diverges as w | Wq). 
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The critical exponents take the values a — 1 and v — 1/d, in agreement both with finite- 
size-scaling theory for first-order transitions (z/ = 1/d and a/v — d) [94,95] and with the 
hyperscaling law dv = 2 — a ioi critical points. ^'^ 

This sort of phase transition seems to be rare; indeed, we are aware of only two other 
equilibrium-model examples: 

1) One- dimensional q-state clock model with a 6 -term [96-99]. In this model, the first 
derivative of the free energy has a finite discontinuity at the transition points, and the 
correlation length diverges there with exponent v = 1/d = 1. However, the specific heat 
does not diverge at the transition points: it is a discontinuous but bounded function of 
the temperature (see [99] for a computation of the specific heat in the limit q — > oo). 
This behavior does not contradict the Fisher-Berker scahng theory for first-order critical 
points [94, p. 2511]: "first-order critical points are no more than ordinary critical points in 
which either one or both of the relevant eigenvalue exponents attains its thermodynamically 
allowed limiting value A = . Unfortunately the magnetic exponent was not considered in 
ref. [99]. 

2) Six-vertex model (which contains the KDP model of a ferroelectric [100, 101] as a 
special casc).^^ We follow the notations of Baxter's book [27, Sections 8.10 and 8.11]. The 
transition between any of the two fcrroclcctrically ordered phases (regimes 1 and 11 in Baxter's 
notation) and the critical phase (regime III) is characterized by the following properties: 

a) The first derivative of the free energy with respect to the temperature, /gy, has a finite 
jump discontinuity on the transition line. 

b) The second derivative of the free energy with respect to the temperature, /gy, is iden- 
tically zero in the ferroelectrically ordered phases and diverges like (T — Tc)"^^^ as we 
approach the transition line from the critical phase. 

c) The correlation length is infinite in the critical phase (i.e., the correlations decay alge- 
braically to zero) and identically zero in the ferroelectrically ordered phases. 

d) The electric polarization is identically zero in the critical phase and identically 1 in the 
ferroelectrically ordered phases. 

Because of (a) and (b), it is unclear whether we should say Q; = lorQ; = l/2. And because 
of (c), the critical exponent v cannot be sensibly defined. Taken together, the properties 
(a)-(d) are rather unusual. 

It would be interesting to know of other examples of first-order critical points. 

Please note that any two of the four equations a = 1, v = 1/d, a/v = d and dv = 2 — a imply the other 

two. 

^^We are grateful to an anonymous referee for bringing this model to our attention. 

'^^The one-dimensional Ising fcrromagnet with l/r^ long-range interaction, cited by Fisher and Bcrkcr 
[94, p. 2510] as an example of this phenomenon, does not, strictly speaking, qualify, as the (exponential) 
correlation length is -|-oo for all /3 > (by Griffiths' second inequality, {o^Oy) > tanh Jj;j, w /3/|x — yp), 
so that a critical exponent ly cannot be defined. Miillcr [102] studied a onc-dimcnsional model of U{N)- 
or 5?7(-/V)- valued spins with a complex nearest-neighbor interaction, generalizing the work on the g-state 
clock model with a 0-term. He showed that for U (N) , the model undergoes a sequence of first-order phase 
transitions; however, it is not clear to us whether the correlation length diverges at those transition points. 
Finally, Oliveira and coworkers [103, 104] have found that some non- equilibrium models exhibit a first-order 
critical point with v = 1/d. 
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The plan of this paper is as follows: In Section 2 we review the needed background 
concerning the Potts-modcl partition function and the combinatorial polynomials that can 
be obtained from it. In Section 3 we review the small-w and large-w expansions for the 
g — > Potts model. In Section 4 we summarize how we computed the transfer matrices. 
In Sections 5 and 6 we report our results for square-lattice and triangular-lattice strips, 
respectively. In Section 7 we analyze the data to extract estimates of the critical point wq, 
the free energy f{w), the central charge c{w), and the thermal scaling dimensions XT,i{w) 
for each of the two lattices. Finally, in Section 8 we discuss some open questions. In the 
Appendix we discuss how we performed the analysis of the small-w series expansions obtained 
in Section 3. 



2 Basic set-up 

We begin by reviewing the Fortuin-Kasteleyn representation of the g-state Potts model 
(Section 2.1). Then we discuss the various polynomials that can be obtained from the 
g-state Potts model by taking the limit q ^ (Section 2.2). After this, we define the 
specific quantities that we will be studying in our work on spanning forests on the square 
and triangular lattices, and review the basic principles of finite-size scaling in conformal 
field theory (Section 2.3). Finally, we review briefly the Beraha-Kahane-Weiss theorem 
(Section 2.4). 



2.1 Fortuin-Kasteleyn representation 

Let G = {v. E) be a finite undirected graph with vertex set V and edge set E; let 
{Je}eeE be a set of couplings; and let g be a positive integer. At each site i E V we place 
a spin (Tj e {1,2, ... ,q}, and we write a = {crjjjgv to denote the spin configuration. The 
Hamiltonian of the g-state Potts model on G is 

H{(t) = - E JeS{ai,aj), (2.1) 

e=ij£E 

where 5 denotes the Kronecker delta. The partition function Z = q-(3H{(t) written 
in the form 

^g(9,v) = E n [l + ^'e5(a„a,)] , (2.2) 

<T e=ijeE 

where 

= e^-^^ - 1 (2.3) 

and V = {ve}e(zE- A coupling Je (or Ve) is called ferromagnetic if Jg > {ve > 0), antiferro- 
magnetic if — oo < Jg ^ (—1 < Vg < 0), and unphysical if ^ [—1, oo). 

At this point q is still a positive integer. However, we now assert that ^(^(g, v) is in fact 
the evaluation at g e Z+ of a polynomial in q and {ve} (with coefficients that are indeed 
either or 1). To see this, we proceed as follows: In (2.2), expand out the product over 
e E E, and let A C. E be the set of edges for which the term VeS{ai,aj) is taken. Now 
perform the sum over configurations cr: in each connected component of the subgraph [V, A) 
the spin value cxj must be constant, and there are no other constraints. Therefore 

ZG{q,v) = E^'^^^n^- (2.4) 

ACE eeA 
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where k{A) is the number of connected components (including isolated vertices) in the sub- 
graph {V,A). The subgraph expansion (2.4) was discovered by Birkhoff [105] and Whit- 
ney [106] for the special case Ve = —1 (see also Tutte [107,108]); in its general form it is 
due to Fortuin and Kasteleyn [109,110] (see also [111])- Henceforth we take (2.4) as the 
definition of Zciq, v) for arbitrary complex q and v. When Ve takes the same value v for all 
edges e, we write ZG{q,v) for the corresponding two-variable polynomial. 

Let us observe, for future reference, that (2.4) can alternatively be rewritten as 

Zaiq.v) = gl^l E ^^^^^ n - ' (2-5) 

ACE e€A ^ 

where we use the notation IS"! to denote the number of elements of a finite set S, and 

c{A) = \A\+k{A)-\V\ (2.6) 



is the cyclomatic number of the subgraph (V^, ^4), i.e. the number of linearly independent 
circuits in {V, A). 

Remark. In the mathematical literature, these formulae are usually written in terms of 
the Tutte polynomial Tq defined by [112] 

TG{x,y) = J](a;-l)'=(^)-^(«)(|/-l)^(^) (2.7) 

ACE 

where k{G) = k{E) is the number of connected components in G, and c{A) is the cyclomatic 
number defined in (2.6). Comparison with (2.4) shows that 

TG{x,y) = (x-l)-^(«)(|/-l)-I^IZG((x-l)(y-l),y-l). (2.8) 

In other words, the Tutte polynomial Tg{x, y) and the Potts-model partition function Zoiq, v) 
are essentially equivalent under the change of variables 

X = 1 + q/v (2.9a) 
y ^ 1 + v (2.9b) 



q = {x-l){y-l) (2.9c) 
V ^ y-1 (2.9d) 

The advantage of the Tutte notation is that it allows a shghtly smoother treatment of the 
q ^ limit (see Remark 1 at the end of the next subsection). The disadvantage is that 
the use of the variables x and y conceals the fact that the particular combinations q and v 
play very different roles: g is a global variable, while v can be given separate values Ve on 
each edge. This latter freedom is extremely important in many contexts (e.g. in using the 
series and parallel reduction formulae). We therefore strongly advocate the multivariable 
approach in which ZG{q, v) is considered the fundamental quantity, even if one is ultimately 
interested in a particular two-variable or one-variable specialization. 
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2.2 q ^ limits 

Let us now consider the different ways in which a meaningful g — > hmit can be taken 
in the g-state Potts model. 

The simplest limit is to take g — with fixed couplings v. From (2.4) we see that this 
selects out the subgraphs ACE having the smallest possible number of connected compo- 
nents; the minimum achievable value is of course k{G) itself (= 1 in case G is connected). 
We therefore have 

limg-'=(«)ZG(g,v) = Cg(v), (2.10) 

where 

Cg{^) = E (2.11) 

ACE eeA 
k{A) = k{G) 

is the generating polynomial of "maximally connected spanning subgraphs" (= connected 
spanning subgraphs in case G is connect ed).^^ 

A difi^erent limit can be obtained by taking g — > with fixed values of w = v/g. From 
(2.5) we see that this selects out the subgraphs ACE having the smallest possible cyclomatic 
number; the minimum achievable value is of course 0. We therefore have [113, 114] 

limg-I^IZG(g,gw) = Fg{w), (2.12) 

where 

FgM = E (2.13) 

ACE eSA 
c(A) = 

is the generating polynomial of spanning forests , i.e. spanning subgraphs not containing any 
circuits. 

Finally, suppose that in Cg(v) we replace each edge weight Ve by Xve and then take A ^ 0. 
This obviously selects out, from among the maximally connected spanning subgraphs, those 
having the fewest edges: these are precisely the maximal spanning forests (= spanning trees 
in case G is connected), and they all have exactly \V\ — k{G) edges. Hence 

limA'=(^)-I^ICG(Av) = Tg(v), (2.14) 

where 

Toiy) = E (2.15) 

ACE eeA 
k{A) = k(G) 
c{A) = 

is the generating polynomial of maximal spanning forests (= spanning trees in case G is 
connected). Alternatively, suppose that in Fg(w) we replace each edge weight Wg by Xw^ 

^"^A subgraph is called spanning if its vertex set is the entire set V of vertices of G (as opposed to some 
proper subset of V). All of the subgraphs arising in this paper are spanning subgraphs, since we consider 
subsets A of edges only, and retain all the vertices. 

^^We trust that there will be no confusion between the generating polynomial T'g(v) and the Tutte 

polynomial Tg{x, y). We have used here the letter T because in the most important applications the graph 
G is connected, so that Tg(v) is the generating polynomial of spanning trees. 
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and then take A — > oo. This obviously selects out, from among the spanning forests, those 
having the greatest number of edges: these are once again the maximal spanning forests. 
Hence 

lim A'^^^^-'^IFgIAw) = Tg(w) . (2.16) 

A— ►oo 

In summary, we have the following scheme for the g — > limits of the Potts model: 

Cg(v) 



q ^ 0, -V fixed v infinitesimal 



Zo(?,v) To(v or w) (2.17) 

g 0, w = v/g fixed w infinite 



Finally, maximal spanning forests (= spanning trees in case G is connected) can also be 
obtained directly from Zaiq, v) by a one-step process in which the limit g — > is taken at 
fixed X = v/g"^, where < a < 1 [110,113-115]. Indeed, simple manipulation of (2.4) and 
(2.6) yields 

Zoiq, g"x) = g"l^l J2 g"^(^)+(i-")'=(^) J] Xe . (2.18) 

AQE eeA 

The quantity ac{A) + (1 — a)k{A) is minimized on (and only on) maximal spanning forests, 
where it takes the value (1 — a)k{G). Hence 

limg-"l^l-(i-")'=(^)ZG(g,g"x) = Tg{^) . (2.19) 



Remarks. 1. Let us rewrite the formulae of this subsection in terms of the Tutte 
polynomial Tcix.y) [cf. (2.7)-(2.9)]. The limit q,v ^ with w — v/q fixed corresponds to 
X — l + l/w,y — 1; simple algebra using (2.8) and (2.12) gives 

Tg(1 + l/w, 1) = w''^^^-\^\Fg{w) . (2.20) 

In particular, several of the evaluations of Tg(1 + l/'u;,l) have interesting combinatorial 
interpretations: 

• Tg(1, 1) (i.e., w — oo) counts the number of maximal spanning forests in G {— spanning 
trees if G is connected). 

• Tg(2, 1) (i.e., w — 1) counts the number of spanning forests in G. 

• Tg(1 + l/k,l) (i.e., w = k) for k = 1,2,3, .. . is, up to a prefactor, the number of 
possible "score vectors" in a tournament of constant-sum games with scores lying in 
the set {0, 1, ... , k} [116, Propositions 6.3.19 and 6.3.25]. 
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• Tg{0, 1) (i.e., w — —1): If G is a directed graph having a fixed ordering on its edges, 
\Tg{0, 1)1 counts the number of totally cyclic reorientations r of G such that in each 
cycle of r the lowest edge is not reoriented. For a planar graph with no isthmuses, 
\Tg{0, 1)1 counts the number of totally cyclic orientations in which there is no clockwise 
cycle. See [116, Examples 6.3.30 and 6.3.31]. 

We emphasize, however, that from a physical point of view there is nothing special about 
the particular values w = cxo, 1, 2, 3, . . . , —1. Rather, it is important to study Fg{w) as a 
function of the real or complex variable w. The "special" values of w are those lying on the 
phase boundary B; they are determined only a posteriori. 

2. The polynomial Fg{w) also equals the Ehrhart polynomial of a particular unimodular 
zonotope determined by the graph G: see [117, Section XI.A] for details. 

3. Suppose that G — (V, E) is a connected planar graph; then we can define a dual 
graph G* — {V*, E*) by the usual geometric construction.^^ Moreover, there is a one-to-one 
correspondence between the edges of G and their corresponding dual edges in G*; we can 
therefore identify E* with E, and assign the same weights v = {ve}eeE to the edges of G 
and G*. We then have the fundamental duality relation [2] 

ZG*{q,^) = g^-'^' ( ) ZGiq,qM (2.21) 

\e€E J 

[where g/v of course denotes the vector {^/vejeGe]- In particular, we have 



Cg*{v) = (n^e)i^G(l/v) (2.22) 

\eeE J 

Fg*{w) = (1[wAGg{IM (2.23) 

\eeE J 

Tg*{v) = (n^^e]TG(l/v) (2.24) 



4. Suppose that G is connected; and let us consider G as a communications network with 
unreliable communication channels, in which edge e is operational with probability pe and 
failed with probability 1 — Pe, independently for each edge. Let Rg{p) be the probability 
that every node is capable of communicating with every other node (this is the so-called 
all-terminal reliability). Clearly we have 

Mp) = E Up- Yl - Pe) ^ (2-25) 

ACE eeA eeE\A 

(V, A) connected 



^^More precisely, consider a particular plane representation of G, and define G* by placing one vertex in 
each face of G and then drawing an edge of G* through each edge of G. The dual graph G* is not necessarily 
unique; nonisoniorphic plane representations of G (which can arise if G is not 3-connectcd) can give rise to 
nonisomorphic duals (see e.g. [118, p. 114] for an example). In any case, each of the dual graphs G* satisfies 
the relations (2.21)-(2.24) below. 
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where the sum runs over all connected spanning subgraphs of G. The polynomial Rg{v) 
is called the (multivariate) reliability polynomial [119] for the graph G. Modulo trivial 
prefactors, it is equivalent to Cg(v) under the change of variables 

^ Pe 

Pe = (2.26b) 

The reliabihty polynomial is therefore one of the objects obtainable as a ^ — > limit of the 
Potts model. 

5. Brown and Colbourn [120], followed by Wagner [121] and Sokal [122], have studied 
the possibility that the complex roots of the reliability polynomial might satisfy a theorem 
of Lee- Yang type. To state what is at issue, let us say that a graph G has 

• the univariate Brown-Colhourn property if Gg{v) ^ whenever |1 + < 1; 

• the multivariate Brown-Colhourn property if Cg(v) ^ whenever |1 + i>e| < 1 for all 
edges e; 

• the univariate dual Brown-Colhourn property if Fciw) ^ whenever Kew < —1/2; 

• the multivariate dual Brown-Colhourn property if Fg(w) ^ whenever Rewe < —1/2 
for all edges e. 

Here the word "dual" refers to the fact that a planar graph G has the (univariate or mul- 
tivariate) Brown-Colbourn property if and only if its dual graph G* (which is also planar) 
has the (univariate or multivariate) dual Brown-Colbourn property: this is an immediate 
consequence of the identities (2.22)/(2.23). 

Brown and Colbourn [120], having studied the univariate reliability polynomial Rg{p) 
in a number of examples, conjectured that every loopless graph has the univariate Brown- 
Colbourn property.^^ (Of course, they didn't call it that!) Subsequently, Sokal [122] made 
the stronger conjecture that every loopless graph has the multivariate Brown-Colbourn 
property. Moreover, Sokal [122] proved this latter conjecture for the special case of series- 
parallel graphs (which are a subset of planar graphs). And since the class of series-parallel 
graphs is self-dual, it follows immediately that every bridgeless series-parallel graph has 
the multivariate dual Brown-Colbourn property.^® Finally, Sokal's conjecture would imply 

^^A loop is an edge connecting a vertex i to itself. Loops must be excluded because a loop e multiples Cq 
by a factor 1 + Ve and therefore places a root at Ve = — 1, violating the Brown-Colbourn property. 

"'^''A graph is called series-parallel if it can be obtained from a forest by a sequence of series and parallel 
extensions (i.e. replacing an edge by two edges in series or two edges in parallel). The proof that every loopless 
series-parallel graph has the multivariate Brown-Colbourn property is an almost trivial two-line induction; 
it can be found in [122, Remark 3 in Section 4.1]. Earlier, Wagner [121] had proven, using an ingenious 
and complicated construction, that every loopless series-parallel graph has the univariate Brown-Colbourn 
property. 

^^A bridge is an edge whose removal increases (by 1) the number of connected components of the graph. 
Bridges must be excluded because a bridge e multiplies Fq by a factor \ + We and therefore places a root at 
We — —1, violating the dual Brown-Colbourn property. Please note that a planar graph G is loopless (resp. 
bridgeless) if and only if its dual graph G* is bridgeless (resp. loopless). 



17 



that every bridgeless planar graph (series-parallel or not) has the multivariate dual Brown- 
Colbourn property. These conjectures, if true, would constitute a powerful result of Lee-Yang 
type, constraining the complex zeros of the corresponding partition functions. 

Recently, however, Royle and Sokal [123] have discovered — to their amazement — 
that there exist planar graphs for which all these properties fail! Indeed, the multivariate 
Brown-Colbourn and dual Brown-Colbourn properties fail already for the simplest non- 
series-parallel graph, namely the complete graph on four vertices {K^. The univariate 
properties fail for certain graphs that can be obtained from Kt^ by series and/or parallel 
extensions. In fact, Roylc and Sokal [123] show that a graph G has the multivariate Brown- 
Colbourn property if and only if it is series-parallel. 

In addition, Chang and Shrock [23] , building on one of the Royle-Sokal examples, have 
devised families of strip graphs in which the limiting curve B of zeros of Cg{v) penetrates 
into the "forbidden region" |1 + f | < 1. 

Nonetheless, the strip graphs studied in this paper do seem to possess at least the uni- 
variate dual Brown-Colbourn property: all the roots we find lie in the region Kew > —1/2. 



2.3 Quantities to be studied 

The graphs G to be considered in this paper are m x n strips of the square or triangular 
lattice, with periodic boundary conditions in the first (transverse) direction and free bound- 
ary conditions in the second (longitudinal) direction. We therefore denote these strips as 
mp X n-p, and call this cylindrical boundary conditions}^ We shall also use the letter L as 
an alternative name for the strip width m. All these graphs are planar. 

In this paper we will be focussing on Fg(w), the generating polynomial of spanning 
forests. Since we will be assigning the same weight w to all nearest-neighbor edges, we have 
a univariate polynomial Fg{w). We will refer to Fg{w) as the "partition function" (since 
that is basically what it is). Our principal goal is to study the behavior of Fg{w) in the 
thermodynamic limit m,n ^ oo. 

Let us now define the free energy (or "entropy") per site for a finite lattice^° 

fm,n{w) = — logF^px„p(w) (2.27) 

mn 

and its limiting values for a semi-infinite strip 

/„(w) = lim — logF^pxnp(w) (2.28) 
n-^oo mn 



and for the infinite lattice 



/(«;)= lim J-logF„pxnp(H ■ (2-29) 

m,n— >oo mn 



^^This accords with the terminology of Shrock and collaborators [124] for the various boundary conditions: 
free (mp x tif), cylindrical (mp x np), cyclic (mp x np), toroidal (mp x np), Mobius (mp x nxp) and Klein 
bottle (mp X riTp)- Here F denotes "free", P denotes "periodic", and TP denotes "twisted periodic" (i.e. 
the longitudinal ends are identified with a reversal of orientation). 

^°Note that our "free energy" is the negative of the usual free energy. 
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Here we are assuming, of course, that the indicated hmits exist and that in (2.29) the hmit 
is independent of the way that m and n tend to infinity.^^ In particular, in (2.29) we can 
allow n to tend to infinity first and then take m — > oo, so that 

f{w) = hm Uw) . (2.30) 

m— >oo 

Let us also note that for w negative or complex these formulae may contain some ambiguities 
about the branch of the logarithm. We nevertheless expect fm{w) and f{w) to be well-defined 
analytic functions in the complex w-plane minus certain branch cuts; but for simplicity 
we shall mostly focus on the real part of the free energy, which does not suffer from any 
ambiguities. 

For each fixed width m, the partition function Fj^pxripiw) for strips of arbitrary length 
n can be expressed in terms of a transfer matrix: 

i^mpxnp(«;) = u{wfT{w)^-Wid{w) (2.31a) 

M 

= J2Mw)Xk{wr-^ (2.31b) 

k=l 

where M = dimT is the dimension of the transfer matrix. This is explained in detail in 
Ref. [5] and is summarized very briefiy in Section 4 below. The elements of the transfer 
matrix T{w) and of the left and right vectors u{w) and Vid('u;) are polynomials in the com- 
plex variable w. Therefore, the eigenvalues Xkiw) and the amplitudes ak{w) are algebraic 
functions of w. We have numerically checked for m < 9 that none of the amplitudes ak{w) 
vanishes identically. This is important in order to compute the limiting curves B, as one 
would get the wrong curve B if there were an identically vanishing amplitude for an eigenvalue 
that happened to be dominant in some region of the w-plane. 

Let us denote by A*(i/;) the eigenvalue of T{w) having largest modulus, whenever it is 
unique. (Typically there is a unique dominant eigenvalue at all points in the complex w- 
plane with the exception of a finite union of real algebraic curves B.) It follows from (2.31b) 
that the strip free energy (2.28) exists at all such points w — except possibly at isolated 
points where the amplitude ccfe corresponding to the dominant eigenvalue vanishes — and 
equals 

hiw) = ^\ogX,{w). (2.32) 

[In particular, Re/i(w;) = (1/L) log |A^(-u;)|.] We then expect fhiw) to converge as L ^ cxo 
to the bulk free energy /(w); but the rate at which it converges depends on whether the 

^-'^For iTj > 0, it is not hard to prove rigorously that the hmit (2.29) exists, at least if we insist that m and 
n tend to infinity in such a way that the ratio m/n stays bounded away from zero and infinity. The proof is 
based on the submultiplicativity of Fa{w) for disjoint regions of the lattice, together with a standard result 
on subadditive functions [125, Proposition A. 4]. 

This proof, by itself, says nothing about w negative or complex. But the convergence can in some cases be 
extended to part of the complex w-plane, by using a normal- families argument [126,127]. Indeed, suppose 
that D C C is a connected open set having a nonempty intersection with the positive real axis, on which the 
partition function Fmxniw) is nonvanishing for all (or all sufficiently large) m,n. Then the trivial bound 
\Fmxn{w)\ < (1 + |w|)^™" guarantees that Re fm,n{w) is uniformly bounded above on compact subsets of 
D, from which it follows [127, Example 2.3.9] that the analytic functions {fm,n} form a normal family on 
D. Then a Vitali-like argument [92, Lemma 3.5] shows that the convergence for w > extends to all of D. 
Simon [128, p. 343] calls this reasoning "log exp Vitali" . 
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model at w is critical or not. If w is a noncritical point, we expect an exponentially rapid 
convergence: 

hiw) - f{w) + 0(e-^/«) , (2.33) 

where ^ < cx) is the correlation length of the system. If w is a critical point, then we expect 
that its long-distance behavior can be described by a conformal field theory (CFT) [30-32] 
with some central charge c{w); the general principles of CFT then predict [129, 130] that^^ 

AW = /H + f^ + .... (2.34) 

where G is a geometrical factor depending on the lattice structure, 

^ r 1 square lattice ,^ 

\ \/3/2 triangular lattice 

and the dots stand for higher-order corrections. These higher-order corrections always in- 
clude a 1/L'^ term (with a nonuniversal amplitude) coming from the operator TT, where T 
is the stress-energy tensor; sometimes irrelevant operators may give additional corrections 
in-between and 1/L^. 

Let the eigenvalues of the transfer matrix (at some particular value of w and some 
particular strip width L) be ordered in modulus as |A*| > |Ai| > IA2I > ■ ■ ■, and let us 
define correlation lengths ^ > • • • by 



= log 



A. 



(2.36) 



A, 

Then CFT predicts [93] that, at a critical point, the correlation lengths behave for large L 
as 

= 2nG^ + . . . , (2.37) 

where Xi is the scaling dimension of an appropriate scaling operator, G is the geometrical 
factor (2.35) [131], and the dots stand again for higher-order corrections. 

Remcirks. 1. In order to compare our results more directly with those of mathematicians 
(e.g., [132]), it is convenient to introduce also the quantities 

n^{w) = e^''^'"^"') = lA^w;)!^/"* (2.38a) 
n{w) = lim nm{w) = e^''^^'"^ (2.38b) 

m— »oo 

2. For the square lattice at w — 1, Calkin et al. [132] have proven the bounds 

1.29335 !^ log 3.64497 < /(sq,w = l) < log 3.74101 ^ 1.31936. (2.39) 

Their proof uses an m x m lattice with free boundary conditions in both directions, but the 
same bound for cylindrical boundary conditions is an easy corollary. Weaker bounds of the 



^^Note that the 1/Lp' correction in (2.34) has the opposite sign from [129, equation 1]. This change 
compensates the global change of sign introduced in our definition of the free energy (2.27), so that the 
central charge c(w) in (2.34) has the conventional sign. 
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same type were proven earlier by Merino and Welsh [133]. The upper bound in (2.39) comes 
from the inequahty [132, Theorem 5.3] 

/(sq, w) < — log 1(1 + w)Af *^'^''")(w)| for all m > and w>0, (2.40) 

where x^^^'^'"^^ (^ j is the dominant eigenvalue of the transfer matrix for a square-lattice strip 
of width m with free boundary conditions. (The proof given in [132] for w = 1 generalizes 
immediately to any w > 0.) The upper bound in (2.39) was obtained by evaluating the right- 
hand side of (2.40) for m = 8. We have slightly improved this upper bound by computing 
the dominant eigenvalue at = 1 for m = 10:^^ 

1.29335 ^ log 3.64497 < /(sq,'w; = l) < log 3.73264 ^ 1.31711. (2.41) 

3. For the triangular lattice one can compute a similar upper bound for /(tri,w) by 
mimicking the derivation of (2.40). One obtains the following rigorous upper bound: 

/(tri, w) < — log 1(1 + w)'^X^'^'"'\w)\ for all m > and w>0. (2.42) 
m 

We have obtained a numerical bound by computing the leading eigenvalue of the transfer 
matrix at w = 1 for m = 9.^^ A trivial lower bound can be obtained in terms of the 
entropy per site for spanning trees on the triangular lattice, /o(tri), given in (3.16) below 
[114, 133-135]. Putting these bounds together, we have 

1.61533 w /o(tri) < /(tri,w-l) < log 5.77546 1.75362. (2.43) 

2.4 Beraha—Kahane— Weiss theorem 

A central role in our work is played by a theorem on analytic functions due to Beraha, 
Kahane and Weiss [86,87,90,91] and generalized slightly by one of us [92]. The situa- 
tion is as follows: Let D be a domain (connected open set) in the complex plane, and let 
«!, . . . , aM, Ai, . . . , Am {M > 2) be analytic functions on D, none of which is identically zero. 
For each integer n > 0, define 

M 

U{z) = J2Mz)>^k{zr- (2-44) 

k=l 

We are interested in the zero sets 

2{fn) = {zeD: Uz) = 0} (2.45) 

^^We have obtained the symbohc form of the transfer matrices for square-lattice strips with free boundary 
conditions up to widths L < 9f. The dimensions of such matrices were obtained in closed form in [14, 
Theorem 5], and they are listed on column SqPree(L) in Table 1. For L = IOf we have obtained numerically 
the transfer matrix at w = 1. Details of these transfer matrices will be given elsewhere. 

^^We have obtained the symbolic form of the transfer matrices for triangular-lattice strips with free bound- 
ary conditions up to widths i < Sp. The transfer-matrix dimension is the Catalan number Cz,; they are 
listed on column TriFree(L) in Table 1. For i = 9f we have obtained numerically the transfer matrix at 
w = 1. Details of these transfer matrices will be given elsewhere. 
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and in particular in their limit sets as n — > oo: 

liminf = {z ^ D: every neighborhood U 3 z has a nonempty intersection 

with all but finitely many of the sets Z{fn)} (2.46) 

limsupZ(/„) = {z & D: every neighborhood U 3 z has a nonempty intersection 

with infinitely many of the sets Z(/„)} (2.47) 

Let us call an index k dominant at z if |Afc(2;)| > |A;(;s)| for all / (1 < / < M); and let us 
write 

-Dfe = {z E D: k is dominant at z} . (2.48) 
Then the limiting zero sets can be completely characterized as follows: 

Theorem 2.1 [86, 87, 90-92] Let D be a domain in C, and let ai, . . . , q^m, ^i, . . . , Am 

(M >2) be analytic functions on D, none of which is identically zero. Let us further assume 
a "no-degenerate-dominance" condition: there do not exist indices k ^ k' such that Xk = utXk' 
for some constant cu with |ci;| = 1 and such that Dk (— Dk' ) has nonempty interior. For each 
integer n > 0, define fn by 

M 

fn{z) = 5^a,(^)Afc(^)V 

k=l 

Then liminf Z{fn) = limsup 2^(/n), and a point z lies in this set if and only if either 

(a) There is a unique dominant index k at z, and a^iz) = 0; or 

(b) There are two or more dominant indices at z. 

Note that case (a) consists of isolated points in while case (b) consists of curves (plus 
possibly isolated points where all the A^ vanish simultaneously) . Henceforth we shall denote 
by B the locus of points satisfying condition (b) . 

Remark. For the strip lattices considered in this paper, we have looked for isolated limiting 
points using the determinant criterion [5]. For L < 5, we have found that there are no such 
limiting points: the dominant amplitude does not vanish anywhere in the complex w-plane. 
For L > 6, wc were unable to compute the needed determinants due to memory limitations. 
But our computations of zeros for finite-length strips (Figures 5 and 9) give no indication 
of any isolated limiting points. We therefore conjecture that there are no isolated limiting 
points for any of the strip lattices studied in this paper. 

3 Small-i(; and large-^i; expansions 

In this section we discuss the small-ti; and large-w expansions for the spanning-forest 
model. 
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3.1 Small-w; (high-temperature) expansion 

Let us begin by considering spanning forests on an L x L square lattice with periodic 
boundary conditions in both directions. For small k, it is easy to count the number of ways of 
making a A;-edgc spanning forest; the result is a polynomial in the volume V = L"^ provided 
that one restricts attention to L > A; + 1 so as to avoid clusters that "wind around the 
lattice" . For example, doing this for k < 8 one finds 



= 1 + 2Vw + 
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(3.1) 



valid for L > 9. Taking the logarithm, one finds that each coefficient is proportional to V: 
all higher powers of V have cancelled. Dividing by V, one obtains the small-w expansion of 
the bulk free energy: 



/(sq,w) 



2w — w H — w 
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-w 



37 254 

—w^-\ 

3 7 



2 5 

Of course, this derivation assumes that the limit L ^ oo can validly be interchanged 
with expansion in w. This can presumably be proven with additional work, by methods 
of the cluster expansion [136]. It can presumably also be proven, by the same methods, 
that this small-w expansion is convergent in some disc l^l < e. Note, finally, that we 
have for simplicity used here periodic boundary conditions in both directions, in contrast 
to the cylindrical boundary conditions used elsewhere in this paper; but this change should 
make no difference, because the bulk free energy is expected to be independent of boundary 
conditions. This can be proven rigorously for w > 0, by "soft" methods; and by cluster- 
expansion methods it can presumably be proven also in some disc \w\ < e. 

For the triangular lattice one can perform a similar small-w expansion of the spanning- 
forest partition function. On an L x L triangular lattice with periodic boundary conditions 
in both directions, the result for A; < 5 is 

'3V^ 
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1 + 3Vw + 



3V 
3 



-2V 



w 



+ 



+ 



3V 
4 

3V 



- 3V{2V - 1) 
-3V'{3V-4) 



w 
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+ 0{w') , 



(3.3) 



valid for L > 6. Taking the logarithm, one finds again that each coefficient is proportional 
to V. Dividing by V, one obtains the small-w expansion of the bulk free energy for the 
triangular lattice: 



/(tri,w) 
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(3.4) 
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If one wants to extend these expansions to higher order, direct enumeration becomes 
increasingly tedious and cumbersome. A vastly more efficient approach is the finite-lattice 
method [137]. For the square lattice, one can write the small-iy expansion of the infinite- 
volume free energy as [137, 138] 

/H ^ Jhn^^logZ^-^H = Yl MQ,r)\ogZ^JTrH + 0{w''-') , (3.5) 

where Z^^Kw) is the partition function for a square-lattice grid of size q x r with free 
boundary conditions, and the sum is taken over all rectangles belonging to the set 

B{k) = {{q,r): q < r and q + r < k} . (3.6) 

The weights ak{q,r) are defined as 

(^k{Q,r) = \ Wk{q,q) for g = r (3.7) 

where 



Wk{q,r) 



2Wk{q,r 


) for g < r 


Wk{q,q) 


ioT q = r 





for g > r 
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for g + r = A; 


-3 


ior q + r — k — 1 


3 


for q + r — k — 2 


-1 


foT q + r = k — 3 





otherwise 



(3.8) 



The error term in (3.5) is given by the smallest connected graph with no vertices of order 
or 1 that does not fit into any of the rectangles in B{k) [137]. For the square lattice, this 
graph is a rectangle of perimeter 2k — 2, hence the error term vj^^"'^. The main limiting factor 
on this method is the maximum strip width L^ax we are able to handle (due essentially to 
memory constraints). In particular, if we set the cut-off k to 2Linax + 1, then formula (3.5) 
gives the free-energy series correct through order ^4Lmax-i 25 

The extension of this method to the triangular lattice is straightforward: the free energy 
can be approximated using (3.5), and the weights are given by (3.7)/(3.8) [139]. However, 
the method is less efficient than for the square lattice, as the error term is now larger than 
that of (3.5). If Lmax is the maximum width we can compute and if we set the cut-off 
k = 2Z/niax + 1 as before, then the smallest graphs not fitting into any of the rectangles 
in B{k) are trapezoids that have inclined sides of lengths L^ax and Lmax — 1, respectively, 
together with one vertical side of length 1 and one horizontal side of length 1. The perimeter 
of such graphs is 2L^ax + 1, which implies that (3.5) gives the free energy for the triangular- 
lattice model correct through order ^;2^max 26 ^jj-^jg scries can be improved slightly by noticing 
that the correct term of order ^y^Lmax+i ^.^^^^ obtained by adding the quantity — 2'u;^^™'"'"'"^ 
to the finite-lattice-method result. 



^'^Wc have confirmed this fact empirically by doing runs for different values of Z/max and checking to what 

order they agree. 

^^Once again, we have confirmed this fact empirically by doing runs for different values of Lmax and 
checking to what order they agree. 

^''We discovered this fact empirically by comparing the series obtained for two consecutive values of imax- 
In all cases, we found that the difference was precisely — 2w^^>"i>=+i + 0{w'^^'^''^'^^). It would be interesting 
to find a combinatorial proof of this fact (if indeed it is true in general). 
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We have computed the partition functions Z^^l{w) for both the square and triangular 
lattices, by using two complementary methods based on the transfer-matrix approach (see 
Section 4). We obtained the explicit symbolic form of the transfer matrix for square-lattice 
strips of widths L <9f and for triangular-lattice strips of widths L < Sp; we then computed 
the partition functions Zg^'lj.{w) using the standard formulae [5]. For larger widths (up to 
-^max = 12f), we were unable to compute the transfer matrix explicitly; rather, we used two 
independent programs (written in C and Perl, respectively) to build the partition functions 
layer-by- layer. To handle the large integers that occur in the computation of the partition 
function for large lattices, we used modular arithmetic and the Chinese remainder theorem. 
For the square (resp. triangular) lattice, we obtained the first 47 (resp. 25) terms of the 
free-energy series 

oo 

f{w) = Y,hw''. (3.9) 

k=l 

The results are displayed in Table 2. It is interesting to note that for all the free-energy 
coefficients fk we have computed, the quantity kfk is an integer; it would be interesting 
to understand combinatorially why this is so. From these series expansions one can easily 
obtain the corresponding small-w expansions for the derivatives of the free energy (by differ- 
entiating the series term-by-term with respect to w). The analysis of these series is deferred 
to Appendix A. 



3.2 Large-w; (perturbative) expansion 

On any graph G = (V, E), the generating polynomial of spanning forests can be written 



as 



Fg{w) = w\^\-'J2 



Cfe(G) 



k=l 



W 



k-1 



(3.10) 



where Ck{G) is the number of /c-component spanning forests on G (in particular, ci{G) is the 
number of spanning trees of G) . Taking logarithms, we find 



■^log Fg{w) = ^ logw + -^logci(G') + -^log 



oo 



. (3.11) 



In particular, if we take G to be a large piece of a regular lattice (with a not-too-eccentric 
shape), then in the infinite- volume limit we expect the coefficients of this series to tend 
term-by-term to limits. If this indeed occurs, we obtain a large-w expansion of the form 



/H 



logw + f, + + A + ... . 

w 



(3.12) 



Let us now demonstrate, in the case of the first nontrivial term, that this convergence 
indeed occurs; we shall explicitly compute the limiting value /o, which is the entropy per 
site for spanning trees on G. By the matrix-tree theorem [140, Corollary 6.5], ci{G) equals 
l/|y| times the product of the nonzero eigenvalues of the Laplacian matrix of G. For an 
L X ■ ■ ■ X L simple-hypercubic lattice in d dimensions with periodic boundary conditions in 
all directions, the eigenvectors of the Laplacian are Fourier modes, and the eigenvalues can 
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be written explicitly [134] ; we have 



Ci(L X • • • X L) 



1 

1? 



n 



pe{0,27r/L,...,2(L-l)Tr/L}d 



2^(1 - COS Pi) 



1=1 



Therefore, in the infinite-volume limit we obtain^^ 



lim — logci(L X • • • X L) = /o 



(27r)° 



log 



2^(1 - cospj) 



i=l 



(3.13) 



(3.14) 



(This integral is infrared-convergent in any dimension d > 0.) Analogous formulae hold for 
other regular lattices [114,134,135,142]. For the standard two-dimensional lattices one can 
carry out the integrals explicitly, yielding [114, 134, 135] 
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{6k + 5y 



[V^'(l/6) - ^^'(5/6)] 



127r 

1.615 329 736 097... (3.16) 



where G is Catalan's constant and il^'{z) = {(P /dz^)\ogr(z) is the first derivative of the 
digamma function. See also [142] for high-precision computations of /q on the simple hyper- 
cubic lattice in dimensions 3 < d < 20, as well as on the body-centered hypercubic lattice 
in dimensions 3 and 4. 

We defer to a later paper [143] the computation of /i, /2, . . ., which can be obtained by 
perturbative expansion in a fermionic theory that contains a Gaussian term and a special 
four-fermion coupling (see [24] for details). It is worth mentioning that this perturbation 
expansion can sometimes be infrared-divergent, refiecting the fact that the coefficients of 
(3.11) may in fact diverge as L — > cxo. For example, in dimension d = 1 with periodic 
boundary conditions, we have the exact formula 



L-l 



L-1 1 (L-l)(L-2) _2 
1 H + -w ^ + 



so that 

1 ^ , , L-1 logL L-l 

-logFLp(w) = — ^logw + -J- + 



, (L-l)(L-5) 2 



24L 



(3.17) 



(3.18) 



We see that /i = 1/2, but /2 is infrared-divergent (as are /a and subsequent terms). We do 
not yet know whether f2, fs, ■ ■ ■ are infrared-finite in dimension d — 2. 

^^The same formula holds for the infinite- volume limit of free boundary conditions, but the proof is more 
complicated: see e.g. [141] and the references cited in [141, footnote 1]. 
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We can estimate the first terms fi by computing f{w) numerically at a set of (large) 
values of w, and trying to fit it to a polynomial Ansatz in 1/w: 

Refiw) - log 1^1 - /o = ^ 4 ■ (3-19) 

k=l ^ 

For this calculation, we began by computing the finite-width free energies fhiw) at w = 
±2, ±3, . . . , ±10 on strips up to L = 14p (square) and L = 13p (triangular). We then 
estimated the infinite- volume free energy f{w) by extrapolating the finite-width data using 
the Ansatz (7.35) below (see Sections 7.4 and 7.8 for details). Finally, we used a polynomial 
fit (3.19) with /cmax = 5 to obtain the subleading coefficients fk for k — 1,...,5, fitting 
separately the data with w > and w < 0. More precisely, each sign of w and each value 

of Uiinin = 2, . . . , 6, WC fit thc poiuts with \w\ = "ifmini "W^min + Ij • • • j "^min + 4 tO the Ausatz 

(3.19) with A';jnax = 5. The observed variations in the parameter estimates are due to thc 
neglected higher-order terms k > fcmax in (3.19). These fits are displayed in the rows labelled 
"Non-Biased" in Tables 3 (square lattice) and 4 (triangular lattice). 

It is clear that thc coefficient /i is close to 0.25 = 1/8 for thc square lattice and to 
0.08333 1/12 for the triangular lattice. This motivates the following conjecture: 

Conjecture 3.1 For any regular lattice in dimension d > 2 with coordination number r, 
the large-w series expansion of the spanning-forest free energy takes the form 

f{w) = log|w| + /o + -Li + 0{w-^) (3.20) 

zr w 

where fo gives the entropy per site for spanning trees on the lattice in question. 

Note that the formula (3.20) is false in dimension d — 1: as seen from (3.18), we have 
/i = 1/2, not 1/4. 

We have redone the above fits by fixing fi to its conjectured value l/(2r). Then, for each 
sign of w and each value of Wmin = 2, . . . , 7, we fit the points with \w\ = Wmin, ■ ■ ■ , Wmm + 3 
to this biased Ansatz. The results are shown in rows labelled "Biased" in Tables 3 and 4. 
The next coefficient /2 appears to take the same value for the fits with positive and negative 
values of w: /2(sq) fti 0.0098 and /2(tri) ^ 0.0025. 

Remark. We have very recently proven Conjecture 3.1, and have also understood why 
d = 1 is an exception. These results will be reported elsewhere [143] . 

4 Transfer matrices 

The first step in our analysis is to compute the transfer matrix T and the vectors u and 
Vid that appear in (2.31a). We can proceed in three alternative ways: 

(a) Compute the transfer matrix for the general g-state Potts model, symbolically as a 
polynomial in q and v. (See [5] for the theory, and [14,15] for the computations.) Then 
perform the limit q,v ^ with w = v/q fixed. 

(b) Take the limit q,v ^ with w = v/q fixed right at the beginning, and compute the 
transfer matrix for the spanning-forest problem, symbolically as a polynomial in w. 
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(c) Same as (b), but compute the transfer matrix numerically for a specified (but arbitrary) 
value of w. 



For L < 9 we used methods (a) and (b), and verified that they give the same answer (this is an 
important check on the correctness of our programs). For L = 10 we could use only method 
(b). For L > 11 we were unable to perform the computation symbolically; we therefore 
performed the computations numerically (using machine double- precision) for 11 < L < 16 
at selected real values of w. 

The dimension of the transfer matrix for a square-lattice strip of width L and cylindrical 
boundary conditions is given by SqCyl(L) in Table 1. It coincides with the dimension of 
the transfer matrix for the full Tutte polynomial on square-lattice strips with the same 
boundary conditions [14], and equals the number of equivalence classes modulo translation 
and reflection of non-crossing partitions of {1, 2, . . . , L}. This number is denoted by iVz,sq,PF,L 
in Refs. [14,15]; see (4.2) below for an exact formula. 

The dimension of the transfer matrix for a triangular-lattice strip of width L and cylin- 
drical boundary conditions is given by TriCyl'(L) in Table 1. It coincides with the dimension 
of the transfer matrix for the full Tutte polynomial on triangular-lattice strips with the same 
boundary conditions [15], and equals the number of equivalence classes modulo translation 
of non-crossing partitions of {1,2,..., L}. In Ref. [15] TriCyl'(L) is denoted A^z,tri,PF,L, and 
an exact formula is provided: 



IViCyl'(L) = 1 



L + 1 



d\L,l<d<L 



(4.1) 



where d\L means that d divides L, and 0(n) is the Euler totient function (i.e. the number of 
positive integers k < n that are relatively prime to n). The sequence TriCyl'(L) is given as 
sequence A054357 in [144] ; it equals the number of bi-colored unlabelled plane trees having 
L edges [145]. 

As in previous work on other cases of the Potts model [7, 15], we can reduce the di- 
mension of this transfer matrix down to TriCyl(L) = SqCyl(L) by noting that in the 
translation-invariant subspace of connectivities, the transfer matrix does commute with re- 
flections (see [7, beginning of Section 4] for a detailed explanation). Thus, in a new basis 
with connectivities that are either odd or even under reflection, the transfer matrix takes 
a block-diagonal form: the block that is even under reflection has dimension TriCyl(L) = 
SqCyl(L) and a non-vanishing amplitude; the block that is odd under reflection has a zero 
amplitude and hence docs not contribute to the partition function. 

The dimensions SqCyl = TriCyl and TriCyl' are related by the following equation (proved 
in [15]): 



SqCyl(L)=TriCyl(L) = l 



TriCyl'(L) + 



(lv2j) 



(4.2) 



The numerical values of SqCyl(L) and TriCyl' (L) up to L = 20 are shown in Table 1. Both 
grow asymptotically like 4^L~^/^ as L — > oo. 
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5 Square-lattice strips with cylindrical boundary con- 
ditions 



For each lattice width L up to L = 10, we have computed symbohcally the transfer 
matrix T and the vectors u and Vij. Then we computed the zeros of the partition function 
FmpxriFiuj) for strips of aspect ratio p = 5 and p = 10; and for L < 8 we also computed the 
accumulation set B of partition-function zeros in the limit p — > oo. For L = 9, 10 we were 
unable to compute the full limiting curve B, but we did compute some selected points along 
it. 

The limiting curves B resulting from these computations are shown in Figures 3-6 (su- 
perposed in Figure 7), and their principal characteristics are summarized in Table 5. One 
interesting feature of B is the point (s) where it crosses the real w-axis. For odd width L, it 
turns out that there is only one such point, which we shall denote wqq{L). For even width L, 
by contrast, B contains a real segment [wq^{L), Wo+(-^)]- This segment contains a multiple 
point of i3, where an arc of B lying at \m.w ^ crosses the real axis; we shall denote this 
point Wqq{L). It is a curious fact that for the even- width square lattices considered here 
(L = 2,4,6,8,10), we find u'o+(I/) = —1/4 exactly. Finally, we define wb{L) to be the 
complex-conjugate pair of endpoints of B with the largest real part. 

For L < 5 we computed the limiting curve B using the resultant method explained 
in [5, Section 4.1.1]. For L = 6, we computed the endpoints using this method, and the 
rest of B using the direct-search method explained in [5, Section 4.1.2]. In all these cases 
we can be sure that we did not miss any endpoints or connected components of B. For 
L > 7, hy contrast, we were unable to compute the resultant; we therefore located B using 
the direct-search method. Here we could easily have missed some small components of B. 
Our reports in Table 5 for the number of connected components (t^C), endpoints (t^E) and 
multiple points (#Q) must therefore be viewed for L > 7 as lower bounds. For L > 9 the 
computation of the limiting curve B is very time-consuming; we therefore computed only a 
few important points. 



5.1 L = 2 

In this case the connectivity basis is two-dimensional. The transfer matrix T and the 
vectors Vid and u are given by 



(5.1a) 



2w + l Bw"^ + Aw + 1 
= (l,2w + l) (5.1b) 
= (0,1) (5.1c) 



The zeros of the polynomials F2pxnF{w) with n = 10, 20 are displayed in Figure 3(a). In the 
same figure we also show the limiting curve B. The curve B is connected: it is the union of 
a horizontal segment running from wq- — —1/2 to Wo+ = — 1/4 and an arc running between 
the complex-conjugate endpoints wb — —1/4 ± (-\/3/4)i. The segment and the arc cross at 
the multiple point at wqq ^ —0.3660254038. 

It is worth noting that, even though the curve B contains the point w = —1/2 as an 
endpoint, the value w = —1/2 is not a zero of the partition function for any finite strip 
2p X np other than the trivial case n — 1. As a matter of fact, the partition function takes 
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the value 

i^2pxnp(«^ = -1/2) = -{n - l)/4"-i . (5.2) 

This observation supports the conjecture that there are no roots in the half-plane Rew < 
—1/2, i.e. that the family 2p x np of strip graphs possesses the univariate dual Brown- 
Colbourn property for all finite lengths n. 



5.2 L = 3 

The connectivity basis is three-dimensional; the transfer matrix T and the vectors Vid 
and u are given by 

T = w^{7w + 1) w^{Aw + 1) (5.3a) 

\ 3tt; + 1 3(8w2 + 5w + l) 16w^ + 15w^ + 6w + 1 J 

= (l,3(2w + l),3w^ + 3w + l) (5.3b) 

= (0,0,1) (5.3c) 

The zeros of the polynomials -F3pxmp('"^) with n — 15, 30 are displayed in Figure 3(b), along 
with the limiting curve B. The curve B has three connected components. One of them nms 
between the complex-conjugate endpoints —0.2432796623 ± 0. 1389560739 i and intersects 
the real w-axis at wqq ^ -0.2868497019. The other two run from w a; -0.3020539886 ± 
0.1587843374 i to wb ^ -0.0574443351 ± 0.48062531611 



5.3 L = 4 



The connectivity basis is six-dimensional; the transfer matrix T and the vectors Vjd and 



u are given by 



T = 



w 






4/^11,6 



4w^ 

4 

W 

2/^12,6 



12^6 



u 



Vid = 



\ 4:W +1 

(1,4C2,2C2,4L>3,3,2C2^C2L'2,2) 

(0,0,0,0,0,1) 



where we have used the shorthand notations 



kw + 1 



nw 



1 



8w;° 
2w4C5 


2C3D10 5 



2w^C3 

26,6 



45w^ + 52w^ + 28w^ + 8«; + 1 



(5.4a) 



(5.4b) 
(5.4c) 



(5.5a) 
(5.5b) 
(5.5c) 



The zeros of the polynomials Fi^xripiw) with n — 20, 40 are displayed in Figure 3(c), along 
with the limiting curve B. The curve B has three connected components. One of them is the 
union of a horizontal segment running from ^ —0.2776248333 to ^0+ = —1/4 and an 
arc running between the complex-conjugate endpoints w ~ —0.2252016334 ±0.2287233280 i. 
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The segment and the arc cross at the multiple point at wqq ^ —0.2670015604. The other two 
components arc complex-conjugate arcs running from w ~ — 0.2068057238±0. 2160687519 i to 
wb ~ 0. 0632994010±0. 5099839130 i. Finally, there is a pair of very small complex- conjugate 
bulb-like regions emerging from the endpoint Wq- ~ —0.2776248333 — which is therefore a 
T point — and going back to the real w-axis at the multiple point w ^ —0.2775806860: see 
the blow-up picture in Figure 4. 

5.4 L = 5, 6, 7, 8, 9, 10 

The transfer matrices for L > 5 are too lengthy to be quoted here. Those for L < 9 can 
be found in the Mathematica file f orests_sq_2-9P .m that is available with the electronic 
version of this paper in the cond-mat archive at arXiv.org. The file for L — 10, which is 
13.6 MB long, can be obtained on request from the authors. 

We have plotted for each L (= 5,6,7,8,9, 10) the zeros of Flpx{pL)f{w) for aspect ra- 
tios p = 5, 10 as well as the hmiting curves B {p — oo). See Figure 3(d) for L — 5, 
Figures 5(a,b,c,d) for L = 6, 7, 8, 9, and Figure 6 for L = 10. 

The principal features of the limiting curves B arc summarized in Table 5. For L = 5, 
there are ten endpoints located &t w ^ -0.2645566722 ± 0.0812460198 i, -0.2562049781 ± 
0.0804259250 -0.1745330113 ± 0.2687656693 i, -0.1684195011 ± 0.2629234881 i, and 
0.1533657968 ± 0.5306112949 i (the latter is wb)- The hmiting curve crosses the real axis at 
woQ ^ -0.2620754678. 

For L — 6 there is a multiple point at wqq ^ —0.2563792782. This point belongs to 
a horizontal segment running from ^ —0.2609570768 to Wo+ = —1/4. There are ten 
more endpoints at w « -0.2441302665 ± 0.1392115396 i, -0.2401194057 ± 0.1380752982 i, 
-0.1359169474 ± 0.2964285650 i, -0.1338793838 ± 0.2939739183 i, and 
0.2262944917 ± 0.5460127254 i (the latter is wb). 

For L = 7 there are fourteen endpoints located at w ^ -0.258797±0.056312 i, -0.256989± 
0.056335i, -0.222450±0.180872i, -0.220899±0.18026H, -0.103677±0.316136i, -0.104346± 
0.317122 i, and 0.2879810252 ± 0.5579895735 i (the latter is wb). The limiting curve crosses 
the real w-axis at woq ^ -0.2559077691.^9 

For L = 8, we find a multiple point at wqq ^ —0.2532620078, which belongs to a 
horizontal segment running from wq- ~ —0.2556827597 to Wo+ = —1/4. There arc fourteen 
additional endpoints at w ^ -0.250381 ± 0.101684i, -0.249500 ± 0.101496i, -0.201140 ± 
0.212546 i, -0.200485 ± 0.212232 i, -0.077462 ± 0.333326 i, -0.077248 ± 0.332940 i, and 
0.3415981731 ± 0.5675966263 i (the latter is wb). 

For L = 9 we were unable to compute the full limiting curve B, but we were able to 
compute the points corresponding to fixed values of | Im^l = 0, 0.01, 0.02, etc. These points 
give us a rough idea of the shape of the limiting curve. In carrying out this computation, it 
is important to keep track of the quantity 9 = \ Imlog(A*/A^)| as we move along the limiting 
curve B (here A* and A^ denote the dominant equimodular eigenvalues), since 9 vanishes at 

^^As explained at the beginning of this seetion, we loeated the endpoints nsing the resuhant method for 
L < 6 and the direet-seareh method for L > 7. The direet-search method is quite efficient for locating 
a real endpoint, but is extremely tedious for locating a complex endpoint (since we have to search a two- 
dimensional space). This explains why the precision obtained for the complex endpoints with L > 7 (error 
« 10~^) is inferior to that obtained both for the endpoints with L < Q and for the real endpoints with L >7 
(error « 10~^°). However, we have put an extra effort in computing the rightmost complex endpoints Wb 
for L > 7 with error w 10~^° (see Section 7). 
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endpoints [5]. By careful monitoring of the angle 6, we can obtain a lower bound on the 
number of endpoints and connected components of the curve B. In particular, for L = 9 we 
have found that there are at least 18 endpoints and 9 connected components (see Table 5). 
We have also computed the point where the the limiting curve crosses the real w-axis Wqq ~ 
-0.2534832041, and the rightmost endpoints wb ^ 0.3890914478 ± 0.57549594941 

For L = 10 the computation of the limiting curve B using arbitrary-precision Mathe- 
MATICA scripts is beyond our computer facilities. However, we have been able to compute 
some points along this curve (those corresponding to | Im u"! = 0, 0.01, 0.02, etc.) by using the 
double-precision Fortran subroutines of the arpack package [146].^° We find a multiple point 
at Wqq ~ —0.2519570283. This point belongs to a horizontal segment running from wq- ~ 
—0.2534268314 to 1^0+ = —1/4. Using the same method as in L = 9 we have concluded that 
B is formed by at least 9 connected components and contains at least 20 endpoints. Finally, 
we have computed the rightmost endpoints wb ~ 0.4317571213 ± 0.5821234989 i. 

6 Triangular-lattice strips with cylindrical boundary 
conditions 

We performed for the triangular lattice the same calculations as for the square lattice. 
Thus, for each strip width L < 10 we computed the transfer matrix and the associated left 
and right vectors; from these we obtained the partition-function zeros for strips with aspect 
ratio p = 5. 10 as well as the limiting curves B corresponding to the limit p — > 00. The 
limiting curves B resulting from these computations are shown in Figures 8-10 (superposed 
in Figure 11), and their principal features are summarized in Table 6. Once again, the full 
curves B are computed by the resultant method for L < 5, by a combination of the resultant 
method (endpoints only) and direct search for L — 6, and by the direct-search method for 
L>7. 



The connectivity basis is two-dimensional, and the transfer matrix T, and the vectors Vid 
and u are given by 



We display in Figure 8(a) the zeros of the polynomials -F2pxnF('"^) ^ — 10, 20, along with 
the corresponding limiting curve B. The curve B is connected: it is the union of a horizontal 
segment nmning from Wq_ = —1/4 to Wo+ ~ —0.2017782928 and an arc running between 
the complex- conjugate endpoints wb ~ —0.1178608536 ± 0.2520819223 i. (The latter three 
points are the roots of the polynomial 1 -|- 8w -|- 28w^ + 6Aw^.) The segment and the arc 
cross at the multiple point at wqq ^ —0.2251972448. 

■^"The arithmetic precision of the arpack package is thus less than that of Mathematica. However, we 
have checked in a difficult but manageable case {L = 9) that the results obtained by the two methods agree 

to at least ten decimal digits. We have also checked the performance of the ARPACK routines at some specific 
points {wo+, Wo- and wqq) for L = 10, and the disagreement with Mathematica is again less than 10^^°. 



6.1 L = 2 




f Aw^ 2w'^(6w + 1) 
\4:W + 1 12^2 + 6w + l 

{l,2w + l) 

(0,1) 



(6.1b) 
(6.1c) 



(6.1a) 



32 



6.2 L = 3 



The connectivity basis is three-dimensional, and the transfer matrix T, and the vectors 
Vid and u are given by 

/ 8w^ 6w^{12w + l) 9w\6w + l) \ 

T = 4«;2 2w^{19w + 3) w^{30w^ + lOw + 1) (6.2a) 
\ 6W + 1 3{20w^ + 8w + l) 50w^ + 33w^ + 9w + 1 J 

= {l,3{2w + l),3w^ + 3w + l) (6.2b) 

= (0,0,1) (6.2c) 

We display in Figure 8(b) the zeros of the polynomials -F3pxnF('"^) with n — 15, 30, along with 
the corresponding limiting curve B. The curve B has three connected components. One of 
them runs between the complex-conjugate endpoints w ^ —0.1763241163 ± 0.1013423678 i 
and intersects the real w-axis at Wqq ~ —0.1921127637. The other two run from w ~ 
-0.1662139535 ± 0.0960130540 i to wb ^ -0.0100833905 ± 0.2849353892 i. 



6.3 L = 4 



The connectivity basis is six-dimensional; the transfer matrix T, and the vectors vjd and 
u are given by 



(6.3a) 



u 



/ 16w^ 


16w^Ci2 


96^5 




2w^Cj^ 


36w^C6 


\ 


Sw^ 


12w^C8 


2w^C24 




2w^CgCi2 


3U'^£)36,11 









8^4 




3w^Cs 






Aw^ 


2u;2T4,2 




2tv^n^^ 












2w'^Ci2 


2w'^D52M 


2w^Ci 


2w^D24^9 




\ Cs 


4-^28,10 


2-D32,io 


^6,4 


20^024^8 


^6,6 


/ 


(1,4^2, 


2C2,4Z}3,3 


,2ClC2D 


2,2) 









V: 



id 



(0,0,0,0,0,1) 



where we have used the shorthand notations (5.5a)/(5.5b) and 



T2,4(W 
T4,2{W 

T4a{w 

T4,5iw 
T4fi{w 

Te,4iw 

The zeros of the polynomials F4 



{8w + l){27w + 2) 
28w + 5 
80^2 + 28w; + 3 
{6w+ l)(16w + 3) 
w^{96w^ + 5Aw^ + 12w 
A{80w^ + A7w^ + llw + 
192w^ + 164w^ + 62w^ + 12w + 1 
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(6.3b) 
(6.3c) 

(6.4a) 
(6.4b) 
(6.4c) 
(6.4d) 
(6.4e) 
(6.4f) 
(6.4g) 



P xnp 



[w) with n = 20, 40 are displayed in Figure 8(c), along 
with the limiting curve B. The curve B has three connected components separated by two 
very small gaps. One of the components is the union of a horizontal segment running from 
Wo- ~ —0.1846154722 to ^0+ = —0.1833753245 and an arc running between the complex- 
conjugate endpoints w fa —0.1343195918 ± 0.1422271304 i. The segment and the arc cross 
at the multiple point at wqq ^ —0.1839945026. The other two components are complex- 
conjugate arcs running from w a; -0.1325113192 ± 0.1408071373i to wb ^ 0.0588726934 ± 
0.3024953798 i 
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6.4 L = 5, 6, 7, 8, 9, 10 



The transfer matrices for L > 5 are too lengthy to be quoted here. Those for L < 9 can 
be found in the Mathematica file f orests_tri_2-9P .m that is available with the electronic 
version of this paper in the cond-mat archive at arXiv.org. The file for L — 10, which is 
31.9 MB long, can be obtained on request from the authors. 

We have plotted for each L {— 5,6,7,8) the zeros of -Flpx(pL)f('"^) for aspect ratios 
p = 5, 10 as well as the limiting curves B {p — oo). See Figure 8(d) for L — 5, Figure 9(a) 
for L = 6, Figure 9(b) for L = 7, and Figure 9(c) for L = 8. 

The principal features of the limiting curves B are summarized in Tabic 6. For L = 5, 
there are ten endpoints located at w ^ -0.1732373428 ± 0.0559213543 i -0.1725608663 ± 
0.055734331H, -0.1046922573 ± 0.1666741778 i, -0.1043151240 ± 0.1662893120 i, and 
0.1110085784 ± 0.3142261926 i (the latter is wb)- The limiting curve crosses the real axis at 
woQ^ -- 0.1805863920. 

For L = 6 there is a multiple point at wqq ~ —0.1788458357. This point belongs to a 
very short horizontal segment running from wq- ~ —0.1788596197 to Wo+ ~ —0.1788320527. 
There are ten more endpoints at w ^ -0.1572400315 ± 0.0919540390 i, -0.1570191889 ± 
0.0918929205 i, -0.0812379488 ± 0.1829264297 i, -0.0813151207 ± 0.1830202393 i, and 
0.1534513521 ± 0.3227719569 i (the latter is wb). 

For L = 7 there are fourteen endpoints located at -0.173950 ± 0.039380?, -0.173885 ± 
0.039369 i, -0.141195±0. 117208 i, -0.141154±0.117209 i, -0.061954±0.194837i, -0.061949± 
0.194832 i, and 0.1894308555 ± 0.3293185497i (the latter is wb). The limiting curve crosses 
the real axis at Wqq ~ —0.1778368691. 

For L — S there is a multiple point at Wqq ^ —0.1772011941. This point belongs to a 
very short horizontal segment running from wq- ~ —0.1772035681 to Wo+ ~ —0.1771988202. 
There are fourteen additional endpoints at -0.165087 ± 0.068518 i, -0.165105 ± 0.06852H, 
-0.126193±0.13589H, -0.126182±0.135887i -0.045467±0.203914i, -0.045470±0.203918i, 
0.2207308275 ± 0.3345156821 i (the latter is wb). 

For L = 9 we were unable to compute the full limiting curve B, but as in the square- 
lattice case, we were able to obtain a rough estimate by computing the points corresponding 
at fixed values of |lmw| = 0,0.01,0.02, etc. In particular, the point where B crosses the 
real axis is Wqq ~ —0.1767753501. We have also computed the rightmost endpoints Wb ~ 
0.2484623416 ± 0.3387563998 i. Using the same method as in the square-lattice case, we 
have concluded that B is formed by at least 9 connected components and contains at least 
18 endpoints (see Table 6). 

For L = 10, we have computed the points of the limiting curve corresponding to | Imtfl = 
0, 0.01, 0.2, etc by using the ARPACK subroutines. We find a multiple point at Wqq ~ 
—0.1764764073. This point belongs to the horizontal segment running from wq- ~ —0.1764769384 
to 1^0+ ~ —0.1764758762. Using the same method as for L = 9, we conclude that the limit- 
ing curve is formed by at least 9 connected components and contains at least 20 endpoints. 
Finally, we have computed the rightmost endpoints wb ~ 0.2733731381 ± 0.34229360971 

7 Extrapolation to infinite width 

In this section we analyze the finite-strip data from Sections 5 and 6 and study the 
limiting behavior as the strip width L tends to infinity. Our goal is to determine the nature 
of the phase transition at w — wq and to extract numerical estimates of the critical point wq, 
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the free energy f{w), the central charge c{w), and the thermal scaling dimensions XT,i{w) 
for each of the two lattices. 

7.1 Procedure I: Estimates of critical points 

One goal of this paper is to obtain the limiting value as L — > oo of the quantities Wo+{L), 
wo-{L) and wqq{L). We expect that all three quantities converge to the same limit wq. 

Wo = lim woj{L) for j = +,-,Q . (7.1) 

This limit is also expected to be independent of the choice of boundary conditions (e.g. free, 
cylindrical, cyclic, toroidal, etc.). We have estimated Wq by fitting our finite- width data to 
the Ansatz 

WojiL) = wo + AjL-^ . (7.2) 

At least for the square lattice, this Ansatz is theoretically justified due to the vicinity of the 
critical theory at = 1/4, and we should have A = 1/z/. Thus, for each value Lmin we fit the 
data for the triplet L — Lmm, -Linin-l-2, Linin-|-4 to the Ansatz (7.2) and extract the estimates 
of Wo, Aj and A. (For j = Q we can also use the triplets L = Lmin, Lj^in + 1, L^in + 2.) The 
observed variations in the estimates as a function of Lmin arise from higher-order correction- 
to-scaling terms that are neglected in (7.2). 

It is also of interest to estimate the exponent uj associated to the width of the interval 
[wo-,iWo+] for even L: 

wo+{L) - wo-{L) ^ AL-"" . (7.3) 

We fit to the Ansatz (7.3) using pairs L — Lmin, L^in + 2. 

By definition, the endpoints wo±{L) correspond, at any given finite L, to the collision (in 
modulus) of the two dominant eigenvalues. Now, for even widths L on the square lattice, we 
have found that Wq^{L) = 1/4 exactly. Therefore, at w = Wq = -1/4 we have ^{^{L) = 
exactly for all even L (and not merely in the limit L — * oo). From (2.37) we can conclude 
that the leading thermal scaling dimension xt,i is zero, and hence that u — 1/2. We expect 
this latter conclusion to hold (by universality) also for the triangular lattice. In view of this 
prediction, we can expect to obtain more accurate estimates for wq by fixing A in the Ansatz 
(7.2) to its predicted value A = l/u = 2. 

More generally, we have also tried to extrapolate the whole limiting curve B to L ^ oo: 
for fixed Imw — 0.01, 0.02, . . ., we extrapolated the values of Rew using an Ansatz of the 
form (7.2). We performed several different extrapolations, e.g. including all data points or 
including only even or odd widths. We obtained consistent estimates for | Im-u;| < 0.33 on 
the square lattice and | Im-U7| < 0.23 on the triangular lattice (see Table 7). The results are 
given by the black dots of Figures 7 and 11. 

By the same method we can obtain the limit as L — > oo of the endpoints with largest 
real part, which we have denoted wb{L). As a first guess, we used the Ansatz 

wb{L) = WB + AL~^ (7.4) 

separately for the real and imaginary parts of wb{L). Then, after examining the results, we 
conjectured better Ansatze (see Sections 7.3 and 7.7 below). 



35 



7.2 Procedure II: Estimates of the free energy 

The second goal of this paper is to estimate the bulk free energy 



lim fL{w) 



(7.5) 



and its derivatives as a function of w. For each of the two lattices (square and triangular), 
we have computed the five largest eigenvalues (in modulus) of the transfer matrix T{w) for 
selected real values of w and for all L < 14. The w values are taken between w — —2 and 
w — 2 in steps of 0.1. Prom the dominant eigenvalue Xi,{w) we can extract the strip free 
energy 



[cf. (2.32)]. In Figures 12 and 21 we have plotted the real part of the strip free energy, 



for the square and triangular lattices, respectively. The imaginary part of the free energy, 
not shown in the plot, arises from the following qualitative properties of the eigenvalues: 

(a) For odd L: 

— For w > woQ, there is a unique dominant eigenvalue X*{w), which is real and 
positive. 

— At w — WoQ, the dominant eigenvalue Xi,{w) > becomes equimodular with an 
eigenvalue — A*(w) < 0. 

— For w < wqq, there is a unique dominant eigenvalue Xi,{w), which is real and 
negative. 

(b) For even L: 

— For w > wo+, there is a unique dominant eigenvalue A*(w), which is real and 
positive. 

— At w — wo+, the two leading eigenvalues collide. The common value A*(w) is still 
real and positive. 

— For wq- < w < ^0+, there is a complex-conjugate pair of dominant eigenvalues. 

— At w — Wo-, the two leading eigenvalues again collide. The common value Xi,{w) 
is real and positive. 

— For w < Wo-, there is a unique dominant eigenvalue Xi,{w), which is real and 
positive. 

The only exception to the above pattern (at least for L < 8) is the square-lattice strip with 
L — A at w — Wo- (see Figure 4). In this case, the point Wq- happens to be a T point, 
not a regular endpoint. At wo-, three eigenvalues attain the same modulus: one is real and 
positive, while the other two form a complex-conjugate pair. 

The limit L — oo of the finite- width free energy fiiw) can be extracted by using the 
Ansatz (2.34) or (2.33), depending on whether the model is critical or not. As we do not 



fhiw) = -\ogX^{w) 



(7.6) 




(7.7) 
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know a priori to which regime a given value of w belongs (this is especially true for the 
triangular- lattice model, for which we do not know the exact value of Wq), we can try first 
the preliminary Ansatz 

fdw) = fiw) + A(«;)L-^(-) . (7.8) 

We fit to triplets L = Lmin, -^^min + 2, -Lmin + 4 in order to minimize the effects of even-odd 
oscillations. The variations of the estimates for f{w), A{w) and A(w) as a function of 
will give us an idea about the critical or non-critical nature of the model. If the model is 
critical, we expect that the estimates for A{w) will converge to 2 as Lmin — > oo; if the model 
is non-critical, we expect that the estimates for A{w) will grow without bound. Once we 
determine the regime to which a given value of w belongs, we can use the more appropriate 
Ansatz (2.34) or (2.33) to obtain more accurate estimates of the physical parameters, notably 
the central charge c{w) for the critical models. 

Finally, by studying the L-dependence of the gap between the dominant eigenvalue A*(u') 
and the subdominant eigenvalues Xi{w), we can estimate the thermal scaling dimensions XT,i 
using the Ansatz (2.37). 

7.3 Square lattice: Estimates of critical points and phase bound- 
ary 

The square-lattice data that we wish to extrapolate arc collected in Table 5 and depicted 
graphically in Figure 7. The quantities Wq_{L) and Wq^{L) arc defined only for even L. The 
quantity Wqq{L) is defined for all L, but it exhibits such strong even-odd oscillations that 
it makes sense to fit it separately for odd and even L. We therefore perform the fits to the 
Ansatz (7.2) using L — Lmin, -^min + 2, Lmin + 4. The results are displayed in Table 8. The 
estimates of wq are well converged, and from Lmin = 5, 6 we can estimate 

u;o(sq) = -0.2501 ± 0.0002 . (7.9) 

This result agrees very well with the theoretical prediction (1.11): 

wo(sq) = -1/4. (7.10) 

It is also noteworthy that for all even values of L < 10, we have found that Wq^{L) = — 1/4 
exactly. Therefore, the exponent u) defined by (7.3) is identical to the exponent A extracted 
from wq_{L). 

The exponent A in Table 8 seems to be approaching A. — 1/v ^ 2 [aX least roughly). 
This behavior is consistent with the interpretation that at w = there is a first-order 
phase transition (see Section 7.4 for overwhelming evidence that the first derivative of the 
free energy is in fact discontinuous at w = Wq): such a phase transition can be described in the 
renormalization-group framework by a discontinuity fixed point characterized by u — 1/d, 
where d is the dimensionahty of the lattice (d — 2 in our case) [94, 147]. The behavior 
A = 2 is also consistent with an ordinary critical point having xt = (see Section 7.6 for an 
analysis of the correlation length in this model). Note, finally, that the only possible values 
of A consistent with conventional finite-size-scaling theory are A < d (corresponding to a 
critical point or a first-order transition) and A = oo (corresponding to a non-critical point). 
Estimates of A that are slightly larger than 2 but decreasing with L (or, at any rate, not 
strongly increasing with L) therefore suggest that A = 2. 
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Numerical investigations of first-order pfiase transitions [148, and references therein] sug- 
gest that the higher-order corrections to the L~'^ behavior are simply integer powers of L~'^. 
This observation motivates the Ansatz 

woj{L) = Wo + A2L-^ + A^L-^ (7.11) 

The results of fitting the data to this Ansatz are shown in Table 9. In this table we observe 
that the estimates for wo(sq) converge slightly faster to values close to the exact result than 
do the estimates in Table 8 (at least for small L). 

Let us next extrapolate the endpoints Wb{L), handling separately their real and imag- 
inary parts. The real part of Wb{L) seems to diverge as L — ^ oo: indeed, as we increase 
L-mim the estimate for Rew^ using the Ansatz (7.4) appears to grow without bound (e.g., 
wb ~ 1-19 for Lmin = 2, wb ~ 2.75 for L^i^ — 3, and wb ~ 9.18 for L^i^ — 4). We therefore 
tried to guess the behavior of Rews(L) as a function of L, and found that it fits very closely 
to the Ansatz 

RewB(L) = CilogL + ca . (7.12) 

The estimates for Ci and function of Lmin are shown in Table 10(a). For Lmin > 5, 

the estimate of Ci is clearly rising with Lmin and is not yet slowing down. This suggests that 
the true value of Ci is quite a bit larger than the value ~ 0.405 observed at our largest Lmin; 
how much larger is difficult to say, especially since it is very difficult to fit a slowly- varying 
function (such as logL) with data points in the (extremely narrow) interval 2 < L < 10. 
One expects, in any case, that the true Ci may exceed the value 0.405 by an amount 
considerably larger than the observed variations in Table 10(a). If forced to guess, we might 
estimate 

Cl = 0.41 ±0.03 (7.13a) 
C2 = -0.52 ±0.06 (7.13b) 

but these estimates should be taken with a grain of salt! 

A more precise Ansatz can be tried if we recall the recently-discovered relationship [24] of 
this model with the A^- vector model with = — 1 and a sign change in the coupling constant. 
This relation implies that the spanning-forest model is perturbatively asymptotically free, 
and yields the theoretical prediction 

RewB{L) = CilogL + CgloglogL + C2 + ^( ^^ogL^ ) ^^'^^^ 

with Cl = 3/(27r) ^ 0.477465 and cg = -l/(27r) ^ -0.159155 [24].^^ The estimates obtained 
by performing a one-parameter fit to (7.14) with Ci and C3 fixed to their theoretical values 
(and the correction term O (log log L/ log L) neglected) are displayed in the last column of 
Table 10(a). Prom these values we can estimate that 

C2 = -0.535 ±0.005 (7.15) 

Please note, however, that the correction terms omitted in (7.14) — loglogL/logL, 1/logL 
and so forth — are extremely slowly varying functions of L. Hence, as we have access only 

^^Indeed, it was our numerical observation of the behavior (7.12) that led us to conjecture that the 
spanning-forest model is asymptotically free in two dimensions — a conjecture that played an important 
role in catalyzing the work reported in [24] . 
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to a very narrow range of L values, we cannot rule out the possibility that the actual value 
of C2 may differ from the above estimate by many times our estimated error. 

The imaginary part of Wb{L), by contrast, seems to converge to a finite limit: the 
estimates using the Ansatz (7.4) are displayed in Table 10(b). We conclude that Imw^ = 
0.69 ± 0.03 and A 1/2. Guessing that A — 1/2 exactly, we can refine the estimate for 
Im wb by fitting the data to the Ansatz 

Imws(L) = ImwB + AL'^^^ (7.16) 

The results arc displayed in Table 10(b). The estimate for lmwB{sq) is still decreasing with 
Lmin, and our best estimate is 



Imws(sq) = 0.70 ±0.02. (7.17) 

Finally, let us try to extrapolate the limiting curve B to L oo; the result is a curve 
that can be interpreted as the phase boundary in the complex w-plane. We began 
by extrapolating Rew at the selected values Imw — 0.01,0.02, . . ., using the Ansatz (7.4) 
for Rew. We performed several different extrapolations, e.g. including all data points or 
including only even or odd widths. We obtained consistent estimates for | Imwl ^ 0.33 (see 
Table 7). The results are given by the black dots of Figure 7. 

Next we attempted to complete this curve at larger Rew. The behaviors (7.12)-(7.15) and 
(7.16)-(7.17) suggest that the curve B^ continues all the way to Rew = +oo, asymptotically 
approaching Imw ^ ±0.70 as Rew — > +oo. This behavior is in agreement with the theoret- 
ical prediction arising from renormalization-group computations in the four-fermion model 
describing spanning forests at small g = 1/w [24]. These computations show [24] that the 
fixed point at = (a free fermion theory) is marginally repulsive for g > (i.e., perturba- 
tively asymptotically free) and marginally attractive for g' < 0, with a renormalization-group 
fiow given by 

+ hg^ + hg' + hg' + ... (7.18) 



dl 

or equivalently 

dw 63 64 bs 

— = -62 ^ 3 - ■■■ 

with 62 > 0, where / is the logarithm of the length rescaling factor. The phase boundary B^o 
is a special RG fiow curve: namely, it is a sepamtrix that divides the complex ?«-planc into 
two "phases" . Initial conditions belonging to the "phase" containing w = are attracted to 
the "high-temperature" fixed point at w = and are therefore noncritical. Initial conditions 
belonging to the "phase" containing w — —00 are attracted to the free- fermion fixed point 
at Rew = — cxo and arc thus critical. From (7.19) we see that all the RG fiow curves (not 
only the separatrix) tend to a constant value Imw as Rew ^ +00. It is a nonperturbative 
question to determine which one of these RG fiow curves is the separatrix. But our transfer- 
matrix calculations yield an approximate answer to that question — namely, Imw ~ ±0.70 
— which we can combine with the perturbative calculations to estimate the shape of the 
curve Boo at large Rew. To do this, let us write w = a + ip, substitute into (7.19), separate 
real and imaginary parts, and divide to obtain a differential equation for the unparametrized 
fiow curves: 

«P _ »3a2+/32 + »4(a2+^2)2 -t- 05(a2+/32)3 • • • ^ 7 90^1 

rim h h a L a^-jS^ l a^-Sa/S^ \' ■■^^) 

aa - &3^^2^ - 04 (^2+^2)2 - 05 („2+^2)3 - • • • 
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Plugging into (7.20) the Ansatz 



Ai A2 As 



we obtain 



Ai = 63/62 
A2 — 64/62 

As = h/b2 - I3^h/{3b2 



(7.21) 



(7.22a) 
(7.22b) 

(7.22c) 



Finally, we need to determine the RG coefficients 62, 63, . . . . This can be done by applying 
the recently-discovered mapping [24] of the spanning-forest model onto the A^- vector model 
with = — 1 and a sign change in the coupling constant (i.e., Wforests = —/^w- vector)- Using 
the known coefficients of the A"- vector RG beta function through four loops [149-151], we 
obtain 

(7.23a) 
(7.23b) 
(7.23c) 
(7.23d) 
and hence 



62 = 


-Wo 


N=-l 


= 3/(27r) 


bs = 




N=-l 


= -3/(27r)^ 


64 = 




N=-l 


^ 2.34278457/(27r)^ 


b5 = 




N=-l 


^ 1.43677/(27r)^ 



A2 
A3 



0.78092819/(27r)2 
6.92 706/(27r)^ if /?o = 0.7 

For numerical purposes, we found it convenient to use the variant Ansatz 

-81,-62,-63, B4 



+ 



+ 



+ 



a — {a — CKo)^ (o; — ccq)^ (ck — ckq) 



+ 



(7.24a) 
(7.24b) 
(7.24c) 



(7.25) 



If we impose = B4 — and fit to the theoretical prediction (7.24) of the first three 
derivatives at a = 00, we obtain ao = -0.479224, Bi = -l/(27r), B2 = -0.069155; the 
resulting fit (7.25) evaluated at a = —0.0674 has value 0.316354 and derivative 0.923329, 
compared to the correct value 0.33 and derivative 0.6849 deduced from Table 7. This is 
quite good agreement between extrapolated perturbation theory and our nonperturbative 
estimates. We therefore used the Ansatz (7.25), imposing the first three derivatives at 
a — 00 along with the estimated value and derivative at a = —0.0674; we obtain 



"0 

B, 
B2 
Bs 
Ba 



-1.06145 
-l/(27r) ^ 
-0.161818 
-0.134478 
-0.284334 



-0.159155 



(7.26a) 
(7.26b) 
(7.26c) 
(7.26d) 
(7.26e) 
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The corresponding curve (7.25) is depicted as a black dotted-dashed curve in Figure 7. 
This curve is presumably not quite right, because it intersects the L = 10 curve, whereas the 
monotonicity of the finite-L curves suggests that the curve B^o lies to the right of all of them. 
Therefore, the true curve B^o probably lies slightly to the right of the dotted-dashed curve 
in Figure 7. Nevertheless, we suspect that this latter curve is a fairly good approximation 
to the true curve Boo, especially at large positive Rew. 

7.4 Square lattice: Behavior of the free energy and its derivatives 
as w ^ wq 

We have computed the real part of the strip free energy. Re fL{w) = (1/L) log |A*(w)|, 
for widths 2 < L < 14 and w — —2, —1.9, . . . , 1.9, 2 as well as w = — 1/4. These results are 
plotted in Figure 12. We have extrapolated these values to L = oo (sohd black hue) using 
the Ansatz 

Rehiw) = Re/H + A{w)L-^^'^^ . (7.27) 

We show also the [20,20] Fade approximant to the small-w series for this quantity (dashed 
violet curve), and the large-w series (3.20) through order (dot-dashed pink curves). It is 
interesting to note that the small-w series agrees very well with the free-energy data in the 
regime for w > wo{sq) — —1/4 (in fact, the two curves he one over the other, at least out to 
w = 2); but it says nothing whatsoever about the behavior for w < Wq. On the other hand, 
the large-w series agrees well with the free energy for \w\ > 1, for both signs of w; indeed, it 
gives a fairly good description of the whole regime — oo < w < Wq, except very close to the 
transition point. 

The real part of the free energy clearly has a minimum at w — —1/4. From the plot it 
appears that the free energy is continuous at w = —1/4, while its derivative with respect 
to w (the "internal energy") has a jump discontinuity there. To verify this latter point, we 
computed the derivative of the free energy with respect to w at the same list of values of 
w, for widths L < 10, by numerical differentiation.^^ The results are shown in Figure 13. 
In this figure we have also shown the [20, 20] Fade approximant to the small-w series for 
the internal energy f'{w) [dashed violet curve] and the first derivative of the large-w series 
(3.20), namely f'{w) = 1/w — l/(8w^) [dot-dashed pink curves]. 

It is clear that f'{w) has a jump discontinuity at w = Wq = —1/4: when w t Wq it 
decreases towards a negative limit, and when w I wq it increases towards a positive limit. 
What is less clear from Figure 13 is whether these limiting values are finite or infinite. In 
order to address this question, we have plotted in Figure 14(a) a blow-up of the region very 
near w = Wq = —1/4, showing f'i^{w) together with the same small-'w; and large-w curves as 
in Figure 13. The curves f^iw) show a clear parity effect due to the existence for even L of 
the interval [wo-,tfo+] where the two dominant eigenvalues form a complex-conjugate pair, 
and due to the existence for odd L of the point wqq where the two dominant eigenvalues 
cross in modulus. More precisely, for even L the free energy has sqiiarc-root branch points at 
w — Wo±{L), so that f'L{w) diverges to ±oo at those two points, at a rate ~ \w — 'Wo±{L)\~^/^. 

^^For L < 9 we have done this computation using Mathematica with high-precision arithmetic (200 
digits). For L = 10 we have used the ARPACK package [146] with double precision. In the former case we have 
used the three-point interpolation formulas for central differentiation [152] with mesh width Aw = 0.5 x 10^^° 
(resp. Aw = 10~^°) for the first (resp. second) derivative of the free energy. In the latter case, we have used 
seven-point interpolation formulas with mesh widths Aw = 0.5 x 10~^ and Aw = 10~^, respectively. This 
is enough to get results with at least nine-digit precision. 
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These divergences are seen clearly in Figure 14(a) (note that wo+{L) = —1/4 for all L, while 
Wo^{L) for L = 6, 8, 10 is marked with a brown dot- dot- dashed vertical line). For odd L, by 
contrast, the free energy switches at w = Wqq{L) between two distinct eigenvalues, so that 
fiiw) has a jump discontinuity there. This discontinuity is also seen clearly in Figure 14(a). 

We now wish to infer from Figure 14(a) the infinite- volume behavior of of the internal 
energy f'{w) in the two regimes w > wo and w < wq. To do this, we must disregard 
the behavior too near the points w = wq±{L), where f'^iw) is dominated by the finite-L 
divergence, and focus instead on extracting the limit L — > oo at fixed w. Fortunately, a 
fairly clear picture emerges: 

1) For w > Wq, f'i{w) is decreasing (resp. increasing) in L for even (resp. odd) L; the two 
data sets seem to be converging to a common limiting curve that lies very near the [20,20] 
Fade approximant to the small-w series. Furthermore, this limiting curve very likely remains 
finite as w I Wq. Indeed, for odd L there is not the slightest indication that f'l^iw) could 
be diverging; the data seem, rather, to indicate a limiting value around 4, and certainly no 
more than 5. For even L the data are less clear, but at least for w ^ —0.248 we see that 
fiiw) is rapidly decreasing with L and apparently approaching a limiting value close to the 
Fade- approximant curve; furthermore, this convergence is rapid when w is well above wq 
and less rapid when w is nearer to Wq. It is therefore reasonable to expect that a similar 
convergence occurs for all w > Wq, but which for w very near Wq requires very large L to be 
seen definitively. We would guess a limiting value f^{wo) = lim^j^^^ /'(^) the range 4-5. 
However, it is at least conceivable that the infinite-volume internal energy f'{w) is weakly 
divergent ets w I wq, e.g. like f'{w) ~ \og{w — wq) or f'{w) ~ loglog(w — wq). 

2) For w < Wq the behavior is even clearer: both the odd-L and even-L curves stay 
close to each other (at a value —4.7) until they peel away due to the discontinuity at 
w — Wqq{L) for odd L or due to the proximity of the divergence aX w — Wq_{L) for even L; 
furthermore, these "peel away" points get closer to as L grows. This strongly suggests 
that a similar convergence occurs for all w < wq- Finally, the apparent limiting value ^ —4.7 
is not far from the curve corresponding to the large-w expansion (namely, about 20% above 
it). Once again, there is no evidence that f'{w) could be diverging as w f wq. Rather, the 
data for both odd and even L suggest a limiting value f'_{wo) = lim^j^Q /'(^) ~ —4.7. 

Since wo-{L) varies with L, it is useful to plot fiiw) versus w — wq-{L) in order to 
compare the behavior at "comparable" values of w. Such a plot is shown in Figure 14(b). 
We see clearly that the region of divergence grows narrower as L grows. In particular, at 
any fixed value of -u; — Wq^{L) < 0, f'^iw) appears to tend to a limiting value around —4.7: 
this is seen clearly for w — Wq_{L) < —0.003, and would presumably be seen also for smaller 
\w — wo-{L)\ if we were to go to larger L. This plot confirms, from a shghtly different point 
of view, our conclusion that /i(wo) = lini,/,|,/,o f'{w) ~ —4.7. 

In summary, there is strong evidence that the infinite- volume internal energy f'{w) tends 
to a finite value as u; | Wq. The evidence is less clear concerning the behavior as w | Wq, 
but the most likely scenario is also a finite limit, with weak divergences also possible. In any 
case, the existence of a discontinuity — whether finite or infinite — in the first derivative of 
the free energy is a clear signal of a first-order phase transition at wq- 

We can also compute the second derivative of the free energy with respect to w (i.e., 
the analogue of the specific heat in the usual Potts model), again by numerical differentia- 
tion. The results are shown in Figure 15, along with the [20, 20] Fade approximant for the 
corresponding small-u; series for f"{w) [dashed violet curve] and the second derivative of 
the large-w series (3.20), namely f"{w) — —l/w"^ -\- l/{Aw^) [dot-dashed pink curves]. This 
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derivative is clearly getting large and negative as we approach wq from both sides, and it 
may well be diverging to — oo. To check this point more carefully, we made a blow-up plot 
of f"{w) very near w = wq, using a logarithmic vertical scale [see Figure 16(a)]. Just as for 
/', we find a different behavior for odd and even L. The curves f'liw) for odd L have a jump 
discontinuity &X, w — wqq{L). Their value &X, w — wq — —1/ A grows with L and it is not 
clear how these values behave in the limit L — > oo. For fixed even L, we find that f'l{w) 
grows rapidly as w t '"^o-(-^) ^-nd as w | Wqj^{L) — —1/4; this rapid growth is once again 
explainable as arising from the square- root branch points in fiiw), causing a divergence 
f'li'^) ^ \w — Wq±{L)\~^/'^ . However, the magnitude of f'liw) at fixed w is decreasing as L 
grows, and it is unclear how f"{w) = Uhil^^ f'l{w) behaves. Once again, let us consider 
this question separately on the two sides of Wq: 

1) For w > Wq, it is plausible that the finite-L curves are converging from opposite sides 
(i.e., from above when L is even and from below when L is odd) to a limiting curve f'{w) 
that is close to the [20,20] Fade approximant to the small-w series (i.e., the dashed violet 
curve in Figure 15). Analysis of this Fade approximant (see below) then suggests that f"{w) 
diverges sts w I wq, with an exponent near a ~ 1. 

2) For w < Wq, let us plot f'liw) versus w — wo-{L) in order to compare the behavior at 
"comparable" values of w. Such a plot is shown in Figure 16(b). Careful inspection shows 
that the ratios of f"{L) for different L — i.e., the vertical distances between curves on this 
logarithmic plot — are decreasing as w j wq. (The strange behavior for L = 4 very near 
wo-{L) arises from the fact that this point is not an endpoint, but rather a T point, as 
shown in Figure 4.) This suggests that the infinite-volume specific heat f"{w) may diverge 
as w I in a manner somewhat similar to that shown by the L = 10 curve. But the 
evidence is admittedly weak, and it is fair to say that we cannot really be sure how f"{w) 
behaves as w t ^o- 

Finally, further information about the behavior of f{w), f'{w) and f"{w) as w I wq can 
be obtained by analyzing the small-w series for these quantities. This is done in detail in 
Appendix A; here we simply summarize the results. We find that the behavior of these three 
quantities is consistent with 

f{w) ~ + with q; = 0.94 ±0.10 (7.28a) 

f'{w) ~ {w + lY^'' with a = 0.90 ± 0.04 (7.28b) 
f"{w) ~ {w + iy with a = 0.91 ± 0.02 (7.28c) 

We thus have strong evidence that f"{w) diverges aswl wq, but it is not clear whether f'{w) 
diverges there or not. Taking our estimates literally, the corresponding critical exponent is 
1 — a = 0.10±0.04 > 0, so that /' is finite aswlwo; but we cannot make a firm conclusion as 
the error bar is very large! Taking into account the uncertainties of series analysis, it seems 
fair to say that (7.28c) is consistent with the theoretically predicted value a — 1, especially 
if we recall the possibility of multiplicative logarithmic corrections. Note that, depending on 
the form of the latter corrections, /" could be either integrable or nonintegrable at w = Wq, 
i.e., /' could be either finite or infinite there. 
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7.5 Square lattice: Estimates of the free energy and central charge 

Baxter's exact solution [28] for the Potts- model free energy on the curve (1.3) implies 
the exact value 

/(sq,^ = -l/4) = 21og(^^|^^ ^ -0.6031055757. (7.29) 

In order to make a high-precision test of this prediction, we computed the strip free energies 
fziw — —1/4) for all widths L up to 15. For 2 < L < 9 we used our symbohc transfer 
matrices, while for 10 < L < 15 we used a numerical transfer-matrix algorithm. These data 
are reported in Tabic 11. The free energies fiiw = —1/4) show strong even-odd oscillations; 
moreover, they appear to be monotonically increasing in L for even L, and monotonically 
decreasing for odd L. If this monotonicity persists for larger L, then from L = 14, 15 we can 
deduce the bounds 

-0.6067349394 < /(sq,w = -l/4) < -0.6019335900. (7.30) 

Separate extrapolations of the even and odd subsequences using the Ansatz (7.27) yield the 
results reported in Table 12. Putting them together, we conclude that 

/(sq,'fi; = -l/4) = -0.60310 ± 0.00010 , (7.31) 

in excellent agreement with the exact result (7.29). We also find a rate of convergence 
compatible with A = 2 (see below for a more detailed analysis). 

We have similarly extended up to width L = 16 the data for the points w — ±1 oi 
interest to mathematicians: see Table 11. In both cases. Re f{w) appears to be monotonically 
increasing in L. If this monotonicity really holds for all L, it would imply that 

Re/(sq,'u; = 1) > Re/i=i6(sq, -u; = 1) 1.3079471010 (7.32a) 
Re/(sq,'fi; = -l) > Re /L=i6(sq, w = -1) fa 1.0432380373 (7.32b) 

The lower bound (7.32a) hes between the rigorous lower and upper bounds (2.41). Extrap- 
olation of these sequences using the Ansatz (7.27) for triplets L — L^in, L^m + 1, -^min + 2 
yields the results reported in Table 12. The convergence for w — —1 is slower (A 2) than 
ior w — 1 (A fti 3.6). Our best estimates /(sq, w — ±1) are 

/(sq,'u; = l) = 1.30819 ±0.00010 (7.33a) 
/(sq,'u; = -l) = 1.04870 ±0.00010 (7.33b) 

These values are of course larger than the conjectured lower bounds (7.32); and the estimate 
for /(sq, w = 1) lies between the rigorous lower and upper bounds (2.41). 

The different values found for the correction exponent A can be explained in the following 
way. For w > —1/4, we expect that the system is governed by the high-temperature fixed 
point, which is non-critical. As a consequence, we expect an exponentially rapid convergence 
of the strip free energy to its infinite- volume limit. This explains why for w = 1 the estimates 
of A grow with I/mm, apparently without bound. For w < —1/4, by contrast, we expect 
that the system will be governed by the fixed point located inside the Berker-Kadanoff 
phase [46, 47] , which is critical. Conformal field theory (CFT) then predicts that 

h{w) = f(w) + ^c{w)L-' + ... (7.34) 
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where c(u') is the central charge (see Section 2.3); we therefore expect A = 2. Moreover, 
the central charge associated to the critical point of the g = Potts model in the Berker- 
Kadanoff phase is predicted from (1.9) to be c = —2. Finally, at w = —1/4 there is a phase 
transition that we expect also to be a critical point; we therefore again expect A = 2. Here, 
however, the central charge is predicted from (1.6) to be c = —1. 

Let us therefore use our finite-size data to estimate the central charge c{w). As a first 
approach, we have fit the data to the Ansatz (7.34) with no higher-order corrections. The 
results are displayed in Figure 17(a). We observe that there are two clear regimes separated 
by w = Wo = —1/4, in qualitative agreement with the above discussion. 

For w > Wo, the estimates for the central charge are very close to c = when wq < w < 
0.4. This agrees well with the prediction that the system is non-critical (the convergence 
of fhiw) to the bulk free energy f{w) is exponentially fast, so that c = 0). For w ^ 0.4 
we observe sizable downward deviations from c = 0, even for large values of Lmji,; but the 
estimates of c increase monotonically and are at least consistent with convergence to the 
value c = 0. These deviations from c = can be explained as due to crossover from the 
marginally repulsive (i.e., asymptotically free [24]) ferromagnetic fixed point aX, w — -|-oo, 
which by (1.4) is a c = —2 theory. In this regime, we expect large finite-size corrections that 
are not adequately handled by our Ansatz (7.34). 

For w < Wo, we find that the estimates for the central charge are quite far away from the 
predicted value c— —2 (at least up to Lmin = 13). However, these estimates are increasing 
monotonically with Lmin, and it is at least plausible that they are approaching c = — 2 as 
-^min oo. Furthermore, the estimates are closer to c = —2 for larger negative values of w 
(i.e., for points deeper inside the Berker-Kadanoff phase). Indeed, for large negative values 
of w we are closer to the renormalization-group fixed point (1.2_), which lies at w = — oo; 
therefore, the corrections to scaling are expected to be weaker for larger negative values 
of w. All these findings refiect the marginally attractive nature of the RG fixed point at 
w = — oo [24]. Indeed, the theory predicts very slow convergence to the infinite- volume 
limit, as a consequence of the log log L/ log L and 1/logL corrections to scaling arising 
from the marginally irrelevant operator. It is quite plausible that the behavior exhibited in 
Figure 17(a) is a manifestation of such logarithmic corrections to scaling. 

It is instructive to plot the estimates of c{w) versus 1/w, as shown in Figure 17(b). The 
smooth behavior of the plot near 1/w = suggests that the fixed points lying at w = ±oo are 
in fact identical and have central charge c = —2, in agreement with the prediction [24] that 
this is a theory of a pair of free scalar fermions. The L-dependence of the estimates of c{w) 
— tending towards — 2 as L — >• oo when 1/w < 0, and towards as L — >• oo when 1/w > — 
is likewise in agreement with the prediction [24] that this fixed point is marginally repulsive 
ior g = 1/w > and marginally attractive for g < 0. 

Somewhat stabler estimates for the central charge can be obtained if we fit the numerical 
data to the Ansatz 

hiw) = f{w) + lc{w)L-^ + AL-' . (7.35) 

As explained in Section 2.3, the term is predicted by conformal-field-theory arguments; 
but we cannot rule out the existence of additional correction terms, lying between L^^ and 
L"*^, arising from irrelevant operators with critical exponents between y = —2 and y = —4. 
In practice, inclusion of the term usually tends to accelerate the convergence of the 
estimates for the central charge [153,154]. The resulting plot of c{w) versus w is shown 
in Figure 17(c). We observe that the estimates for c{w) converge more rapidly than with 
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Ansatz (7.34), although we are still far away from the predicted value c = — 2 in the regime 
w < —1/4. For example, at w = —2 we obtain c ~ —2.487 at L^ain = 12 with the Ansatz 
(7.35), compared to c ~ —2.543 at Lmin = 13 with the Ansatz (7.34). The corresponding 
plot of c{w) versus 1/w is shown in Figure 17(d). 

The mapping of the spanning- forest model to the A?^- vector model at N — —1 [24] suggests 
that a more correct Ansatz for the free energy in the critical phase would be 

hM - fM + ^c(»)L- + AL-^'^ + . (7.36) 

The plot of c{w) obtained from a fit to (7.36), restricted to the regime w < —1/4, is shown 
in Figure 18.'^'^ We obtain values for the central charge closer to the expected result c = —2 
than for the simpler Ansatze (7.34)/(7.35). For instance, for w = — 2 we get c = —2.161 at 

^min — 11- 

For w = —1 we have extended the above analysis to include finite-size data up to L = 16 
(see Table 13). The estimates for c{w = —1) continue to increase monotonically, but we are 
still far away from the predicted value c = —2. We obtain values slightly closer to c = — 2 
using the Ansatz (7.35) than with the original Ansatz (7.34), but we are still far away. As 
a side benefit, we obtain an estimate for f{sq,w — —1) that is slightly more stable than 
(7.33b): 

/(sq, w = -1) = 1.04883 ± 0.00008 . (7.37) 

If we use the improved Ansatz (7.36), the estimates for c{w — —1) still increase monotonically 
for Ljniii > 5 (see Table 14), but these estimates are now much closer to the predicted value 
c = —2 (e.g., for Lmin = 13, we get c ~ —2.046). Finally, this Ansatz provides a more stable 
and precise estimate for /(sq, w = — 1): 

/(sq, = 1.048755 ± 0.000008 . (7.38) 

For w — —1/4 we already know the exact value of the bulk free energy f{w), namely 
(7.29). We therefore used this value and performed a one-parameter fit to (7.34) to estimate 
c{w = —1/4): see Table 15. For even L, the estimates of c increase monotonically from 
c ~ —1.79 at Lmiii = 2 to c ^ —1.36 at Lmin = 14. This suggests that the central charge for 
even widths would eventually converge to the predicted value c = —1 in the limit Lmin oo. 
For odd L, by contrast, the estimates of c decrease monotonically from c 0.556 at Lmin = 3 
to c 0.5036 at Lmin = 15. This latter value strongly suggests an effective central charge 

Ceffective(sq,W = -1/4) = 1/2 . (7.39) 

Similar results are obtained by the two-parameter fit to the Ansatz (7.35) with the bulk free 
energy fixed to its exact value (7.29). For even L, we obtain values of c{w = —1/4) slightly 
closer to the predicted value c = — 1 (the estimates increase monotonically from c ~ —1.57 at 
-^min = 2 to c ~ —1.31 at Lmin = 12). For odd L the estimates are stable around c ~ 0.501. 

The free energies for even L are related to the ground state of the corresponding con- 
tinuum theory, since frustration effects are absent. Similarly, it is natural to assume that 
the free energies for odd L are related to the ground state of a modified continuum theory 

^^For w = —0.3 we found strong even-odd oscillations. Therefore, for this particular value we performed 

the fits using the data with L = Lmin, Lmin + 2, Lmin + 4, Lmin + 6. 
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in which the fields are subjected to a twist in the spatial direction (antiperiodic boundary 
conditions). If this is indeed the case, conformal field theory [30-32] predicts that for L odd, 

Ceflfcctive = C - 24/itwist • (7.40) 

Assuming the conjecture (7.39) to be correct, the conformal weight of the twist operator 
must then be 

^twist = -1/16 . (7.41) 

The twist operator was not discussed in the original paper of Saleur [47] on the CFT analysis 
of antiferromagnetic Potts models. However, recent work [155] on the antiferromagnetic 
critical line (1.3) has identified the twist operator with the fundamental disorder operator 
[156] of parafermions, where 5 is defined by (1.1). In the limit g ^ this confirms the 
conjecture (7.41). 

7.6 Square lattice: Estimates of the thermal scaling dimensions 

The first step in obtaining the thermal exponent v is to obtain the inverse correlation 
lengths for a square-lattice strip of width L [cf. (2.36)]: 

(7.42) 

where A*(w) is the largest eigenvalue (in modulus) of the transfer matrix associated to a 
square-lattice strip of width L and cylindrical boundary conditions, and \i{w) is the i-th 
subdominant eigenvalue. In Figures 19(a) and (b) we have plotted the first two inverse 
correlation lengths ^~^{w] L) for all widths L < 10 in the range —2 <w<2. 

The second step is to extract from this finite-width data the thermal scaling dimensions 
XT,i- If w is a noncritical point, we expect as usual an exponentially rapid convergence: 

= C'W + 0(e-^^^) (7.43) 

where Q < Ai{w) < oo. If w is a critical point, then the behavior of the inverse correlation 
length can be described by the CFT prediction (2.37). In order to fit the data in an unified 
fashion, we began by using the Ansatz 

C'Ki)=C'W + ^^. (7.44) 

For a critical system, we expect .^^^^ = 0, and then XT,i is the true scaling dimension for that 
system, modulo higher-order finite-size-scahng corrections that are neglected in (7.44). For 
a noncritical system, we expect that > 0. In this case, the term in (7.44) tries to 
mimic the exponentially small corrections, so that the estimates of XT,i should tend rapidly 

to zero as L grows. 

We began by fitting our finite- width data to the Ansatz (7.44). We found, as expected, 
three different regimes, according as w < —1/4, w — —1/4 oy w > —1/4. 

Let us start by discussing the first correlation length ^i: see Figure 19(a) for the finite- 
width data as well as the extrapolated infinite- width values (black solid circles). Throughout 
the regime w < —1/4, the estimated value of ^^"^ is small in magnitude, in agreement with 
the prediction that this phase is critical. Indeed, for w ^ —0.5 we find that ii^{w) ^ 0.095 



A^(w) 
\i{w) 
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and the estimates decrease steadily with Lmin- For —0.5 ^ w ^ —1/4, the estimates of 
are larger, but they again decrease as Lmin grows. We interpret these deviations of ^ 
from zero as arising from corrections to scaling, which not surprisingly become stronger as 
the transition point wq = —1/4 is approached. We take all this as an indication that the 
true value is ii^{w) = throughout the region w < —1/4. We therefore repeated the fit, 
imposing ^ = 0. The resulting estimates for the scaling dimension xt,i are displayed in 
Figure 20(a). We see that the scaling dimension is close to the value xt,i{w) — 2 predicted by 
Saleur [47] for the Bcrkcr-Kadanoff phase w < —1/4 [cf. (1.10)]. Furthermore, as we proceed 
deeper into the Berker-KadanofF phase (i.e., towards w more negative), the estimates also 
become closer to xt,i = 2. 

For w > —1/4, the estimated value of is strictly positive (corresponding to a non- 
critical theory), at least up to w ~ 0.8. Around w — the inverse correlation length becomes 
very large, as expected; as w increases beyond this point, ^ decreases. For w > 0.8 our 
estimates of ^ become small (< 0.18); but we do not take this fact as an indication of a 
critical system. Rather, it simply reveals the fact that the correlation length has become 
comparable to the strip widths we are considering, and that the Ansatz (7.44) is not adequate 
to describe accurately the corrections to scaling in this regime (indeed, there is no evidence 
of a phase transition around w ~ 0.8). We expect, in fact, that the correlation length is 
finite for all w > 0, tending to infinity as w — > +00. Our estimates for ^j"^ are depicted in 
Figure 19(a) as black sohd dots. 

The case two = —1/4 is rather special. We have already seen that for even widths, 
this point is an exact endpoint of the limiting curve. Therefore, ^^^{w = -1/4; L) = 
exactly for all even L. In particular, from even L we conclude that x^^i = 0. The values of 
Ci^{w = — 1/4; L) for odd L are displayed in Table 16, where we also show (in columns 3-5) 
the results of fitting the data with odd L to the Ansatz (7.44). We see that the estimate 
for ^ is a small negative number, which decreases in modulus as Lmin is increased. This 
suggests that .^"^ = exactly, so we repeated the fit with this value fixed. The results 
are shown in the last two columns of Table 16. As we increase L^^^, the estimate for Xt,i 
decreases; but it is not clear whether it is tending to zero or to a nonzero value around 
xt,i ~ 0.3. What is clear is that Saleur's conjecture (1.7) for the thermal exponent along 
the critical antiferromagnetic curve, which reduces to xt,i = 2 at = 0, is incorrect, as our 
estimates are clearly smaller than this value. If it is indeed the case that xt,i = (as we 
suspect), then the correlation- length exponent is 

1/ = = - , (7.45) 

2 - XT,i 2 ' ^ ^ 

in agreement with our findings from the finite-size shift of the critical point (Section 7.3 and 
Table 8) and with the existence of a first-order phase transition at w — wq- 

The second inverse correlation length is displayed in Figure 19(b). If we try to 
estimate the thermodynamic limit of ^2 using the Ansatz (7.44), we find that 0.56 < < 
0.74 in the regime w < Wq. At first sight this could imply that this correlation length is finite 
when L 00. However, because the correlation length ^1 actually diverges in this regime, 
we expect that ^2 will do so as well. On the other hand, we expect that as the index i grows, 
the finite-size corrections to $,i{L) will become larger. We therefore repeated the fit, fixing 
= for w < Wq and attempting to extract Xt,2 from the coefficient of the linear term in 
1/L. The results are displayed in Figure 20(b). We see that in the Berker-Kadanoff phase 
the value of this scaling dimension is close to xt,2 = 4. 
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For w > Wo, we find a pattern similar to what was observed for ^i: namely, the inverse 
correlation length grows as we approach w = 0, and then it decreases, getting stable aronnd 
~ 2.5. As we concluded that remains finite in this regime, ,^2 should stay finite too; 
but the finite-size corrections might be larger than for ^1 (as happens in the w < Wq regime). 

The values of a,t w = wq = —1/4 are displayed in Table 17. We observe two different 
behaviors depending on the parity of the width L. For even L we find that the estimates 
for Xt^2 decrease as Lmin increases. Furthermore those estimates are not far away from those 
obtained for Xt,i with odd L (see Table 16). For a fixed value of L„^i^, the estimates for Xt^2 
(even L) are systematically larger than the corresponding ones for xt,i (odd L); but it is 
unclear whether the hmiting values coincide or not.^^ On the other hand, the estimates for 
xt,2 for odd L grows as Lmin increases, and we can conclude that xt,2 ~ 1.25. 

As the finite-size effects in the quantities ^^^(L) seem to grow with the index i, and we 
have seen that these corrections are important ah^cady for i = 2 (at least for the widths 
considered here), we did not try to analyze the correlation lengths i > 3. However, it is 
worth noting the possibility that some of the operators corresponding to i > 3 are (like 
the ones corresponding to i = 1, 2) more relevant than the operator with xt — 2 that was 
conjectured (incorrectly, we now know) by Saleur [46,47] to be the leading thermal operator 
[cf. (1.7)]. 



7.7 Triangular lattice: Estimates of critical points and phase bound- 
ary 

The triangular-lattice data that we wish to extrapolate arc collected in Table 6 and 
depicted graphically in Figure 11. The quantities Wo-{L) and Wo+{L) are defined only for 
even L, so we must perform fits to the Ansatz (7.2) using L = Lmm, L^in + 2, Lmin + 4. 
The quantity wqq{L) is defined for all L; and for the triangular lattice (contrary to what 
happened for the square lattice) there do not seem to be significant even-odd oscillations. 
We can therefore perform fits to (7.2) using L = L„^in, L„^in + 1) -^min + 2. 

The results of these fits are displayed in Table 18. The estimates of wq are very well 
behaved (particularly those based on wqq), and from Lmin = 6, 7, 8 we can estimate 

wo{tri) = -0.1753 ± 0.0002 . (7.46) 

This value is far away from — 1/6 ~ —0.1667, and so rules out the naive conjecture that 
wo{C) — —1/r whenever £ is a regular two-dimensional lattice of coordination number r. 
As in the square-lattice case, the exponent A seems to approach A = l/i/ = 2as Lmin 
increases. This fact agrees with the interpretation that there is a first-order phase transition 
ai w = Wq. In Table 19 we present the results of fitting the data to the modified Ansatz 
(7.11): the estimates converge slightly faster, and we again find 

woitri) = -0.1753 ± 0.0002 . (7.47) 

If we fit the width of the interval [wq-, ^0+] for even L to the Ansatz AL~^, we find 

ouui ~ 6.0 , (7.48) 

^^Please notice that if we merge the ^ data for odd L (displayed in Table 16) with the data for even 
L (displayed in Table 17), we obtain a data set with notable even-odd effects that prevent us from doing any 
joint fit. 
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which is three times larger than Usq = A ?a 2. In other words, on the triangular lattice the 
interval [■u;o-,u'o+] is extremely narrow, and gets rapidly narrower as L grows. We do not 
know the reason for this different behavior on the two lattices. 

Next let us extrapolate the endpoints wb{L), handling separately their real and imaginary 
parts. Once again, the real part of wb{L) seems to diverge as L — > oo like ci logL + C2- We 
have used the Ansatz (7.12), and the estimates for ci and C2 as a function of L^in are shown 
in columns 2-3 of Table 20(a). As for the square lattice, the estimates of Ci are rising with 
Lmin and not yet slowing down. This suggests that the true value of Ci is quite a bit larger 
than the value ~ 0.236 observed at our largest Lmm- If forced to guess, we might estimate 

ci = 0.24 ±0.02 (7.49a) 
C2 = -0.27 ±0.04 (7.49b) 

The relation with the = — 1 vector model [24] predicts that 

RewB{L) = cilogL + cgloglogL + C2 + o( i^s]^ ) (7.50) 



logL 

with ci = V^/(27r) fa 0.275664 and C3 = -1/(2^377) -0.091888 [24]. Indeed, the ratio 

44 = 44 - 1.732051 (7.51) 

ci(tri) C3(tri) ^ ^ 

is characteristic of asymptotically free theories, and simply comes from comparing the lattice 
free propagators with their continuum limits (see e.g. [157, eq. (5.7)]). 

In the last column of Table 20(a), we display the estimates of the parameter C2 obtained by 
performing the one-parameter fit to the Ansatz (7.50) with ci and C3 fixed to their theoretical 
values (and the correction terms 0(log logL/ logL) neglected). Our estimate is 

C2 = -0.285 ±0.005 (7.52) 

However, the very-slowly-varying nature of the expected correction terms and the narrow 
range of L values available here make plausible that the true value for C2 differs from the 
above estimate by many times our estimated error. 

The imaginary part of u'b(L), by contrast, seems to converge to a finite value as L — > cxd: 
the estimates using the Ansatz (7.4) are displayed in Table 20(b). We find 

Im wb = 0.394 ±0.004 (7.53) 

with A = 0.63 ±0.05. 

Finally, we tried to extrapolate the whole limiting curve B to the infinite- volume limit, by 
the same methods used for the square lattice (Section 7.3). We began by extrapolating Hew 
at the selected values Imw = 0.01,0.02, . . ., using the Ansatz (7.4) for Rew. We obtained 
consistent estimates of Re w for | Imw| < 0.23 (see Table 7). The results are shown by black 
dots in Figure 11. Next we attempted to complete this curve at larger Rew, using once again 
the theoretical predictions from the perturbative renormalization group at small g = 1/w. 
Unfortunately, for the triangular lattice the three-loop and four-loop computations have not 
yet been performed; therefore, all we have available are the universal values for the first two 
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beta-function coefficients. Using the relation /^continuum = V^Ptn [157, eq. (5.7)] to match 
normahzations, we obtain 



62 = -wo/Vs = V3/(27r) (7.54a) 



N=-l 



h = ^'i/SU^i = -l/(27r)' (7.54b) 

and hence 

Ai = h/h = -l/(2\/37r) . (7.55) 

Using the Ansatz (7.25) with /3o = 0.394 and -B3 = -B4 = and imposing the first derivative 
at a = 00 along with the estimated value and derivative at a = 0.0198, we obtain 

ao = -0.550842 (7.56a) 
Bi = -1/(27371) ^ -0.091888 (7.56b) 
B2 = -0.122843 (7.56c) 

The corresponding curve (7.25) is depicted as a black dotted-dashed curve in Figure 11. 

7.8 Triangular lattice: Behavior of the free energy and its deriva- 
tives as w ^ wq 

We have computed the real part of the strip free energy Re/L(u') for widths 2 < L < 13 
and w — —2, —1.9, . . . , 2, as well as w — —0.175. These results are displayed in Figure 21. 
The solid black line represents the extrapolation of the free energy to L = (X) using the 
Ansatz (7.27). We also show the [10,10] Pade approximant to the small--u; series for this 
quantity (dashed violet line), and the large-?!; series (3.20) through order w^^ (dot-dashed 
pink curves). Once again, the small-w series expansion agrees very well with the free-energy 
data in the regime w > Wq. On the other hand, the large-w expansion gives very accurately 
the free energy in the whole regime w < Wq. The free energy has clearly a minimum at a 
value close to the estimate Wo{tYi) ~ —0.1753. 

The first derivative of the free energy clearly has a jump discontinuity at wq, as is seen 
on Figure 22. Again, we have made a blow-up picture (Figure 23) in an effort to determine 
whether this discontinuity is finite or infinite. As for the square-lattice case we find a clear 
parity effect: for odd L, f'l{w) has a jump discontinuity w — Wqq, while for even L we 
find divergences f'l{w) ~ — Wq^I'^/"^ due to the square-root branch points w = Wq±{L). 
Note that for the triangular lattice the intervals [ifo- j ""^0+] ioi even L are extremely narrow 
[cf. (7.48)]; this makes it much easier, compared to the square lattice, to extract the limiting 
behavior as L ^ 00 with w fixed, as it suffices to exclude from consideration a very narrow 
region surrounding the interval [t^o- , tWo+] ■ In fact. Figure 23 strongly suggests that the 
discontinuity in f'{w) is finite at -u; = Wq. For w ], Wq, we observe that /' is very well 
approximated by the [10, 10] Pade approximant to the small-w series for /'. Thus, we can 
conclude that /^(wq) = lim^j^^^ /'(w) ~ 3.389 (this value comes from evaluating the Pade 
approximant ai w — wq ^ —0.1753). On the other side, the curves f^iw) for w < wq- 
(even L) or w < Wqq (odd L) stay rather close to (and slightly above) the first derivative 
of the large-w series (3.20), i.e. f'{w) — 1/w — l/{12w'^). We conclude that fL{wo) = 
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limyjiyjg f (w) ~ —8.14.^^ Thus, the jump in the internal energy is apparently finite and 
takes the value A/'(wo) ~ 11.53. The existence of a discontinuity in the first derivative of 
the free energy implies that there is a first-order phase transition at w — Wq also for the 
triangular-lattice model. 

Since wo-{L) and wq-^-{L) vary with L, it is useful to plot fiiw) versus w — wq-{L) or 
w — wo-^.{L) in order to compare the behavior at "comparable" values of w. Such plots are 
shown in Figure 24(a,b). For both w < Wq and w > wq, we see that the region of divergence 
grows rapidly narrower as L grows, so that at any fixed value oi w — Wq^{L) < [resp. 
w — Wo+{L) > 0], fi{w) clearly tends to a limiting value around —8.0 (resp. 3.4). This plot 
confirms, in a most striking way, our conclusions in the preceding paragraph. The evidence 
for a finite limit of f'^iw) as w — > from either side is, in fact, vastly more compelling for 
the triangular lattice than for the square lattice. 

Wc have also computed the second derivative of the free energy: see Figure 25 for an 
overview plot, and Figure 26 for a blow-up plot near wq using a logarithmic vertical scale. 
Once again, for odd L, f'l{w) has a jump discontinuity sX w — woq{L), while for even 
L we find divergences f'l{w) ~ — wo±{L)\''^/'^ due to the square-root branch points at 
w — Wo±{L). However, the divergences of f'l{w) for w ~ Wo± are extremely narrow, which 
allows us to extract with fair confidence the limit L — > oo at w fixed Wq, as follows: 

1) For w slightly above Wq, the specific heat f'l{w) stays very small (and almost in- 
dependent of L) until w goes fairly close to wo+{L). Moreover, the width of this "diver- 
gence" region becomes smaller as L grows. This suggests that the infinite-volume specific 
heat f"{w) remains bounded (and in fact) small as w | wq, in sharp contrast to what 
was observed for the square lattice. Indeed, the behavior for w > Wq is well described for 
L > 5 by the [10, 10] Fade approximant to the small-w series, from which we can con- 
clude that /+(wo) = lim^j^^Q f"{uj) ^ —1.187 (based on evaluating the Fade approximant at 
w^wo^ -0.1753). 

2) For w slightly below Wo, by contrast, we find large negative values for f"{w), although 
these large values do seem to remain finite as w ^ Wq. Indeed, the data points for w < wq_ 
(even L > 6) or w < Wqq (odd L > 7) remain very close to the second derivative of the 
large-w expansion (3.20), namely f"{w) = — l/w^ -|- l/{6w^). This suggests that f"{w) does 
not tend to — oo as w t wq, but rather to the finite value f'l{wo) = lim^-f^g f"{w) ~ —63.5 
(based on evaluating the large-w series at = ~ —0.1753). 

In Figure 27(a,b) we plot f'l{w) versus w — Wo-{L) or w — «;o+(-^) in order to compare 
the behavior at "comparable" values of w. These plots confirm our conclusion that f'l{w) 
stays almost independent of L except in a narrow region around Wo±{L), which further- 
more becomes narrower as L grows, so that the infinite- volume specific heat f"{w) remains 
uniformly bounded as w — >■ Wq from cither side. 

These findings on the behavior of the free-energy derivatives merit some discussion. We 
have seen that, on the triangular lattice, both /' and /" apparently remain finite as w Wq 
from either side; at Wq these quantities simply jump from one finite value to another. This 
behavior is very different from what was found for the square lattice (Section 7.4), where /" 
very likely diverges at wq (at least when approached from above), and /' may conceivably 
diverge as well. At least for fixed w < wq (i.e., inside the Berker-Kadanoff phase) one would 

^^We have empirically found that for large values of L and for w | wq, the quantity 5f'{w) = f'[^{w) — 
[l/w — l/(12u>^)] behaves approximately linearly in w. In fact, a rough fit of our data for L = 9. 10 gives 
6f'{w) « 0.727+ 2.580'u;. The value fLiwo) ~ —8.14 comes from evaluating both the large-w series and the 
fit for 6f'{w) at w = wo ~ -0.1753. 
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expect some universal behavior (as indeed happens for the central charge and the thermal 
scaling dimensions xx^i)- So the discrepancy between the square and triangular lattices is 
puzzling; see Section 8.3 for further discussion. 

7.9 Triangular lattice: Estimates of the free energy and central 
charge 

The real part of the triangular- lattice free energy at w = ±1 is monotonic in L, at least 
as far as we have been able to compute it (see Table 21). If this property continues to hold 
for larger L, we can conclude that 

Re /(tri, w = 1) > Re /L=i5(tri, w ^ 1) ^ 1.7010224200 (7.57a) 
Re/(tri,'u; = -1) > Re/i=i3(tri, w = -1) ^ 1.5273746626 (7.57b) 

The lower bound (7.57a) falls between the rigorous lower and upper bounds (2.43). 

If we fit the finite-L data at w — ±1 to the Ansatz (7.8), we get the results displayed 
in Table 22. At w = —1, the estimates of the exponent A are very close to A = 2 and 
apparently tending to it, exactly as expected for the Berker-KadanofF phase. At -u; = 1, 
the estimated exponent A is slightly larger than 2, and slowly rising; this is compatible 
with the idea that w = 1 is a non-critical point which has, however, a reasonably large 
correlation length, so that the system is strongly affected by crossover from the critical point 
a,t w — -|-oo.^^ Our preferred estimates for Re /(tri, w) are: 

Re/(tri,«; = 1) = 1.70255 ± 0.00010 (7.58a) 
Re/(tri,w = -1) = 1.53413 ± 0.00007 (7.58b) 

These values are indeed larger than the conjectured lower bounds (7.57); and the estimate 
for Re /(tri, to = 1) is smaller than the rigorous upper bound (2.43). 

Finally, we can estimate the central charge of the model as a function of w, using the 
same method as for the square lattice. In this case, the Ansatz (7.34) has to be modified to 
include the geometrical factor (2.35): 

/(Lp X OOF, w) = f{w) + > L-^ + . . . . (7.59) 

Wc began by fitting the data to the Ansatz (7.59) with no higher-order correction terms. 
The results are displayed in Figure 28(a). There are two clear regimes separated by Wo(tri) a; 
—0.1753, and the values of the central charge agree qualitatively with those found for the 
square lattice. In Figure 28(b) we show the same plot versus l/w; the behavior near 1/^ = 
is once again in agreement with the prediction [24] that the fixed point at 1/w = has central 
charge c — —2 and is marginally repulsive for = 1/w > and marginally attractive for 
g<Q. 

In Figure 28(c) we show the estimates for the central charge c{w) coming from the 
modified Ansatz in which an L""^ term is included. The estimates seem to converge faster in 
the region w < »o. l)ut they are still far from the expected value of c = —2. In Figure 28(d) 

^^In Figure 30(a) we observe that the extrapolated value for the inverse correlation length at t« = 1 is very 
close to zero. In fact, our best estimate is a small negative number, ^(w = 1) rj —0.01. 
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we show the corresponding plot versus 1/w. Finally, we show in Figure 29 the estimates 
for c{w) for w < wq using the improved Ansatz (7.36). As for the square- lattice case, these 
results agree better with the expected result c = —2. 

For w = — 1 we have extended the above analysis to include finite-size data up to L = 13 
(see Tables 23 and 24). As for the square-lattice case, we find a good agreement with the 
predicted value c — —2 only for the improved Ansatz (7.36): for Lmin = 10 we find c{w — 
— 1) —2.077. Finally, we obtain a more stable and precise estimate for /(tri, w = —1): 

/(tri, w = -l) = 1.534166 ± 0.000008 . (7.60) 

7.10 Triangular lattice: Estimates of the thermal scaling dimen- 
sions 

As for the square lattice, we began by obtaining the inverse correlation lengths ^^^{w; L) 
for finite- width strips using the definition (7.42). The first and second inverse correlation 
lengths for widths L < 10 are plotted in Figures 30(a) and 30(b), respectively. 

The scaling exponent XT,i is obtained by fitting the finite- width data to the Ansatz 

ii {w] L) = [w) + , (7.61) 

which includes the geometrical factor (2.35). The applicability of this Ansatz is subject to 
the same comments made in Section 7.6. The estimates of the infinite- volume correlation 
length are depicted in Figures 30(a) and 30(b) by black solid dots. 

Let us begin by discussing the first correlation length ^i. As expected, we found clearly 
different behavior in the regimes w <Wq and w > Wq. 

For w < Wq, the estimated infinite-volume inverse correlation length is compatible with 
the value ^ = 0. Indeed, we find estimates that are small negative numbers, whose absolute 
value decreases as we increase I/min; we therefore expect that they will converge to zero as 
L — > oo. We therefore repeated the fit, fixing — 0. The results for xt,i are displayed 
in Figure 31(a). Inside the Berker-Kadanoff phase w < wq, we again find that the scaling 
dimension is close to Saleur's prediction xt,i = 2. The scaling dimension at the transition 
point w = Wo is difficult to obtain, as we do not know the exact value of wq and the 
correction-to-scaling effects are very large near this point. 

For w > Wq, by contrast, we find a clear non-zero estimate for ^, at least up to w ~ 0.5. 
For w > 0.6, our estimates for are rather small (,^f ^ < 0.07); but as before, we believe 
this indicates only that has become comparable to the strip widths and that we are unable 
to handle the strong corrections to scaling in this regime. (Indeed, we have not found any 
evidence of a phase transition 0.5.) 

For ^2 the overall picture is very similar, although the finite-size corrections are larger (as 
happened in the square-lattice case). For w < —0.5, the data are compatible with ,^2^^ = 0: 
indeed, our estimates for ,^2^"'^ are small negative numbers, whose absolute value decreases as 
we increase Lmin- In the interval —0.5 ^ w < Wq, we find strong correction-to-scaling effects; 
but again we expect that the correlation length will diverge, as it does for ^1. We therefore 
repeated the fit, fixing = 0. The values obtained for xt,2 are displayed in Figure 31(b). 
In this regime the estimates are close to Xt,2 = 4, as for the square lattice. In the regime 
w > Wq, the behavior is similar to that of ^1: ^2^{L) tends to a finite value, which goes to 
zero as w — >• 00. 
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8 Discussion 



8.1 Behavior of dominant-eigenvalue-crossing curves 3 

We have computed the symbohc transfer matrices for square- and triangular-lattice strips 
of widths 2 < L < 10 with cylindrical boundary conditions. From these matrices one can 
compute (a) the generating polynomial of spanning forests Fg{w) for a cylindrical strip of 
width L and arbitrary length n, along with its complex zeros, and (b) the accumulation 
points of such zeros in the limit n — > oo. According to the Beraha-Kahane-Weiss theorem 
(Section 2.4), these zeros accumulate at certain points (the isolated limiting points, which 
seem not to exist in the models studied here) and along certain curves (the limiting curves 
B). By studying the behavior of these limiting curves as a function of the strip width L — > oo, 
we hope to shed light on the thermodynamic limit of the model, in which L, n — > oo. 

For strips of widths L < 6, we were able to compute the limiting curves B using the 
resultant method [5]; this method allows a complete determination of B, including all its 
endpoints. For widths 7 < L < 10, we used the direct-search method, which gives only 
lower bounds on the numbers of endpoints, as it could miss some. For L = 7, 8 we were 
able to obtain a fairly good description of B, but for L > 9 we were only able to obtain the 
points at the discrete set | Imw] = 0,0.01,0.02, etc. These hmiting curves B are shown in 
Figures 7 and 11 for the square and triangular lattices, respectively, and their characteristics 
are summarized in Tables 5 and 6. In all cases we can verify the identity 

endpoints = (2 x components) + (2 x double points) + (T points) 

— (2 X enclosed regions) , (8.1) 

which can be derived by simple topological/graph-theoretic arguments. Furthermore, for 
all but one of the strips considered in this paper, the limiting curves B do not contain 
any enclosed region or T point (the only exception is the L — 4 square-lattice strip, which 
contains two enclosed regions and a single T point, as shown in Figure 4). Finally, the 
number of double points is given simply by 



#Q = mod(L-l,2) = r ^Jt^^^^f (8.2) 

10 if L IS odd 



if L is odd 

Therefore, the topological structure of the sets B can be summarized very simply: for odd 
L, the set B consists of a collection of arcs separated by extremely small gaps (these gaps 
are in fact nearly invisible in Figures 7 and 11 for L > 5); for even L, it is the same, except 
that the component touching the real axis consists of a complex arc and a small segment of 
the real axis intersecting at a double point at w = wqq. 

Furthermore, the number of endpoints (#E) and the number of components (#C) seem 
to grow with the strip width in a very regular way. In fact, for all the strips considered in 
this paper except the L = 4 square lattice, we find 

#E = 2L (8.3a) 
#C = L-mod(L-l,2) (8.3b) 

We conjecture that (8.2) and (8.3) hold for larger L as well. If this conjecture is indeed 
true, it suggests that the endpoints (whose number grows linearly with L) will become dense 
in B. This may imply in turn that the infinite-volume limiting curve B^o will be a natural 
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boundary for the infinite- volume free energy f{w), i.e., a boundary through which analytic 
continuation of f{w) is impossible. 

Finally, we have found that the endpoint Wb{L) behaves with L in a very regular way: 
namely, its real part tends logarithmically to +00 [cf. (7.12)-(7.15) and (7.49)-(7.52)] while 
its imaginary part tends to a finite value [cf. (7.16)-(7.17) and (7.53)]. This suggests that 
the infinite- volume limiting curve Boo will extend to infinity. This behavior is consistent with 
the idea of an asymptotically free theory defined around the critical point at w — +00 [24]. 

8.2 Conformal field theory of the q — > Umit 

In this paper we have obtained several new pieces of information concerning the phase 
diagram of the q ^ limit of the Potts model on the square and triangular lattices, and 
on the interpretation of this model in terms of conformal field theory (CFT). This informa- 
tion provides evidence concerning the behavior of the Potts model also at small positive q, 
confirming some conjectures of Saleur [46,47] and refuting others. 

We have found that the phase diagram for the two lattices is similar: for w > wq, the 
model is non-critical with exponential decay of correlations; for w < Wq, the model is critical 
and can be interpreted as a Berker-Kadanoff phase. Finally, at w = wq the model displays a 
first-order critical point [94] , that is, a first-order transition point for which the correlation 
length diverges as w I wq. The transition at w — wq corresponds to the — > limit of the 
antiferromagnetic transition curve of the Potts model. 

It is well-known that the Potts model in the antiferromagnetic (— 1 < ?; < 0) and un- 
physical {v < —1) regimes may display non-universal (i.e. lattice-dependent) behavior; this 
is true in particular for the zero-temperature antiferromagnetic (chromatic polynomial) limit 
V — —1 [5-7]. Since both the Berker-Kadanoff phase and the w — wo transition point lie in 
the antiferromagnetic regime w < oi the spanning-forest model, it is therefore interesting 
to assess the extent (if any) to which they display universality. 

For the Berker-Kadanoff phase we have indeed found a number of universal features. For 
both the square- and triangular-lattice models, the central charge is found to be consistent 
with the g — > limit of Saleur 's prediction (1.9) for the Berker-Kadanoff phase: 

c(Berker-Kadanoff) = -2 . (8.4) 

Furthermore, the first two thermal critical exponents appear to be given by 

xt,i = 2 (8.5a) 
xt,2 = 4 (8.5b) 

for both lattices. The first exponent agrees with the — > limit of Saleur 's prediction (1.10). 
The second exponent can presumably be derived from a more detailed analysis of the free 
fermion theory that is thought to govern the w ^ —00 limit. 

It is interesting that the results (8.4) and (8.5a) are identical to what is predicted [24] 
for the ferromagnetic fixed point at w = -|-oo. Indeed, all our results are consistent with 
the idea [24] that the fixed points lying at w — ±00 are identical and are given by a theory 
of a pair of free scalar fermions, with central charge c = —2 and marginal first thermal 
exponent Xt^i = 2. Moreover, renormalization-group computations in the four-fermion model 
describing spanning forests at g = 1/w ^ [24] show that the g = fixed point is marginally 
repulsive for g > and marginally attractive for g < 0. Therefore, the Berker-Kadanoff 
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phase a,t —oo < w < Wo renormalizes onto this free fermion theory and hence has c — —2, 
while the phase a.t Wq < w < +00 is noncritical (hence c — 0) and obeys asymptotic freedom 
as w ^ +00. More generally, there is a separatrix in the complex g-plane separating these 
two regimes: it is given by the curve Boo- In Figures 32 and 33 we have depicted the relevant 
phase diagrams, by mapping Figures 7 and 11 from the w- plane to the ^r-plane. Our best 
estimate for the separatrix B^o is given by the black dots and, more roughly, by the black 
dotted-dashed curve, which were obtained by the methods described in Sections 7.3 and 7.7. 
This phase diagram in the complex g-plane (where g is the bare coupling) is typical of what 
is expected for an asymptotically free theory, whose renormalization-group flow satisfies 

f = b,g' + hg' + ... (8.6) 

with 62 > 0. Points lying inside the cardioid B^o are attracted as / — > +cxd to the fixed point 
at g — 0, and hence belong to the complex Berker-Kadanoff phase. Points lying outside the 
cardioid Boo flow away from the fixed point at = 0, and hence belong to a different phase 
(i.e. the noncritical phase). 

Let us now turn to the w = Wo transition point, which corresponds to the g — > limit 
of the antiferromagnetic transition curve and gives information about the probable behavior 
of this curve at small positive q. For the square lattice, we have strong evidence that the 
correlation length ^(w) diverges as w I wq, so that this transition point is indeed a critical 
point. For the triangular lattice, however, the analysis of this point is difficult, since the 
exact value of Wq is unknown and we know fewer terms in the small--u; expansion. The 
following pieces of information are therefore based on the square-lattice case only: 

1) For even L we find a central charge that is consistent with the g — > limit of Saleur's 
prediction (1.6) for the antiferromagnetic critical curve: 

c = -1 . (8.7) 

For odd L, we find [cf. (7.39)] that the effective central charge is consistent with 

Ceff = 1/2 , (8.8) 

which we explain by conjecturing the existence of a twist operator with conformal weight 

/itwist = -1/16 (8.9) 

[cf. (7.40)-(7.41)]. The role of boundary-condition-changing operators in the type of CFT 
proposed by Saleur [46,47] is currently under study [155]. 

2) Our data at w > Wq and even L are consistent with the following behavior as w I 
Wo = —1/4 of the infinite- volume correlation length ^{w) and the infinite- volume free energy 
fH: 

(8.10a) 

f"(w) ~ (w-wo)'^ (8.10b) 
f'{w) ~ log(w-wo) (8.10c) 

f{w) ~ {w — Wq) log{w — Wq) . (8.10d) 

3) We find that the dominant thermal eigenvalue is consistent with 

xt,i = (8.11) 
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and hence v — 1/2. This result is in clear disagreement with Saleur's prediction xt,\ — 2 
arising from the g ^ limit of his conjecture (1.7)/(1.8); it indicates that the thermal 
operators along the antiferromagnetic critical curve need to be restudied [25]. Our data for 
the second thermal eigenvalue indicates that 

XT,2 ^ 1.25 . (8.12) 

An important open question is to obtain the same quantities for the triangular-lattice 
model and compare the to the above ones. This will provide a direct check of universality 
for these two models. 

8.3 Behavior of the free energy and its derivatives 

We have studied the behavior of the free energy and its derivatives near the transition 

point w = Wo using two complementary methods: extrapolation of small-w series expansions 
(Appendix A) and finite-size analysis of the transfer-matrix data (Sections 7.4 and 7.8). 

Series expansions give valuable information concerning the behavior as w | wq, directly in 
the infinite- volume limit. Our results show that for the square lattice /" diverges as w I wq 
with an exponent a ~ 1, while /' is most likely finite (but we cannot rule out a very weak 
divergence). For the triangular lattice, both /" and /' appear to take finite values when 
w I Wq. The finite-size data are consistent with these conclusions. 

For the regime w t Wq, by contrast, only the finite-size data give useful information. For 
the square lattice we have weak evidence for a diverging /". For the triangular lattice, we 
have clear evidence that /" and /' stay finite as w approaches Wq from below. 

These puzzling differences of behavior between the square and triangular lattices have 
implications for universality; and conversely, the sundry (and more or less plausible) univer- 
sality conjectures have implications for the interpretation of our numerical results. On the 
one hand, the analytical results of [24] show that, at least for the spanning-forest model (i.e., 
the limit q,v ^ with w — v/q fixed of the Potts model), the behavior inside the Berker- 
Kadanoff phase is universal: more precisely, at each fixed w < Wo, the spanning-forest model 
is simply the c = —2 theory of a pair of free scalar fermions, perturbed by a four-fermion 
operator that is (in this phase) marginally irrelevant. On the other hand, as noted in the 
Introduction, Saleur [46, p. 669] expects universaUty not only for the Berker-Kadanoff phase, 
but also for the critical theories forming its upper and lower boundaries. The validity of this 
latter conjecture is, however, an open question: in particular, we do not yet understand very 
well the new phase-transition curve that we have found [25] in the triangular-lattice Potts 
model, so we are unable to say for sure whether it belongs to the same universality class as 
the square-lattice antiferromagnetic critical curve, in the sense of having the same central 
charge c{q) and the same critical exponents XTi{q),XT2{(l), etc. 

Now, the behavior as w Wq is controlled by the theory at wq. Therefore, if universality 
does in fact hold (at least in the limit q ^ 0) for the upper boundary of the Berker-Kadanoff 
phase, then (barring subtleties) the square and triangular lattices should exhibit identical 
behavior as w — > Wq. In particular, our clear results on the triangular lattice (/" and /' 
staying finite as w t Wq) would suggest that the same behavior ought to occur for the square 
lattice; if so, our contrary indications for the square lattice would be an artifact of too-small 
strip widths L. Likewise, our clear results on the square lattice (/" diverging as w I Wo) 
would suggest that the same behavior ought to occur for the triangular lattice; if so, our 
contrary indications for the triangular lattice would be an artifact of too-short small-w series. 
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On the other hand, if universahty does not hold for the theory at wo, then we are unable 
at present to draw any firm conclusion about the behavior of /" as w ^ Wq on the square 
lattice, or as w J, Wo on the triangular lattice. 



A Series analysis 

In this appendix we shall analyze the small-w series expansions displayed in Table 2. For 
simplicity, we shall denote the series expansions for the free energy and its derivatives as 
follows: 

oo 

f{w) = Y.fkw" r.{w-wof-'' (A.la) 

k=l 

oo oo 

f{w) = Y.^kw'' -Y.^k + l)fk+iw^ ^{w-wof-^ (A.lb) 

fe=0 k=0 

oo oo 

f"{w) = ^c,i«'^ = ^(fc + l)(A; + 2)/fe+2«;'~(«;-«;o)"" (A.lc) 

fe=0 fe=0 

The rightmost side of each equation in (A.l) shows the expected leading asymptotic behavior 
for w near wq (assuming, for simplicity, that no multiplicative logarithms occur). We have 
used the standard notation a for the critical exponent associated to the "specific heat" f"{w). 

Let us start with the square-lattice case, which is the one that leads to the most interesting 
results, both because the series is twice as long as for the triangular lattice and because the 
putatively exact critical point wq — —1/4 is known. The sign of the coefficients of the free- 
energy expansion has a clear alternating pattern (at least up to k = 47): fk ~ (—1)*^"'"^. This 
behavior is a clear sign that the singularity nearest to the origin lies on the negative w-axis. 

A first (and very rough) approach can be obtained using the ratio method [158, and 
references therein]. Let us suppose that F{w) is a function for which we can obtain the 
series expansion F{w) = YlT=i fk'^'^i ^-nd whose behavior close to w — wq is of the form 

Fiw) ^ a(i-—\ +B , (A.2) 

Estimates of the critical point Wq and the critical exponent A can be obtained from the series 
coefficients {fk} by computing the ratios 

rk ^ 4^. (A.3) 

Jfc-i 

In absence of any competing singularity, a plot of versus 1/ /c is expected to be a straight 
line taking the value l/woail/k — and having slope — (A-|- l)/wo. In Figure 34(a) we show 
the ratio rfc versus 1/k for the free energy and its first and second derivatives. In addition 
to the expected linear behavior in 1/k, we find a strongly oscillatory behavior that makes 
the analysis rather difficult. It is apparent from Figure 34(a) that the limit of the three 

^^It is not clear from Figure 34(a) what is the period of these oscillations (if indeed they are periodic at 
all), but it is clear that the period is not 2. 
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sequences as — > is close to the expected value wo(sq) = —1/4. As a matter of fact, the 
value of Wq can be estimated from the sequence [158] 

{k + e)rk - {k + e- l)rk-i = — [l + 0{l/k^)] (A.4) 

Wo 

where e is any small /c-shift. The oscillations are smaller for these new sequences, and we 
conclude that wq — —0.250(1). We can therefore assume that wq — —1/4 exactly, and try 
to estimate the exponent A using the sequence [158] 

Afc = k-1- kwoTk = A [1 + 0{l/k)] (A.5) 

with Wq fixed to its theoretical value. The estimates A^ as a function of 1/A; for /, /' and /" are 
displayed in Figure 34(b). Again, wc observe oscillations in addition to the expected linear 
behavior in 1/k. Rough estimates of A can be obtained from the corresponding sequences 

2- a = 1.10(8) (A.6a) 
1-a = 0.10(3) (A.6b) 
a = 0.90(3) (A.6c) 

These results are close to the expected value a = 1 for a first-order critical point [94], 
although the error bars cannot be taken at face value. 

A more quantitative study can be performed using differential approximants [158, and 
references therein]. The i^-th order differential approximant to a function F{w) is built 
as follows: we choose polynomials^Qo, ■ ■ ■ , Qk and P of degrees A^o, A^i, ■ ■ ■ , Nk and L, 
respectively, so that the solution F of the inhomogeneous linear differential equation 

agrees with the first coefficients of the series F{w) = Yl!k=ofkW^- The resulting equations 
can be solved if the number of unknown coefficients in the polynomials is smaller than the 
order of the available scries expansion for F. The singularities of the function F[w) are 
located at the zeros {w^} of the polynomial Qk (along with w = and w — ±oo). The 
critical exponent associated to a simple zero of Qk is given by 

A,= ^-l-%^. (A.8) 

WiQK\Pd 

When wg^ is a multiple zero of Qk, this formula should be modified; details can be found in 

reference [158]. 

In practice, we have used a modified version of the program newgrqd.f described in 
reference [158] to obtain the differential approximants. Our program uses Mathematica to 
obtain the polynomials Qj exactly (i.e., with exact rational coefficients), and then MPSolve 
[159, 160] to compute the Nk zeros of Qk to arbitrarily high precision (100 digits in our 

case). We have computed all the differential approximants of first and second order (i.e., 
K = 1,2) satisfying |Aj ~ ^jl ^ 1 ^^d using at least 36 coefficients of the corresponding 
series. For each zero w^ of Qk-, we have computed the corresponding critical exponents A^. 
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Our procedure has the advantage, over the double-precision Fortran program newgrqd.f, 
that roundoff errors are under control. We have checked the accuracy of our results with the 
help of an independent program by Y. Chan, A.J. Guttmann and A. Rechnitzer written in 
C++. 

Once all the approximants have been computed, we need to discard the defective ones. 
We consider an approximant to be non-defective if there is a real zero of Qk sufficiently 
near to the expected value Wo(sq) = — 1/4 and this zero is sufficiently well separated from 
all other zeros of Qk- In practice, we asked that the zero satisfy —0.251 < < —0.249 
for /", -0.256 < < -0.248 for /', and -0.257 < < -0.247 for /, and that no other 
singularity should appear in the region —0.3 < Rew < 0, —0.05 < Imw < 0.05. This is 
essentially what the program tabul.f (described in reference [158]) does. 

The non-defective approximants for the square-lattice spanning-forest free energy and its 
first and second derivatives are displayed in Figures 35(a)-(c), respectively. It is interesting 
to note that the estimates are not scattered uniformly over the corresponding plots; rather, 
they tend to accumulate along certain curves. Unbiased estimates for wq and A are obtained 
by averaging over the data points; crude estimates of their precision can be obtained by 
defining the error bars to be one standard deviation of the corresponding data distribution. 
(But we emphasize that these are not true statistical errors as in a Monte Carlo simulation, 
and the error bars should not be taken too seriously.) 

As we know the exact value of wo(sq) = —1/4, we can also obtain biased estimates of 
the critical exponents A. The idea is to construct differential approximants for which wq is 
an exact simple zero of the polynomial Qk- We have analyzed these approximants following 
the same criteria as for the non-biased ones. 

Let us begin with the second derivative of the free energy, f"{w). The data points in 
Figure 35(c) fall over a narrow interval in both axes, and we see no significant differences 
between the first- and second-order approximants. Our preferred unbiased estimates come 
from merging both types of approximants: 

Wq = -0.25011 + 0.00016 (A.9a) 
a = 0.933 + 0.030 (A.9b) 

The estimate for the critical point agrees very well with the theoretical prediction wo(sq) = 
— 1/4. The estimate for the critical exponent is close to the expected result a = 1; but it 
is barely compatible with it within errors (if one takes the error bars literally). The biased 
estimates for a are also very similar ior K — 1 and K — 2,so that once again our preferred 
estimate comes from merging both data sets: 

a = 0.913 + 0.021 (A.IO) 

This estimate is four standard deviations away from the expected value a = 1. 

The analysis of the first derivative f'{w) is less clear: the spread of the data points in 
Figure 35(b) is much larger than it was for the second derivative. We observe that most of 
the first-order approximants give a singularity with Wi < —1/4, while most of the second- 
order ones lie at somewhat higher values of wq (the latter approximants also have smaller 
dispersion). The estimates obtained from the K — 2 approximants are 

Wo = -0.25004 + 0.00063 (A.lla) 
I- a = 0.10 + 0.22 (A.llb) 
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while those coming from merging both data sets are 



Wo = -0.2506 ±0.0011 (A.12a) 

1- a = -0.08 ±0.31 (A.12b) 

We observe that the error bars are 4-7 times larger than for the second-derivative series. 
The estimates of wq are still compatible within errors with the expected resuh wq = —1/4. 
The estimates of a are compatible with the expected value a = 1, but the error bars are so 
large that we are unable to tell whether the first derivative diverges or not at w — Wq (i.e., 
whether 1 — a is negative or positive). If we look at the biased estimates of 1 — a, we see 
that the dispersion of the second-order approximants is smaller than for the first-order ones. 
Thus, we take as our preferred estimate the value coming from the K — 2 approximants: 

1- a = 0.104 ±0.040 (A.13) 

This value suggests that the first derivative /' does not diverge ai w = wq; but this conclusion 
cannot be taken too seriously as the signal is only 2.5 times larger than the alleged error bar. 

We conclude this analysis with the free energy [Figure 35(a)]. The dispersion of the data 
is even larger than for the first derivative, so we should expect even larger error bars. The 
K = 1 approximants yield estimates of wq quite a bit below —1/4, and estimates of 2 — a 
that are rather near to zero than to 1. Therefore, our preferred estimate comes from the 
K — 2 data set: 

Wo = -0.2501 ±0.0011 (A.14a) 

2- a = 1.07 ±0.39 (A.14b) 

The estimates are compatible with the expected values, but the error bars are huge. The 
analysis of the biased estimates for 2 — o; is similar and our preferred estimate comes from 
the K — 2 approximants: 

2- a = 1.063 ±0.095 , (A.15) 

which is again compatible within errors with the expected value a = 1. 

Let us now summarize the results of this analysis. We have found that the the unbiased 
estimates of wq are in all cases compatible within errors with the expected value wq — —1/4: 

( -0.25011 ±0.00016 for /" 
Wo ^ I -0.25004 ±0.00063 for/' (A.16) 
( -0.2501 ±0.0011 for/ 

This result is an independent confirmation that Wo(sq) = —1/4. The unbiased estimates of 
the critical exponent a are given by 

( 0.933 ±0.030 for/" 

a = I 0.90 ±0.22 for/' (A. 17) 

[ 0.93 ±0.39 for/ 

while the biased estimates are 

0.913 ± 0.022 for /" 

a = { 0.896 ±0.040 for/' (A.18) 

0.937 ± 0.095 for / 
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The foregoing estimates are compatible among themselves within the quoted errors. However, 
the estimates coming from /" are at least a factor of 2 more precise than those coming 
from /' or /. It is intriguing that all the estimates of a are consistently smaller than 
the expected value a = 1. Furthermore, the difference between the above estimates and 
the theoretical prediction is significant, if one takes the error bars seriously. Perhaps this 
indicates a multiplicative logarithmic correction to the leading behavior a — \. 

The story for the triangular lattice is rather short. Prom Table 2 we observe that the 
coefficients have a complicated sign pattern. This is an indication that the dominant 
singularity (or singularities) is/are complex. This fact makes the ratio method useless. 
Furthermore, the differential-approximant method does not give any sensible estimate of Wq 
(or of the corresponding critical exponents) . We computed the differential approximants to 
the second derivative /" (which is a •priori the most favorable observable) using at least 20 
series coefficients and with the same constraints on the A^^ as for the square lattice. We found 
that the zero closest to «;o(tri) a; —0.1753 varies over a range vastly wider than for the square 
lattice (by two orders of magnitude) . We were therefore obliged to take a very lenient view 
of what constitutes a non-defective approximant: we considered an approximant to be non- 
defective if there is a real zero in the interval —0.3 <Wii< —0.1 and there is no other zero 
in the region —0.5 < Retu < and —0.05 <\va.w < 0.05. The non-defective approximants 
are displayed in Figure 35(d). We find that, even with this lenient definition, there are many 
fewer non-defective appoximants than for the square-lattice case (45 ioY K — \ and 5 for 
K — 1^ compared to 91 and 164, respectively). The non-defective approximants seem to 
accumulate along a curve; but the low density of zeros prevents us from drawing any reliable 
conclusion from these data. 
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L 


SqCyl(L) 


TriCyl'(L) 


SqFree(L) 


TriFree(L) 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


3 


3 


3 


4 


5 


4 


6 


6 


10 


14 


5 


10 


10 


26 


42 


6 


24 


28 


76 


132 


7 


49 


63 


232 


429 


8 


130 


190 


750 


1430 


9 


336 


546 


2494 


4862 


10 


980 


1708 


8524 


16796 


11 


2904 


5346 


29624 


58786 


12 


9176 


17428 


104468 


208012 


13 


29432 


57148 


372308 


742900 


14 


97356 


191280 


1338936 


2674440 


15 


326399 


646363 


4850640 


9694845 


16 


1111770 


2210670 


17685270 


35357670 


17 


3825238 


7626166 


64834550 


129644790 


18 


13293456 


26538292 


238843660 


477638700 


19 


46553116 


93013854 


883677784 


1767263190 


20 


164200028 


328215300 


3282152588 


6564120420 



Table 1: Dimensions of the transfer matrices for square- and triangular-lattice strips of 
width L with cylindrical and free boundary conditions. TriFree(L) gives the total number 
of non-crossing partitions of the set {1,2,..., L}. The other three columns give the number 
of equivalence classes of non-crossing partitions of {1,2, ... ,L} modulo some symmetries: 
SqFree(L) [reflection], TriCyl'(L) [translation], and SqCyl(L) [translation and reflection]. 
We show in boldface the dimensions of the largest transfer matrices that we have used in 
symbolic form in the computations reported in this paper. 
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k 


fe/fc(sq) 


fc/fc(tri) 


1 


2 


3 


2 


-2 


-3 


3 


2 


-3 


4 


-6 


9 


5 


22 


3 


6 


-74 


-33 


7 


254 


3 


8 


-918 


393 


9 


3422 


-2325 


10 


-12862 


9327 


11 


48138 


-24483 


12 


-178530 


1815 


13 


655826 


458253 


14 


-2391370 


-3497133 


15 


8688262 


17900097 


16 


-31600918 


-70905543 


17 


115606190 


209837565 


18 


-426864494 


-314875455 


19 


1593065490 


-1242857637 


20 


-6004037966 


13767148419 


21 


22795625582 


-73591894407 


22 


-86925982926 


264618912819 


23 


332053760646 


-455046987303 


24 


-1268578680714 


-2402029948737 


25 


4844247521322 


31594277221653 


26 


-18494593884938 




27 


70644561464090 




28 


-270190727926594 




29 


1035346222307838 




30 


-3976144389096514 




31 


15304261265448246 




32 


-59027790560689238 




33 


228068574553887894 




34 


-882464779625379526 




35 


3418403165602360314 




36 


-13253811969559767270 




37 


51425605359158653378 




38 


- 199663639129405278414 




39 


775682274057446798018 




40 


-3015326560156376960998 




41 


11728909932236608346846 




42 


-45652065541079598767758 




43 


177805114097058031764786 




44 


-692953582585445674377902 




45 


2702276554574318555870842 




46 


- 10544061987158176469650990 




47 


41164706135505931628292550 





Table 2: Small-w (high-temperature) expansions for the bulk free energy f{C,w) — 
YlT=i fk{J^)w^ for the square and triangular lattices. For simplicity, we present here, in- 
stead of the coefficient fk{^), the product kfk{jC), which is always an integer. 



72 



Tvue 


7/? 


J 1 


J 2 


J A 


J 4: 


Jo 




2 




01005 


00349 


00405 


00020 
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Table 3: Fits of the large-w series for the square lattice. We show the estimates of the fits 
of the (infinite-volume) free-energy data (for large values of \w\ > 2) to the Ansatz (3.19) 
with /cmax = 5. The column "Type" shows whether /i has been left free in the fit ("Non- 
Biased") or it has been fixed to its conjectured value 1/8 ("Biased"). Non-biased fits are 
based on data points with |w| = Wmin, • • • , Wmin + 4; biased fits are based on data points with 

\w\ = Wmin, • • • , W^min + 3. 
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Table 4: Fits of the large-w series for the triangular lattice. We show the estimates of the 
fits of the (infinite- volume) free-energy data (for large values of lu^j > 2) to the Ansatz (3.19) 
with /Cmax = 5. The column "Type" shows whether /i has been left free in the fit ("Non- 
Biased") or it has been fixed to its conjectured value 1/12 ("Biased"). Non-biased fits are 
based on data points with \w\ = Wmin, • • • , Wmin + 4; biased fits are based on data points with 

\w\ = Wmin, . . . , Wmin + 3. 
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Table 5: Characteristics of the hmiting curves B for square-lattice strips with cylindrical 
boundary conditions. For each width L, we show the number of connected components 
(#C), the number of endpoints (#E), the number of multiple points (#Q), and the Wq and 
wb values defined in the text. The numbers marked with a t are only lower bounds on the 
exact values. 
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Table 6: Characteristics of the limiting curves for triangular- lattice strips with cylindrical 
boundary conditions. For each width L, we show the number of connected parts (#C), the 
number of endpoints (#E), the number of multiple points (#Q), and the Wq values defined 
in the text. The numbers marked with a t are only lower bounds on the exact values. 
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Table 7: Extrapolated limiting curve B for the square and triangular lattices. For each 
lattice and each fixed value of Imw, we show the extrapolated value of Rew together with 
its subjective error bar. 
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Table 8: Fits for wq on the square lattice, using the Ansatz woj{L) = wq + AL ^ (j = 
+,—,Q). Each fit is based on data points with L = L^i^, L^i^ + 2, Lmin + 4. 
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Table 9: Fits for wq on the square lattice, using the Ansatz Woj{L) — wo + A2L ^ + A4L 
{j — +, —,Q). Each fit is based on data points with L — Lmin, -^min + 2, 1/min + 4. 
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(b) 

Table 10: Fits for Wb on the square lattice, (a) For the real part, we use the Ansatze 
Re Wb{L) = Cl log L + C2 with Ci and C2 free, and Rc Wb{L) = Cy log L -f C3 log log L + C2 with 
fixed Cl = 3/(27r) and C3 = — l/(27r) [see (7.14) and surrounding text], (b) For the imaginary 
part, we first use Im u'b(L) = Imws + AL~^, and then we perform the fit with A = 1/2. 
Two-parameter fits are based on data points with L = Lmin, -^min + 1; three-parameter fits 
are based on data points with L — Lmm, -^min + 1, -^min + 2. 
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Table 11: Free energies /^(ty) [cf. (2.32)] for square-lattice strips of width L and cylindrical 
boundary conditions, for w = — 1, — 1/4, +1. The row labelled oo shows the exact value 
(7.29) in the thermodynamic limit L — oo. 



79 



-^min 


Re/(sq, w = -1/4) 


A{w = -1/4) 


A{w = -1/4) 


3 


-0.6026552580 


0.3514454435 


2.1843042660 


5 


-0.6029925573 


0.3092573584 


2.0851083253 


7 


-0.6030693569 


0.2909915798 


2.0464952763 


9 


-0.6030901419 


0.2826853879 


2.0304080199 


11 


-0.6030977473 


0.2779994113 


2.0219793544 


2 


-0.6018136184 


-1.0158505843 


2.1108840738 


4 


-0.6036426651 


-1.1575753432 


2.2297430523 


6 


—0.6034714802 


— 1.1161060082 


2.2049156257 


Q 
O 


— U.DUdZobioiU 


— i.Uozi4U00o/ 


z. ioyo4zozo / 


10 


-0.6031958219 


-0.9688529236 


2.1266061029 


oo 


-0.6031055757 








-^min 


Re /(sq, w = —1) 


A{w = -1) 


A{w = -1) 


2 


1.0366588370 


-2.6329852431 


2.4336613152 


3 


1.0442846989 


-2.3545416365 


2.2944945872 


4 


1.0467035384 


-2.1844090010 


2.2227750669 


5 


1.0477001325 


-2.0688028151 


2.1789288156 


6 


1.0481526051 


-1.9898682648 


2.1511947863 


7 


1.0483784334 


-1.9340682142 


2.1327570133 


8 


1.0485022480 


-1.8926046519 


2.1197454120 


9 


1.0485758832 


-1.8603185323 


2.1100523544 


10 


1.0486226265 


-1.8342219236 


2.1025132464 


11 


1.0486538611 


-1.8125199233 


2.0964523136 


12 


1.0486756030 


-1.7940748336 


2.0914535840 


13 


1.0486912486 


-1.7781266664 


2.0872465655 


14 


1.0487028225 


-1.7641447178 


2.0836471061 



-^min 


Re/(sq, w — 1) 


A{w = 1) 


A{w = 1) 


2 


1.3107968081 


-0, 


.7156482298 


2.4689333981 


3 


1.3093320071 


-0, 


.7881463786 


2.5852771137 


4 


1.3089601334 


-0, 


.8429753668 


2.6461540521 


5 


1.3087423951 


-0, 


.9168496742 


2.7098055866 


6 


1.3085813800 


-1 


.0307347078 


2.7878824576 


7 


1.3084637131 


-1, 


.1971341766 


2.8782847595 


8 


1.3083793230 


-1, 


.4306376635 


2.9776440535 


9 


1.3083189900 


-1, 


.7531876469 


3.0836969829 


10 


1.3082756267 


-2, 


.1979960385 


3.1951197122 


11 


1.3082442022 


-2, 


.8146824223 


3.3111199451 


12 


1.3082212320 


-3 


.6772081777 


3.4311777392 


13 


1.3082043053 


-4, 


.8963374373 


3.5549124381 


14 


1.3081917409 


-6, 


.6392455365 


3.6820180432 



Tabic 12: Fits of Re/i:(w) on square-lattice strips to the Ansatz (7.27). For each value of w 
(+1, —1, —1/4), we show the estimates of Re f{w), A{w) and A{w) obtained by fitting three 
consecutive data points with L > Lmin- This means Lminj -^min + 1) -^min + 2 for w = ±1, but 
-f'min, -^min + 2, Lmin + 4 for w = —1/4. The row labelled oo shows the known exact value 
(7.29) in the infinite-volume limit. 
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uooo^ooiyi 


-3 
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± 
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-2 
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5 


1 
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-3 
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1 


u^y-yiuuzooz 


-2 


8493574206 
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o^oz / yuoO'y: 


6 
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1 


nc;n849n'?i 4 

UOUo^ZUO 


-3 


0834112435 


1 
1 


u^yziu^oio 


-2 


8185453467 


9 
— z 


o ( ouoyu ( zi 


7 


1 


0501282194 


-3 


0166105063 


1 


0490891975 


-2 


7923749159 


-3 


2583728412 


8 


1 


0497177416 


-2 


9664373949 


1 


0490078464 


-2 


7698463825 


-3 


6800970526 


9 


1 


0494621793 


-2 


9269022644 


1 


0489513545 


-2 


7503180147 


-4 


1376812436 


10 


1 


0492933116 


-2 


8946509060 


1 


0489109249 


-2 


7332535305 


-4 


6268794806 


11 


1 


0491764712 


-2 


8676499122 


1 


0488811967 


-2 


7182077146 


-5 


1448634916 


12 


1 


0490926063 


-2 


8445853935 


1 


0488588065 


-2 


7048231509 


-5 


6897520658 


13 


1 


0490305778 


-2 


8245646824 


1 


0488415841 


-2 


6928174237 


-6 


2602275770 


14 


1 


0489835393 


-2 


8069566713 


1 


0488280902 


-2 


6819676896 


-6 


8553060059 


15 


1 


0489471059 


-2 


7913005767 















Tabic 13: Fits of Kc J'l^w = —1) on square-lattice strips to the Ansatze (7.34)/(7.35). We 
show the estimates of /(sq, w = —1), c{w = —1), and A{w = —1) obtained by fitting two or 
three consecutive data points with L > Lmin- 



min 


f(^sq,w= -1) 


c 


(t« = -l) 


A{w = -1) 


B{w = -1) 


2 


1 


0481659817 


-1 


8395026016 


-0 


2962391773 


-0 


8240960988 


3 


1 


0486170833 


-2 


0667855660 


-0 


1728865530 


-0 


7094169716 


4 


1 


0487861877 


-2 


2333779695 


-0 


0629303046 


-0 


6281602055 


5 


1 


0487985937 


-2 


2540969965 


-0 


0473045479 


-0 


6186355145 


6 


1 


0487835568 


-2 


2152906621 


-0 


0796357366 


-0 


6352167690 


7 


1 


0487715509 


-2 


1704302257 


-0 


1200678725 


-0 


6528623630 


8 


1 


0487645572 


-2 


1343787723 


-0 


1547235661 


-0 


6658127726 


9 


1 


0487606108 


-2 


1073401566 


-0 


1821652813 


-0 


6746153901 


10 


1 


0487582804 


-2 


0867372527 


-0 


2040744739 


-0 


6806453142 


11 


1 


0487568178 


-2 


0704483013 


-0 


2221173626 


-0 


6848921285 


12 


1 


0487558491 


-2 


0571283328 


-0 


2374141788 


-0 


6879523513 


13 


1 


0487551802 


-2 


0459596860 


-0 


2506617512 


-0 


6901835157 



Table 14: Fits of Kc Jl^w = —1) on square-lattice strips to the Ansatz (7.36). We show the 
estimates of /(sq, w = —1), c{w — —1), A{w = —1), and B{w = —1) obtained by fitting 
four consecutive data points with L > L^in- 
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fL = f + (c7r/6)L-2 


/l = / + (C7r/6)L-^ + AL-^ 




c(w = -1/4) 


r(vi — —1 /4-) 

(.JJ ± / t: J 


A (on — _1 /4^ 


2 


-1.7867317647 






3 


0.5559211806 






4 


-1.6242104680 




—1 R47197R97n 

X .0'4: / XZ(/UZ(/U 


5 


0.5204270917 


U.OUlOlO 1 oou 


n 9^01 971 zLl n 


6 


-1.5017190128 




—9 7QRR9^9^9^ 


7 


0.5110679214 


0.5017424517 


0.2392574198 


8 


-1.4367580104 


-1.3338122127 


-3.4497467921 


9 


0.5073837852 


0.5016573071 


0.2428685331 


10 


-1.3996975232 


-1.3212758647 


-4.1061484386 


11 


0.5054907346 


0.5015180462 


0.2516914702 


12 


-1.3757353498 


-1.3111049796 


-4.8730151126 


13 


0.5043623971 


0.5013945821 


0.2626165842 


14 


-1.3585885169 






15 


0.5036237409 







Table 15: Fits of Refiiw = —1/4) on square-lattice strips to the Ansatze (7.34)/(7.35), 
using the exact value (7.29) of f{sq,w — —1/4) in both Ansatze. We show the estimates of 
c{w — —1) and A{w — —1) obtained by fitting one or two consecutive data points of the 
same parity. 









+ (27rxT,i)/L 


{27TXT,l)/L 


L 


^min 




Xt,\ 


3 


0.8447211436 


3 


-0.04367 0.42418 


0.40332 


5 


0.4893646433 


5 


-0.03794 0.41962 


0.38942 


7 


0.3387055144 


7 


-0.02746 0.40793 


0.37735 


9 


0.2573363732 






0.36861 


2n 












Table 16: Inverse correlation length ^f^(-L) for the square lattice w — —1/4. We also 
show the results of fitting the odd- width data to the Ansatz (7.44), either with ^j""^ variable 
(columns 4-5) or ^ = fixed (column 6). The last row shows that for all even widths L, 
we have ii^{L) = 0. 
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^2\L) 




^2-' + i27rXT,2)/L 


(27rxT2)/L 


L 






Xt,2 


Xt,2 


3 


2.1006084199 


3 


0.51802 


0.75563 


1.00297 


5 


1.4675739282 


5 


0.19309 


1.01421 


1.16786 


7 


1.1034345588 


7 


0.06489 


1.15703 


1.22932 


9 


0.8726454559 


9 






1.24997 


4 


0.7229219112 


4 


0.09380 


0.40051 


0.46023 


6 


0.5132143892 


6 


-0.01991 


0.50910 


0.49008 


8 


0.3799322185 


8 


-0.04373 


0.53942 


0.48374 


10 


0.2952007610 


10 






0.46983 



Table 17: Inverse correlation length ,^2^ (L) for the square lattice at w = —1/4. We also 
show the results of fitting the odd- width data to the Ansatz (7.44), either with variable 
(columns 4-5) or = fixed (column 6). The fits are done separately for L odd and L 
even. 






Lm'm 


Wo 


A 


A 




2 


-0.176711 


-0.679548 


3.212914 




4 


-0.175667 


-0.303430 


2.541774 




6 


-0.175334 


-0.183161 


2.204638 




2 


-0.178005 


-0.371851 


2.978114 




3 


-0.175706 


-0.222577 


2.373516 




4 


-0.175505 


-0.205877 


2.299951 




5 


-0.175380 


-0.189033 


2.231829 




6 


-0.175345 


-0.182054 


2.205216 




7 


-0.175323 


-0.175936 


2.183217 




8 


-0.175311 


-0.171338 


2.167416 




2 


-0.173481 


-0.080930 


1.516018 




4 


-0.174992 


-0.120997 


1.925676 




6 


-0.175316 


-0.171962 


2.171007 



Table 18: Fits for wq on the triangular lattice, using the Ansatz woj{L) = wo + AL~ 
[j = +, — , Q). Each fit for wo± is based on data points with L = L^^j^, Lmin + 2, Lmin + 4; 
fits for wqq are based on data points with L = Lmin, L,^in + 1, -^min + 2. 
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o 


-^min 


Wo 


A2 


A, 


Wo- 


2 


-0.175684 


-0.091446 


-0.823270 




4 


-0.175347 


-0.108955 


-0.629315 




6 


-0.175247 


-0.118958 


-0.398856 




2 


-0.176194 


-0.101075 


-0.379755 


3 


-0.175074 


-0.129083 


-0.218428 




4 


-0.175180 


-0.124719 


-0.260998 




5 


-0.175202 


-0.123381 


-0.280737 




6 


-0.175230 


-0.121007 


-0.330011 




7 


-0.175245 


-0.119307 


-0.377189 




8 


-0.175255 


-0.117857 


-0.429037 




2 


-0.174942 


-0.144129 


0.147137 




4 


-0.175066 


-0.137675 


0.075637 




6 


-0.175242 


-0.120065 


-0.330093 



Table 19: Fits for wq on the triangular lattice, using the Ansatz woj{L) = wq+A2L~'^+A. 
{j = +, — , Q). Each fit for wq^. is based on data points with L = L-^m, -^min + 2, Lmin 
and those for wqq are based on data points with L = Lmm, -f/min + 1, -^'min + 2. 
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= Cl log L + C2 


RewB{L) = Cl log L + C3 log log L + C2 




Cl 






2 


0.265812 


-0.302108 


-0.342615 


3 


0.239695 


-0.273416 


-0.304290 


4 


0.233643 


-0.265025 


-0.293266 


5 


0.232791 


-0.263654 


-0.288928 


6 


0.233405 


-0.264754 


-0.286884 


7 


0.234402 


-0.266694 


-0.285815 


8 


0.235446 


-0.268865 


-0.285226 


9 


0.236434 


-0.271036 


-0.284900 


10 






-0.281730 



(a) 





Im Wb{L) 




-^min 


Im 


A 


A 


2 


0.369334 - 


-0.205690 


0.810857 


3 


0.398666 - 


-0.215786 


0.582965 


4 


0.401549 - 


-0.216752 


0.564881 


5 


0.398620 - 


-0.216574 


0.585571 


6 


0.396188 - 


0.217407 


0.605898 


7 


0.394800 - 


-0.218515 


0.619295 


8 


0.394137 - 


-0.219347 


0.626437 



(b) 

Table 20: Fits for wb on the triangular lattice, (a) For the real part, we use the Ansatze 
KewB{L) — cilogL + C2 with ci and C2 free, and RewB(L) = cilogL + caloglogL + C2 
with fixed Ci = ■\/3/(27r) and C3 — — l/(2\/37r) [see (7.50) and surrounding text], (b) For 
the imaginary part, we use Imw^l-^) = Imw^ + AL^^. Two-parameter fits are based on 
data points with L = L^m, L^m + 1; three-parameter fits are based on data points with 

L -^min) -^min ~l~ f ; -^min ~l~ 2. 
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L 


RefL{w = -l) 




2 


1.2070144097 


1.5619770404 


3 


1.3954889413 


1.6442216134 


4 


1.4579588323 


1.6714852756 


5 


1.4861038633 


1.6835902358 


6 


1.5011384591 


1.6899375329 


7 


1.5100964902 


1.6936439914 


8 


1.5158576372 


1.6959785719 


9 


1.5197784405 


1.6975342156 


10 


1.5225658002 


1.6986169734 


11 


1.5246173608 


1.6993972288 


12 


1.5261706619 


1.6999756628 


13 


1.5273746626 


1.7004147101 


14 




1.7007546649 


15 




1.7010224200 



Table 21: Values for the quantities Re/i(-u;) (2.32) on triangular-lattice strips of width L, 
infinitely length and cylindrical boundary conditions for w = ±1. 
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Rc /(tri, w = -1) A{w = -1) A{w = -1) 



2 1.5308739278 

3 1.5330939357 

4 1.5337115052 

5 1.5339268515 

6 1.5340208601 

7 1.5340693211 

8 1.5340972882 

9 1.5341147600 

10 1.5341263362 

11 1.5341343629 



-1.4384536794 

-1.3873709566 
-1.3570788682 
-1.3384597021 
-1.3256428931 
-1.3159780873 
-1.3082485775 
-1.3018238947 
-1.2963375714 
-1.2915577847 



2.1510787148 

2.1033613397 
2.0815319808 
2.0701440440 
2.0631760539 
2.0583757793 
2.0548049206 
2.0520093666 
2.0497401388 
2.0478477572 



-^min 


Re/(tri,'u; = 1) 


A{w = 1) 




A{w = 1) 


2 


1 


7033149783 


-0 


6276028679 


2 


1507031799 


3 


1 


7031440137 


-0 


6317621745 


2 


1593529516 


4 


1 


7030024774 


-0 


6399100527 


2 


1718288870 


5 


1 


7028803279 


-0 


6534158534 


2 


1887282099 


6 


1 


7027905253 


-0 


6704079915 


2 


2069423251 


7 


1 


7027253709 


-0 


6899461332 


2 


2253789191 


8 


1 


7026770177 


-0 


7116977010 


2 


2437647939 


9 


1 


7026401520 


-0 


7355754454 


2 


2620579091 


10 


1 


7026113515 


-0 


7615959157 


2 


2802754383 


11 


1 


7025883787 


-0 


7898291690 


2 


2984470780 


12 


1 


7025697290 


-0 


8203795382 


2 


3166023025 


13 


1 


7025543610 


-0 


8533775965 


2 


3347671232 



Table 22: Fits of Re/i(u') on triangular-lattice strips to the Ansatz (7.27). For each value 
of w = ±1, wc show the estimators of Re /(tri, w), A(w), and A(w) obtained by fitting three 
consecutive data points with L > Lmin- 







fL = f+ (CTI 


V3/12)L-2 


fL = f 


+ (c7rV3/12)L-2 


+ AL-^ 


-^'min 


/(tri,w = -1) 


c{w = —1) 


/(tri, w = —1) 


c{w = —1) 


A{w = -1) 


2 


1 


5462685665 


-2.9926499123 


1.5356134959 


-2.6871786001 


-0.3835825420 


3 


1 


5382772635 


-2.8340398079 


1.5349369676 


-2.6498796414 


-0.4810026139 


4 


1 


5361394741 


-2.7586078037 


1.5346428114 


-2.6232826442 


-0.5986650932 


5 


1 


5353079949 


-2.7127659475 


1.5344818361 


-2.6016275649 


-0.7435428479 


6 


1 


5349033457 


-2.6806402993 


1.5343868622 


-2.5838245327 


-0.9110768209 


7 


1 


5346773842 


-2.6562227934 


1.5343273821 


-2.5690020579 


-1.0976063455 


8 


1 


5345391118 


-2.6367069990 


1.5342881192 


-2.5564469074 


-1.3011455956 


9 


1 


5344487544 


-2.6205664178 


1.5342610311 


-2.5456343907 


-1.5205586626 


10 


1 


5343866971 


-2.6068808277 


1.5342416477 


-2.5361873825 


-1.7550985119 


11 


1 


5343423765 


-2.5950541591 


1.5342273504 


-2.5278319025 


-2.0042150985 


12 


1 


5343097063 


-2.5846792515 









Table 23: Fits of Re/i(w = —1) on triangular-lattice strips to the Ansatze (7.34)/(7.35). 
We show the estimates of /(sq, w = —1), c{w = —1), and A{w = —1) obtained by fitting 
two consecutive data points with L > Lmin- 
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-^min I /(tri,W = -1) 

2 1.5342817466 

3 1.5342534868 

4 1.5342099124 

5 1.5341878304 

6 1.5341775460 

7 1.5341723071 

8 1.5341693702 

9 1.5341675995 
10 1.5341664720 



c{w = —1) 
-2.2815122109 
-2.2650710611 
-2.2155031379 
-2.1729190238 
-2.1422717553 
-2.1196678362 
-2.1021864612 
-2.0881779553 
-2.0766679965 



A{w = -1) 
-0.0987125628 
-0.1064401573 
-0.1347733865 
-0.1625865218 
-0.1846992056 
-0.2023424141 
-0.2168955988 
-0.2292081839 
-0.2398080978 



B{w = -1) 
-0.2264608061 
-0.2336450376 
-0.2545830616 
-0.2715365779 
-0.2828772076 
-0.2905771445 
-0.2960154919 
-0.2999650620 
-0.3028824071 



Table 24: Fits of Re fL{w = —1) on triangular-lattice strips to the Ansatz (7.36). We show 
the estimates of /(tri, w — —1), c{w — —1), A{w — —1), and B{w — —1) obtained by fitting 
four consecutive data points with L > Lmin- 
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Square— lattice Potts model 



> 




Figure 1: Phase diagram for the square-lattice Potts model in the (g, f)-plane. The solid 
curves show the ferromagnetic {v > 0) and antiferromagnetic {v < 0) phase-transition curves. 
The dots • (resp. circles o) indicate the known second-order (resp. first-order) transition 
values Vc{q) for integer q. The dashed and dot-dashed curves represent (1.2_) and (1.3_), 
respectively. The hatched region corresponds to the conjectured Berker-Kadanoff phase. 
The horizontal dashed line corresponds to infinite temperature (f = 0), and the horizontal 
dotted line corresponds to the zero-temperature antiferromagnet {v = —1). The pink vertical 
lines show the Beraha numbers q = 4cos'^{TT/n) (n = 2,3, . . .); at these values the Berker- 
Kadanoff phase is not defined. 
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Triangular— lattice Potts model 








2 

q 



Figure 2: Phase diagram for the triangular-lattice Potts model in the (g, f )-plane. The upper 
solid curve shows the ferromagnetic {v > 0) phase-transition curve. The dots • (resp. circles 
o) indicate the known second-order (resp. first-order) transition values Vc{q) for integer q. 
The lower solid curve shows the hypothesized new phase-transition curve for small positive 
q, as estimated numerically in this paper (see Section 7.7). The dashed and dot-dashed 
curves represent the middle and lower branches of (1.12), respectively. The hatched region 
corresponds to (a portion of) the conjectured Berker-Kadanoff phase. The horizontal dashed 
line corresponds to infinite temperature (f = 0), and the horizontal dotted line corresponds 
to the zero-temperature antiferromagnet {v = —1). 
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Zeros sq lattice — 2p 



Zeros sq lattice — 3p 




-0.5 0.5 -0.5 0.5 

Re(w) Re(w) 

(c) (d) 



Figure 3: Limiting curves for square-lattice strips of width (a) L = 2, (b) L = 3, (c) L = 4, 
and (d) L = 5 with cyhndrical boundary conditions. We also show the zeros for the strips 
Lp X (5L)f (black □) and Lp x (10L)f (red o) for the same values of L. 
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Zeros sq lattice = 4 




-0.27765 -0.2776 -0.27755 

Re(w) 

Figure 4: Blow-up of Figure 3(c) around the T point Wq^ ~ —0.2776248333. The multiple 
point is located at w ^ -0.2775806860. 
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Zeros sq lattice L^^ = 6p Zeros sq lattice L^^ = 7] 




-0.5 0.5 -0.5 0.5 

Re(w) Re(w) 

(c) (d) 



Figure 5: Limiting curves for square-lattice strips of width (a) L = 6, (b) L = 7, (c) L = 8, 
and (d) L = 9 with cyhndrical boundary conditions. We also show the zeros for the strips 
Lp X (5L)f (black □) and Lp x (10L)f (red o) for the same values of L. 



93 



Zeros sq lattice L = 10 



T 1 1 1 1 1 1 1 r 




I I I I I I I I I I 

-0.5 0.5 

Re(w) 

Figure 6: Limiting curve for the square-lattice strip of width L = 10 with cyhndrical bound- 
ary conditions. We also show the zeros for the strips lOp x SOp (black □) and lOp x lOOp 
(red o). 
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Figure 7: Limiting curves for square-lattice strips of widths L: 2 (black), 3 (red), 4 (green), 
5 (blue), 6 (pink), 7 (brown), 8 (red), 9 (green), and 10 (blue). The dotted-dashed vertical 
(brown) line marks the line Re(w) = —1/2. As L increases, the limiting curve moves towards 
the right. We have also depicted our estimate for the L —>■ oo limiting curve in the interval 
—0.33 ^ Imw < 0.33 (black dots) and our very rough estimate beyond this (black dotted- 
dashed curve). 
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Figure 8: Limiting curves for triangular-lattice strips of width (a) L = 2, (b) L = 3, (c) 
L = 4, and (d) L = 5 with cylindrical boundary conditions. We also show the zeros for the 
strips Lp X (5L)f (black □) and Lp x (10L)f (red o) for the same values of L. 



96 





0.4 



Re(w) Re(w) 

(c) (d) 

Figure 9: Limiting curves for triangular-lattice strips of width (a) L = 6, (b) L = 7, (c) 
L = 8, and (d) L = 9 with cylindrical boundary conditions. We also show the zeros for the 
strips Lp X (5L)f (black □) and Lp x (10L)f (red o) for the same values of L. 
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Figure 10: Limiting curve for the triangular- lattice strip of width L = 10 with cylindrical 
boundary conditions. We also show the zeros for the strips lOp x 50f (black □) and lOp x lOOp 
(red o). 
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Figure 11: Limiting curves for triangular-lattice strips of widths L: 2 (black), 3 (red), 4 
(green), 5 (blue), 6 (pink), 7 (brown), 8 (red), and 9 (green). As L increases, the limiting 
curve moves towards the right. We have also depicted our estimate for the L ^ oo limiting 
curve in the interval —0.23 < Imw < 0.23 (black dots) and our very rough estimate beyond 
this (black dotted-dashed curve). 
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Figure 12: Real part of the free energy for square-lattice spanning forests as a function of 
w, for strips of width L = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 (A), 9 (0), 10 (x), 
11 (+), 12 (*), 13 (©) and 14 (©). To make clearer any even-odd effect we have displayed 
in red (resp. blue) the points corresponding to even (resp. odd) L. The black solid curve is 
obtained by extrapolating the finite-width data to L — oo and then joining the points. The 
violet dashed line is the Fade [20, 20] approximant to our longest small-w series. Finally, the 
pink dot-dashed curve corresponds to the \axge-w expansion (3.20) through order w~^. 
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Figure 13: First derivative fL{w) of the square-lattice free energy for strips of width L = 2 
(•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 (A), 9 {<}) and 10 (x). Points with even (resp. odd) L 
are shown in red (resp. blue). The violet dashed curve on the right corresponds to the Fade 
approximant [20, 20] to our longest small-w series. The pink dot-dashed curve corresponds 
to the derivative of the large-w expansion (3.20), through order w^^. The vertical brown 
dot-dot-dashed line marks the point Wq = —1/4. 
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Figure 14: First derivative fi{w) of the square-lattice free energy for strips of widths 2 < 
L < 10 close to ifo(sq) = — 1/4. (a) Curves for even (resp. odd) L are shown in red (resp. 
blue). The violet dashed curve on the right corresponds to the Fade approximant [20,20] 
to our longest small-w series. The pink dot-dashed curve corresponds to the derivative of 
the large-w expansion (3.20), through order The vertical brown dot-dot-dashed line 

marks the point wq = —1/4. The vertical brown dot-dot-dashed lines near the bottom of 
the figure mark the position of the points wo-{L) for L = 6,8, 10. (b) Curves for strips of 
widths L = 2, 4, 6, 8, 10 are plotted versus w — wo-{L). 
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Figure 14: First derivative fi{w) of the square-lattice free energy for strips of widths 2 < 
L < 10 close to wo(sq) = — 1/4. (a) Curves for even (resp. odd) L are shown in red (resp. 
blue). The violet dashed curve on the right corresponds to the Fade approximant [20,20] 
to our longest small-iu series. The pink dot-dashed curve corresponds to the derivative of 
the large-w expansion (3.20), through order w~'^. The vertical brown dot-dot-dashed line 
marks the point wq = —1/4. The vertical brown dot-dot-dashed lines near the bottom of 
the figure mark the position of the points wo~.{L) for L = 6,8, 10. (b) Curves for strips of 
widths L = 2, 4, 6, 8, 10 are plotted versus w — wq-{L). 
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Figure 15: Second derivative f'liw) of the square- lattice free energy for strips of width L = 2 
(•), 3 (■), 4 (A), L = 5 (♦), 6 (o), L = 7 (□), 8 (A), L = 9 {<>), and 10 (x). Points 
with even (resp. odd) L are shown in red (resp. blue). The violet dashed curve on the right 
corresponds to the Pade approximant [20, 20] to our longest small-w series. The pink dot- 
dashed curve corresponds to the second derivative of the large-w expansion (3.20), through 
order w~^. The vertical dot-dot-brown dashed line marks the point wq = —1/4. 
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Figure 16: Second derivative f'liw) of the square-lattice free energy for strips of widths 
2 < L < 10 close to wo(sq) = —1/4. (a) Curves for even (resp. odd) L are shown in red 
(resp. blue). The violet dashed curve on the right corresponds to the Fade approximant 
[20, 20] to our longest small-w series. The pink dot-dashed curve corresponds to the second 
derivative of the large-iy expansion (3.20), through order w~^. The vertical brown dot-dot- 
dashed line marks the point Wq = —1/4. The vertical brown dot-dot-dashed lines near the 
top of the figure mark the position of the points Wq-{L) for L = 6,8,10. (b) Curves for 
strips of widths L = 2, 4, 6, 8, 10 are plotted versus w — Wq_{L). 



105 




-0.01 -5x10"^ 

w - Wo_(L) 

(b) 

Figure 16: Second derivative f'l{w) of the square-lattice free energy for strips of widths 
2 < L < 10 close to lyo(sq) = —1/4. (a) Curves for even (resp. odd) L are shown in red 
(resp. blue). The violet dashed curve on the right corresponds to the Fade approximant 
[20, 20] to our longest small-w series. The pink dot-dashed curve corresponds to the second 
derivative of the large-w expansion (3.20), through order w"^. The vertical brown dot-dot- 
dashed line marks the point Wq = —1/4. The vertical brown dot-dot-dashed lines near the 
top of the figure mark the position of the points Wo^{L) for L = 6,8,10. (b) Curves for 
strips of widths L = 2, 4, 6, 8, 10 are plotted versus w — Wo_{L). 
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Figure 17: Estimates for the square-lattice central charge c{w) obtained by fitting the free 
energy to the Ansatze (a,b) Re fiiw) = Ref{w) + [c{w)tt / 6]L~'^ , and (c,d) Re fiiw) = 
Ref{w) + [c{w)tt/6]L-^ + AL-^. Fits are performed for L,ni„ = 2 (•), 3 (■), 4 (A), 5 (♦), 
6 (o), 7 (□), 8 (A), 9 (0), 10 (x), 11 (+), 12 (*) and 13 (©). Points with even (resp. odd) 
L are shown in red (resp. blue). The black dot at 1/ty = 0, c = —2 marks the theoretical 
prediction. The vertical brown dot-dot-dashed line marks the point wq = —1/4. 
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Figure 18: Estimates for the square-lattice central charge c{w) obtained by fitting the free en- 
ergy to the Ansatz Re = Re f{w) + [c{w)7T/6]L-^+Alog\og L/{L'^ log L)+B/{L'^ logL). 
Fits are performed for L^^i^ = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 (A), 9 {<», 10 (x), and 
11 (-I-). Points with even (resp. odd) L are shown in red (resp. blue). For w = —0.3, to avoid 
parity effects, we have performed the fits using data with L = L^in, -^min+2, Linin+4, 



108 



sq-lattice spanning forests 



C\2 - 



T — I — I — I — I — I — r 



1 — I — I — I — I — I — I — I — I — I — r 





oooooooooooooooo 
" □ n □ □ □ ° ° ° ° ° ° ° ' ° □ 



xxxxxxx^ 



■0 

f 7 ■ 



1 



♦ 

o 





o 


♦ 








□ 






♦ 


■ □ 


• A 


o 






10 


OX 


□ 





o 






A 



□ 


□ 


1 • 




X 


A 



A 




• 










♦♦♦♦♦♦♦♦♦♦♦♦*- 

ooooooooooo 
^aaaaaaaaaa 

Ooooooooooo 



J L 



_L!_L 



-2 



-1 





w 

(a) 



Figure 19: Values of the square-lattice inverse correlation length {w) = log|A^/Aj| for 
(a) j = 1 and (b) j = 2. Symbols indicate strip widths L = 2 (•), 3 (■), 4 (A), 5 (♦), 
6 (o), 7 (□), 8 (A), 9 (<C>), and 10 (x). Points with even (resp. odd) L are shown in red 
(resp. blue). The black solid circles (•) correspond to the extrapolated infinite- volume limit 
of the finite-size data (see text). The vertical dot-dot-brown dashed line marks the point 
Wo = -1/4. 
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Figure 19: Values of the square-lattice inverse correlation length $,^^{w) = log|A^/Aj| for 
(a) j = 1 and (b) j = 2. Symbols indicate strip widths L = 2 (•), 3 (■), 4 (A), 5 (♦), 
6 (o), 7 (□), 8 (A), 9 (<C>), and 10 (x). Points with even (resp. odd) L are shown in red 
(resp. blue). The black solid circles (•) correspond to the extrapolated infinite- volume limit 
of the finite-size data (see text). The vertical brown dot-dot-dashed line marks the point 
Wo = -1/4. 
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Figure 20: Estimates for the square-lattice scaling dimension XTj{w) for (a) j = 1 and (b) 
j = 2, obtained by fitting the inverse correlation length to the Ansatz ^^^^(w) = ^^^(^(w) + 
27iXTj{w)L~^. In the region w < — 1/4 we have fixed = 0; in the region w > —l/A we 
have left it variable. Fits are performed for Lmin = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 
(A) and 9 (<)). Points with even (resp. odd) L are shown in red (resp. blue). The vertical 
brown dot-dot-dashed line marks the point wo = —1/4. 
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Figure 20: Estimates for the square-lattice scaling dimension XTj{w) for (a) j = 1 and (b) 
j = 2, obtained by fitting the inverse correlation length to the Ansatz ^^^^(w) = ^^^(^(w) + 
27iXTj{w)L~^. In the region w < — 1/4 we have fixed = 0; in the region w > —l/A we 
have left it variable. Fits are performed for Lmin = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 
(A) and 9 (<)). Points with even (resp. odd) L are shown in red (resp. blue). The vertical 
brown dot-dot-dashed line marks the point wo = —1/4. 
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Figure 21: Real part of the free energy for triangular-lattice spanning forests as a function 
of w, for strips of width L = 2 (•), 3 (■), L = 4 (A), L = 5 (♦), 6 (o), L = 7 (□), 
L = 8 (A), 9 (0), L= 10 (x), L = 11 (+), 12 (*), and L = 13 (©). To make clearer any 
possible even-odd effect we have displayed in red (resp. blue) the points corresponding to 
even (resp. odd) values of L. The black solid curve is obtained by extrapolating the finite- 
width data to L — s> oo and then joining the points. The violet dashed curve is the Fade 
[10, 10] approximant to our longest small-w series. The pink dot-dashed curve corresponds 
to the large-w expansion (3.20) through order w~^. 
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Figure 22: First derivative fiiw) of the triangular- lattice free energy for strips of width L = 2 
(•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 (A) 9 (^), and 10 (x). Points with even (resp. odd) L 
are shown in red (resp. blue). The violet dashed curve on the right corresponds to the Fade 
approximant [10, 10] to our longest small-w series. The pink dot-dashed curve corresponds 
to the derivative of the large-if expansion (3.20), through order w~'^. The vertical brown 
dot-dot-dashed line marks the point wq ~ —0.1753. 
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Figure 23: First derivative fi{w) of the triangular-lattice free energy for strips of widths 
2 < L < 10 close to Wo(sq) ~ —0.1753. Curves for even (resp. odd) L are shown in red 
(resp. blue). The violet dashed solid curve on the right corresponds to the Fade approximant 
[10, 10] to our longest small-if series (this curve is barely visible, as it nearly coincides with 
those for 2 < L < 5). The pink dot-dashed curve corresponds to the derivative of the large- 
ly expansion (3.20), through order W7~^. The vertical brown dot-dot-dashed line marks the 
point Wo ~ —0.1753. 
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Figure 24: First derivative fiiw) of the triangular-lattice free energy for strips of widths 
L = 4,6, 8, 10. (a) Regime w < Wq plotted versus w — Wq_{L). (b) Regime w > Wq plotted 
versus w — Wo+{L). 
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Figure 24: First derivative fiiw) of the triangular-lattice free energy for strips of widths 
L = 4,6, 8, 10. (a) Regime w < Wq plotted versus w — Wq_{L). (b) Regime w > Wq plotted 
versus w — Wo+{L). 
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Figure 25: Second derivative f'l{w) of the triangular-lattice free energy for strips of width 
L = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 (□), 8 (A), 9 (0), and 10 (x). Points with even 
(resp. odd) L are shown in red (resp. blue). The violet dashed curve on the right corresponds 
to the Pade approximant [10, 10] to our longest small-w series. The pink dot-dashed curve 
corresponds to the second derivative of the large-w expansion (3.20), through order w~^. 
The vertical brown dot-dot-dashed line marks the point Wq ~ —0.1753. 
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Figure 26: Second derivative f'l{w) of the triangular-lattice free energy for strips of widths 
2 < L < 10 close to wo(sq) ~ —0.1753. Curves for even (resp. odd) L are shown in red 
(resp. blue). The violet dashed solid curve on the right corresponds to the Fade approximant 
[10, 10] to our longest small-if series. This curve is barely visible, as it very similar to that 
for L = 5. The pink dot-dashed curve corresponds to the second derivative of the large-iy 
expansion (3.20), through order w~^. The vertical brown dot-dot-dashed line marks the 
point Wq ~ —0.1753. 
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Figure 27: Second derivative f'l{w) of the triangular-lattice free energy for strips of widths 
L = 4, 6, 8, 10. (a) Regime w < Wq plotted versus w — Wq_{L). (b) Regime w > Wq plotted 
versus w — wo+{L). 
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Figure 27: Second derivative f'l{w) of the triangular-lattice free energy for strips of widths 
L = 4, 6, 8, 10. (a) Regime w < Wq plotted versus w — Wq_{L). (b) Regime w > Wq plotted 
versus w — wo+{L). 
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Figure 28: Estimates for the triangular-lattice central charge c{w) obtained by fitting the free 
energy to the Ansatze (a,b) Re fL{w) = Ref{w) + [c{w)\/3n/12]L^'^, and (c,f) Re fL{w) = 
Ref{w) + [c(w)\/37r/12]L-2 + A/L^. Fits are performed for L^m = 2 (•), 3 (■), 4 (A), 5 
(♦), 6 (o), 7 (□), 8 (A), 9 {0), 10 (x) 11 (+), and 12 (*). Points with even (resp. odd) 
L are shown in red (resp. blue). The black dot a.t 1/w = 0, c = —2 marks the theoretical 
prediction. The vertical brown dot-dot-dashed line marks the point wq ~ —0.1753. 
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Figure 29: Estimates for the triangular-lattice central charge c{w) obtained by fitting the 
free energy to the Ansatz Re fL{w) = Ref{w) + [c(w)7r/6]L"^ + A loglog L/(L^ log L) + 
B/{L'^\ogL). Fits are performed for L^,^ = 2 (•), 3 (■), 4 (a), 5 (♦), 6 (o), 7 (□), 8 (A), 
9 (<^), and 10 (x). Points with even (resp. odd) L are shown in red (resp. blue). 
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Figure 30: Values of the triangular-lattice inverse correlation length C^J-^^w) = log |A^/Aj| for 
(a) j = I and (b) j = 2. Symbols indicate strip widths L = 2 (•), 3 (■), 4 (A), 5 (♦), 
6 (o), 7 (□), 8 (A), 9 (<)), and 10 (x). Points with even (resp. odd) L are shown in red 
(resp. blue). The black solid circles (•) correspond to the extrapolated infinite- volume limit 
of the finite-size data (see text). The vertical dot-dot-brown dashed line marks the point 
Wo ^ -0.1753. 
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Figure 30: Values of the triangular-lattice inverse correlation length C^J-^^w) = log |A^/Aj| for 
(a) j = I and (b) j = 2. Symbols indicate strip widths L = 2 (•), 3 (■), 4 (A), 5 (♦), 
6 (o), 7 (□), 8 (A), 9 (<)), and 10 (x). Points with even (resp. odd) L are shown in red 
(resp. blue). The black solid circles (•) correspond to the extrapolated infinite- volume limit 
of the finite-size data (see text). The vertical brown dot-dot-dashed line marks the point 
Wo ^ -0.1753. 



125 



CO 



tri-lattice spanning forests 

T — I — r-"-! — I — I — I — I — r — I — I — I — I — I — I — I — I — I — r 



♦ O 



C\2 



E-i 
XI 



nil??!??!?-' 



J I I I I I I I L I I I I I I I I I L 



-2 



-1 







w 

(a) 

Figure 31: Estimates for the triangular-lattice scaling dimension XTj{w) for (a) j = 1 and (b) 
j = 2, obtained by fitting the inverse correlation length to the Ansatz ^J^{w) = ^j^^^iw) + 
27iXTj{w)L~^ . In the region w < —0.175 we have fixed ^^^^ = 0; in the region w > —0.175 
we have left it variable. Fits are performed for Lmin = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 
(□), 8 (A), and 9 (<)). Points with even (resp. odd) L are shown in red (resp. blue). The 
vertical brown dot-dot-dashed line marks the point Wq ~ —0.1753. 
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Figure 31: Estimates for the triangular-lattice scaling dimension XTj{w) for (a) j = 1 and (b) 
j = 2, obtained by fitting the inverse correlation length to the Ansatz 
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Tj\w)L~^ . In the region w < —0.175 we have fixed ^^^^ = 0; in the region w > —0.175 
we have left it variable. Fits are performed for Lmin = 2 (•), 3 (■), 4 (A), 5 (♦), 6 (o), 7 
(□), 8 (A), and 9 (<)). Points with even (resp. odd) L are shown in red (resp. blue). The 
vertical brown dot-dot-dashed line marks the point Wq ~ —0.1753. 
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Figure 32: Finite-L limiting curves and phase diagram (Figure 7) for the square lattice, 
mapped to the 1/w plane. Our best estimate for the infinite- volume phase boundary Boo is 
depicted in black dots (for —0.33 ^ Imw < 0.33) and as a black dotted-dashed curve (very 
rough estimate). 
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Figure 33: Finite-L limiting curves and phase diagram (Figure 11) for the triangular lattice, 
mapped to the 1/w plane. Our best estimate for the infinite- volume phase boundary Boo is 
depicted in black dots (for —0.23 ^ Imw < 0.23) and as a black dotted-dashed curve (very 
rough estimate). 
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Figure 34: Results of studying the small-if series expansions for the square lattice using 
the ratio method, (a) We show the ratio of two consecutive coefficients [cf. (A. 3)] as a 
function of 1/k for the free energy f{w) (o), its first derivative f'{w) (□), and its second 
derivative f"{w) (A), (b) We show the biased estimate (based on wq = —1/4) for the critical 
parameter A for the same three functions as above. 
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Figure 35: Results of studying the small-w series expansions for the square lattice using the 
differential-approximant method. We show the estimate of the critical exponent A versus 
the location of the dominant singularity Wq for (a) the spanning- forest free energy f{w), (b) 
its first derivative f'{w), and (c) its second derivative f"{w). The vertical dashed line marks 
the theoretical prediction wo(sq) = —1/4. The results for first-order (resp. second-order) 
approximants K = 1 (resp. K = 2) are denoted with red squares □ (resp. black circles o). 
In (d) we show the results for the second derivative /" of the triangular-lattice free energy. 
In this case the vertical dashed line marks our best estimate for wo{tn) ~ —0.1753. 
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